AD-A222  995 


one  r.LE  CORY 


/^C-0  3-0/. <5 


P\ 


JC-2053-A 


STUDIES  OF  ECCM  IMPROVEMENTS 
FOR  FREQUENCY-HOPPING 
CPFSK  SYSTEMS 


DTIC 

SELECTE  M 
JUN  1 5  1990 1  1 


FINAL  REPORT 
May  30,  1990 


PREPARED  FOR 


U.  S.  ARMY  RESEARCH  OFFICE 


CONTRACT 

DAAL03-89-C-0010 


The  views,  opinions,  and/or  findings  contained  in  this  report  are  those  of  the 
authors  and  should  not  be  construed  as  an  official  Department  of  the  Army 
position,  policy,  or  decision  unless  so  designated  by  other  documentation. 


jcrriUBUTIOW  STATEMENT  X 

Approved  lot  public  release^ 
Distnbunoa  Unlimited 


l 


J.  8.  LEE  ASSOCIATES,  INC. 

2001  4EFFER80N  DAVI8  HIGHWAY,  8UITE  601 
ARLINGTON,  VIRGINIA  22202 


Approved  for  public  release;  dlotrlbutlon  unlimited. 


80  06  . 14  176 


t«.  report  security  classification 


nnroi 


2a.  SECURITY  CLASSIFICATION  AUTHORITY 


2b.  DECLASSIFICATION  /DOWNGRADING  SCHEDUL 


4.  PERFORMING  ORGANIZATION  REPORT  NUMBER(S) 

JC-2053-A 


REPORT  DOCUMENTATION  PAGE 


lb.  RESTRICTIVE  MARKING 


3.  DISTRIBUTION  /  AVAILABILITY  OF  REPORT 

Approved  for  public  release; 
distribution  unlimited. 


5  MONITORING  ORGANIZATION  REPORT  NUMBER(S) 


6a.  NAME  OF  PERFORMING  ORGANIZATION 

J.  S.  Lee  Associates,  Inc. 


6c  AODRESS  (Gty,  Start,  and  ZIP  Code) 

2001  Jefferson  Davis  Highway,  Suite  601 
Arlington,  Virginia  22202 


A/Ca  3o/  ■  l  _ 


6b  OFFICE  SYMBOL  7a.  NAME  OF  MONITORING  ORGANIZATION 
(If  applicable) 

U.  S.  Army  Research  Office 


7b.  ADDRESS  (City,  Start,  and  ZIP  Cock) 

P.  0.  Box  12211 

Research  Triangle  Park,  NC  27709-2211 


8a.  NAME  OF  FUNDING /SPONSORING 
ORGANIZATION 

U.  S.  Army  Research  Office 


8c  ADDRESS  (City,  State,  and  ZIP  Coda) 

P.  0.  Box  12211 

Research  Triangle  Park,  NC  27709-2211 


8b.  OFFICE  SYMBOL  9.  PROCUREMENT  INSTRUMENT  IDENTIFICATION  NUMBER 
(If  appUcabk) 

DAAL03-89-C-0010 


10.  SOURCE  OF  FUNDING  NUMBERS 


PROGRAM 
ELEMENT  NO 


PROJECT 

NO 


WORK  UNIT 
ACCESSION  NO 


1 1 .  TITLE  (Include  Security  Classification) 

Studies  of  ECCM  Improvements  for  Frequency-Hopping  CPFSK  Systems 


12  PERSONAL  AUTHOR®  . .  ,  „  „  . 

J.  S.  Lee,  L.  E.  Miller,  R.  H.  French 


sura 


14  D^TE^F  REJgIRT  (Tear.  Month.  Day) 


5  PAGE  COUNT. 

122  +  vvn 


13a.  TYPE  OF  REPORT 

Final 


16.  SUPPLEMENTARY  NOTATION  .  .  ,  .  .  .  . 

The  view,  opinions  and/or  findings  contained  in  this  report  are  those 

of,£he  authqr(s) .and  should  notJbejconst5ugd  as^an  gfficial  Department  of  the  Army  position, 


COSATI  COOES 


GROUP  SUB-GROUP 


IB.  SUBJECT  TERMS  ( Continue  on  reverse  if  necessary  and  identify  by  .block  number) 

Frequency  hopping,  digital  FM  radio,  ECCM,  jamming, 
statistical  communication  theory,  phase  distributions, 
coding,  synthesizers,  diversity  combining 


19.  ABSTRACT  ( Continue  on  reverse  if  necessary  and  identify  by  block  number) 

y  The  performance  of  frequency-hopping  CPFSK  modulations  (such  as  used  in  several  Army 
radio  systems)  in  the  presence  of  worst-case  partial-band  noise  jamming  (WCPBNJ)  is 
derived  for  both  limiter-discriminator  and  differential  detector  forms  of  receiver.  The 
transmitter  is  assumed  to  repeat  the  data  L  times  on  different  hops,  and  the  receiver 
combines  the  L  chips  for  each  data  symbol  in  a  way  to  mitigate  the  effects  of  the  jamming. 
Hard-decision  combining  simply  makes  a  binary  decision  on  each  chip,  then  a  majority 
decision  on  the  whole  symbol;  this  limits  a  jammed  chip  to  one  vote.  AGC  combining 
weights  each  analog  chip  sample  in  inverse  proportion  to  the  received  noise  power  and 
decides  on  the  basis  of  the  sum  of  weighted  chip  samples;  jammed  chips  tend  to  be  weighed 
less.  Analysis  and  computations  show  that  either  form  of  ECCM  combining  can  produce  a 
27  dB  gain  against  WCPBNJ,  and  with  error-control  coding,  a  31-33  dB  gain.  Coding  alone 
can  give  a  28-29  dB  gain.  Soft  decoding  of  hard-decision  diversity  combining  gives  1.4  dB 
more  gain  than  hard  decoding. (continued) 


20.  DISTRIBUTION  /AVAILABILITY  OF  ABSTRACT 
□  UNCLASSIFIED/UNLIMITED  9  SAME  AS  RPT. 


22a.  NAME  OF  RESPONSIBLE  INDIVIDUAL 

Dr.  James  W.  Gault 


□  OTIC  USERS 


OO  FORM  1473, 84  mar 


21.  ABSTRACT  SECURITY  CLASSIFICATION 
Unclassified 


■m.w: 


83  APR  edition  may  ba  used  until  exhausted 
All  other  editions  are  obsolete 


UNCLASSIFIED 


i 


UNCLASSIFIED _ 

SECURITY  CLASSIFICATION  OP  THII 


19.  ABSTRACT  (continued) 

It  is  also  shown  that,  due  to  the  spectral  attenuation 
characteristics  of  the  receiver's  IF  filter,  the  effects  of  the  transient  phase 
difference  between  the  incoming  hopping  waveform  and  the  receiver's  synthesizer 
are  primarily  confined  to  an  interval  of  duration  T  +  ts  +  tau,  where  T  is  the 
symbol  duration  (inverse  of  IF  bandwidth),  ts  is  the  synthesizer  switching  time, 
and  tau  is  the  offset  in  synchronization  between  the  received  signal  and  the 
receiver's  synthesizer.  ..Therefore  a  one-symbol  guard  time  is  sufficient  to  avoid 
effects  of  the  transient  on  the  data;  normally,  a  considerably  larger  guard  time 
is  used  for  the  purpose  of  emission  control.  Therefore,  emission  control  has  more 
of  an  influence  on  the  amount  of  data  per  hop  than  do  transients. 


UNCLASSIFIED 


SECURITY  CLASSIFICATION  OF  THIS  FAGE 


JC-2053-A 


STUDIES  OF  ECCM  IMPROVEMENTS 
FOR  FREQUENCY-HOPPING 
CPFSK  SYSTEMS 


FINAL  REPORT 
May  30,  1990 


PREPARED  FOR 

U.  S.  ARMY  RESEARCH  OFFICE 

CONTRACT 


|  Accesioi.  Ft  .  1 

_ _ _  _  _ _ 

_  _  V 

NTiS 

oi  vu  y 

one 

T"-U  □ 

OU-.T-J  □ 

JuSt1): 

By 

District, o  i  /  1 

A 

vjci.iti.My  Cl’Gcs 

or 

LJiSt 

Ou’L'C.ji 

1 

DAAL03-89-C-0010 


The  views,  opinions,  and/or  findings  contained  in  this  report  are  those  of  the 
authors  and  should  not  be  construed  as  an  official  Department  of  the  Army 
position,  policy,  or  decision  unless  so  designated  by  other  documentation. 


J.  S.  LEE  ASSOCIATES,  INC. 

2001  JEFFERSON  DAVIS  HIGHWAY,  SUITE  601 
ARLINGTON,  VIRGINIA  22202 


Contents 


TABLE  OF  CONTENTS 


Report  documentation  page  and  abstract  .  i 

Lists  of  Figures  and  Tables  .  vii 


1.0  INTRODUCTION 


1 


1.1  CURRENT  WORK  RELATED  TO  FH/CPFSK  SYSTEMS  . 

1.2  RATIONALE  FOR  FURTHER  STUDIES  OF  ECCM  COMBINING  . 

1.3  FH/CPFSK  RECEIVER  IMPLEMENTATION  ISSUES  . 

1.4  ORGANIZATION  OF  THE  REPORT  . 

1.5  SUMMARY  OF  FINDINGS  . 

2.0  ANALYSIS  MODELS  AND  METHODOLOGY  . 

2.1  MODELS  FOR  THE  SYSTEM  STUDIED  AND  ITS  PARAMETERS  .... 

2.1.1  Encoding  and  Transmission  Schemes  . 

2. 1.1.1  Error  control  coding  . 

2. 1.1. 2  Inter-hop  switching  time  . 

2. 1.1. 3  Transmitted  waveform  and  its  parameters  . 

2.1.2  Reception  Scheme  and  Received  Waveform  . 

2. 1.2.1  Received  signal  waveform  . 

2. 1.2. 2  Noise  and  jamming  . 

2.1.3  Demodulator  Options  . 

2. 1.3.1  Differential  phase  in  the  absence  of  noise  . 

2. 1.3. 2  Probability  density  function  for  the  limiter-discriminator 

receiver  output  . 

2. 1.3. 3  Probability  density  function  for  the  output  of 

a  differential  detector  . 

2.1.4  Generic  Chip  Combining  Schemes  . 

2.2  FORMULATIONS  FOR  THE  PERFORMANCE  MEASURES  . 

2.2.1  Uncoded  Bit  Error  Probability  . 

2.2.2  Coded  Bit  Error  Probability  . 


1 

2 


n 

i 

8 

9 

9 

9 

11 

12 

12 

13 

13 

16 

IS 

IS 

20 

24 

25 

26 
26 
27 


2.3  SUMMARY  OF  CALCULATION  METHODOLOGY  .  29 

2.3.1  Characteristic  Function  Method  for  Obtaining  the 

Bit  Error  Probability  .  29 

2.3.2  Characteristic  Function  with  FM  Clicks  Included  .  29 

2.3.3  Characteristic  Function  with  Adaptive  Weights  .  30 

3.0  FH/CPFSK  IMPLEMENTATION  STUDIES  .  31 

3.1  DISCUSSION  OF  IMPLEMENTATION  ISSUES  .  31 

3.1.1  Factors  Affecting  the  Selection  of  a  Hopping  Rate  .  31 

3.1.2  Emissions  Control  and  Hop  Transient  Effects  .  33 

3.1.3  Synchronization  Considerations  .  34 

3.2  ANALYSIS  OF  HOP  TRANSIENT  EFFECTS  .  34 

3.2.1  Model  for  Synthesizer  Transient  Waveforms  .  35 

3. 2. 1.1  Received  and  local  oscillator  transient  phase  functions  .  35 

3.2. 1.2  Mixer  output  transient  phase  difference  .  36 

3.2.2  Effect  of  Receiver  Filtering  on  Transients  .  37 


v 


Contents 


4.0  INVESTIGATION  OF  ECCM  COMBINING  TECHNIQUES  .  43 

4.1  BACKGROUND:  NONCOHERENT  COMBINING  LOSS .  43 

4.1.1  Diversity  Sum  Error  without  FM  Noise  Clicks  .  43 

4.1.2  Effects  of  FM  Noise  Clicks  .  45 

4.2  HARD-DECISION  COMBINING  RECEIVER  .  46 

4.2.1  Analysis  of  the  HD  Receiver  .  46 

4.2.2  Results  Using  the  HD  Receiver  with  Diversity  .  49 

4.3  LIMITER-DISCRIMINATOR  TECHNIQUES  .  54 

4.3.1  Approximations  to  the  Differential  Phase  PDF  and  CHF  .  54 

4. 3. 1.1  Expansion  of  the  CHF  in  terms  of  cumulants  .  57 

4. 3. 1.2  CHF  based  on  approximation  of  the  PDF  .  58 

4.3.2  Adaptive  Gain  Control  Scheme  .  64 

4.3.2. 1  Analysis  .  6J 

4. 3. 2. 2  Numerical  results  .  66 

4.4  DIFFERENTIAL  DETECTION  RECEIVER  TECHNIQUES  .  75 

4.4.1  Characteristic  Function  for  the  Differential 

Detector  in  Jamming  .  75 

4.4.2  AGC  Weighting  Scheme  .  76 

4.4.2. 1  Analysis  .  76 

4. 4. 2. 2  Numerical  results  .  77 

4.5  RECEIVER  COMPARISONS  WITHOUT  AND  WITH  CODING  .  89 

4.5.1  Receiver  Comparisons  Without  Coding  .  89 

4. 5. 1.1  Comparison  of  AGC  results  .  89 

4. 5. 1.2  Comparison  of  limiter-discriminator  results  .  90 

4. 5. 1.3  General  comparison  without  coding  .  92 

4.5.2  Coded  Performance  Results  and  Comparisons  .  92 

4.5.2. 1  Results  for  hard-decision  diversity  combining 

and  hard  decoding  .  94 

4. 5. 2. 2  Results  for  soft-decision  diversity  combining 

and  hard  decoding  .  96 

4. 5. 2. 3  Results  for  hard-decision  diversity  combining 

and  soft  decoding  .  96 

4. 5. 2. 4  Summary  comparison  of  coded  performance  results  .  99 

Appendix  A:  Intersymbol  Interference  Analysis  . 103 

Appendix  B:  Properties  of  the  Truncated  Gaussian  Distribution  . 112 

Appendix  C:  Functions  Related  to  the  Gaussian  Filter  .  116 

References  . 120 


Contents 


LIST  OF  FIGURES 

Figure  no.  Title  page  do. 

1-1  FH/CPFSK  Performance  in  Partial-Band  Noise  Jamming  .  3 

1-2  FH/CPFSK  System  Model  .  4 

1- 3  FH/CPFSK  Diversity  Sum  Performance  Assuming  Perfect  Side 

Information  and  No  Thermal  Noise  .  6 

2- 1  Transmission  Scheme  .  10 

2-2  Reception  Scheme  .  14 

2-3  Demodulator  Options  .  19 

2- 4  Probability  Density  Function  for  Differential  Phase 

Including  FM  Clicks  .  23 

3- 1  Transient  Phase  Difference  Models  .  38 

3- 2  IF  Signal  Amplitude  and  Phase  Transient  Waveforms  .  42 

4- 1  Unjammed  P,  for  FH/CPFSK  for  L= 1,  2,  3,  4,  6  .  44 

4-2  Effect  of  Fullband  Jamming  on  L  hops/bit  FH/CPFSK  .  47 

4-3  Exact  and  Approximate  CPFSK  Error  Probability  Comparison 

for  Limiter-Discriminator  Detection  .  50 

4-4  Hard-Decision  Diversity  Performance  for  L—  1  to  7  .  51 

4-5  Hard-Decision  Diversity  Performance  for  L—1  to  12  .  52 

4-6  Hard-Decision  Diversity  Performance  for  odd  L  .  53 

4-7  Hard-Decision  Diversity  Performance  for  Eil/.N'o  =  15dB  .  55 

4-8  Hard- Decision  Diversity  Performance  for  Eb/X0  =  13 dB  .  56 

4-9  Comparison  of  Exact  and  Gaussian  Approximation  PDF's  .  59 

4-10  Comparison  of  Exact  and  Analytical  Mixture  Approximation  PDF’s  .  61 

4-11  Comparison  of  Exact  and  Simplified  Mixture  Approximation  PDF's  .  62 

4-12  Comparison  of  Exact  and  Approximate  CPFSK  Bit  Error 

Probability  Calculations  .  63 

4-13  Performance  of  FH/CPFSK  in  Partial-Band  Jamming  for  L= 1  Hop/Bit, 

L/D  Receiver,  and  AGC  Weighting  .  68 

4-14  Performance  of  FH/CPFSK  in  Partial-Band  Jamming  for  L= 2  Hops/Bit. 

L/D  Receiver,  and  AGC  Weighting  .  69 

4-15  Performance  of  FH/CPFSK  in  Partial-Band  Jamming  for  L= 3  Hops/Bit, 

L/D  Receiver,  and  AGC  Weighting  .  71 

4-16  Performance  of  FH/CPFSK  in  Partial-Band  Jamming  for  L= 4  Hops/Bit, 

L/D  Receiver,  and  AGC  Weighting  .  72 

4-17  Performance  of  FH/CPFSK  in  Partial-Band  Jamming  for  L= 6  Hops/Bit. 

L/D  Receiver,  and  AGC  Weighting  .  73 


Vll 


Contents 


LIST  OF  FIGURES  (continued) 

] 

4-18 

Diversity  Performance  of  FH/CPFSK  in  Partial-Band  Noise  Jamming 

with  AGC  Weighting  . 

.  74 

I 

4-19 

Performance  of  FH/CPFSK  in  Partial-Band  Jamming  for  L=  1, 

Differential  Detection,  and  AGC  Weighting  . 

.  78 

I 

4-20 

Performance  of  FH/CPFSK  in  Partial-Band  Jamming  for  1  =  1, 

Differential  Detection,  and  Ej/Xo  — 15dB  . 

.  79 

4-21 

Performance  of  FH/CPFSK  in  Partial-Band  Jamming  for  L=  1, 

Differential  Detection,  and  Et/>fo~15dB  . 

.  80 

4-22 

Performance  of  FH/CPFSK  in  Partial-Band  Jamming  for  L— 2, 

Differential  Detection,  and  AGC  Weighting  . 

.  82 

1 

4-23 

Performance  of  FH/CPFSK  in  Partial-Band  Jamming  for  L= 2, 

Differential  Detection,  and  Ei/J'r0-T5dB  . 

.  83 

J 

4-24 

Performance  of  FH/CPFSK  in  Partial-Band  Jamming  for  L= 3, 

Differential  Detection,  and  AGC  Weighting  . 

.  84 

J 

4-25 

Performance  of  FH/CPFSK  in  Partial-Band  Jamming  for  L— 4, 

Differential  Detection,  and  AGC  Weighting  . 

.  85 

2 

4-26 

Performance  of  FH/CPFSK  in  Partial-Band  Jamming  for  L= 5, 

Differential  Detection,  and  AGC  Weighting  . 

.  86 

1 

4-27 

Performance  of  FH/CPFSK  in  Partial-Band  Jamming  for  L=6, 

Differential  Detection,  and  AGC  Weighting  . 

.  87 

J 

4-28 

Diversity  Performance  of  FH/CPFSK  in  Partial-Band  Jamming: 

Differential  Detection,  AGC  Weighting  . 

.  SS 

J 

4-29 

Comparison  of  the  Diversity  Performances  of  FH/CPFSK 

in  Partial-Band  Jamming  . 

.  91 

2 

4-30 

Coded  Performance  with  Hard-Decision  Combining  and  Hard  Decoding 

.  95 

■ 

4-31 

Coded  Performance  with  AGC  Soft-Decision  Combining 

and  Hard  Decoding  . 

.  97 

J 

4-32 

Coded  Performance  wdth  Hard-Decision  Combining  and  Soft  Decoding  .. 

.  98 

1 

4-33 

Coded  Performance  Comparison  with  Optimum  Diversity  . 

.  100 

A-l 

Waveforms  Associated  with  the  Pattern  010  . 

.  107 

A-2 

Waveforms  Associated  with  the  Pattern  0110  . 

.  109 

n 

A-3 

Example  Calculated  Differential  Phase  Waveforms  . 

.  Ill 

| 

TABLES 

J 

4-1 

Requirements  for  10~5  Pe  in  Worst-Case  Jamming  When  Et/>To=20dB 

.  101 

I 

A-l 

ISI  Pattern  Parameter  Values  . 

.  no 

| 

viii 

1 

Introduction 


1.0  INTRODUCTION 

State-of-the-art  Army  tactical  radios,  such  as  JTIDS  (Joint  Tactical  Information 
Distribution  System),  EPLRS  (Enhanced  Position  Location  Reporting  System),  and 
SINCGARS  (Single  Channel  Ground  and  Airborne  Radio  System),  employ  frequency 
hopping  to  render  their  transmissions  “spread  spectrum”  and  therefore  more  difficult  for 
an  opponent  to  intercept  and  to  degrade  through  intentional  electromagnetic  inter¬ 
ference  (jamming).  The  modulation  used  by  these  communications  systems  to  transmit 
data  while  dwelling  at  a  particular  carrier  frequency  is  a  form  of  binary  frequency  modu¬ 
lation  (FM)  known  as  continuous  phase-shift  modulation  (CPSM)  or  minimum-shift 
keying  (MSK),  which  is  a  special  case  of  continuous-phase  frequency-shift  keying  or 
CPFSK.  The  combination  of  frequency  hopping  and  this  modulation  is  denoted  by 
FH/CPFSK. 

1.1  CURRENT  WORK  RELATED  TO  FH/CPFSK  SYSTEMS 

CPFSK  is  the  modulation  of  choice  for  many  systems,  as  opposed  to  other  modula¬ 
tions  such  as  M-ary  noncoherent  frequency-shift  keying  (MFSK),  because  of  its  efficient 
use  of  the  available  spectrum,  that  is.  for  the  same  data  rate,  the  bandwidth  required 
by  CPFSK  to  support  the  reception  of  the  data  with  acceptable  quality  in  terms  of  bit 
error  probability  is  relat’vely  small.  Also,  in  the  care  of  limiter-discriminator  detection 
of  the  hopped  CPFSK  waveform,  it  has  been  estimated  that  CPFSK  enjoys  a  4dB  ad¬ 
vantage  in  performance  over  MFSK  in  the  presence  of  noise  jamming  [1]. 

In  a  previous  report  submitted  to  the  Army  Research  Office  [8],  the  authors 
studied  the  performance  of  FH/CPFSK  when  subjected  to  worst-case  partial-band  noise 
jamming.  Partial-band  noise  jamming,  as  the  name  implies,  features  the  placement  of 
the  jammer’s  effective  available  power,  J,  in  a  fraction  of  the  bandwidth,  W,  over  which 
the  signal  is  hopping.  Rather  than  spreading  the  jamming  noise  power  over  the  entire 
hopping  bandwidth,  with  a  noise  spectral  power  density  of 

Jf,  =  J/W,  (1-1) 

by  concentrating  the  power  in  a  smaller  bandwidth  7W,  where  7  <  1,  the  spectral  power 
density  in  the  jammed  portions  of  the  band  is  Kj/ 7,  a  higher  value  which  can  be 
expected  to  be  more  effective  in  degrading  the  signal.  However,  now  the  signal  is  not 
always  jammed,  so  that  there  is  a  tradeoff  involving  the  choice  of  7  which  results  in  the 
maximum  (worst-case)  value  of  the  average  bit  error  probability  at  the  receiver. 
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Worst-case  partial-band  noise  jamming  can  be  very  effective  in  degrading  the 
FH/CPFSK  system's  performance,  as  illustrated  in  Figure  1-1,  which  is  excerpted  from 
[8]  for  the  case  of  a  receiver  using  differential  detection  to  demodulate  the  data.  The 
figure  shows  the  bit  error  probability  (Pe),  as  a  function  of  the  ratio  of  signal  bit  energy 
to  average  jamming  noise  spectral  power  density,  E J  H j,  assuming  that  the  signal 
energy- to-noise  density  ratio  prior  to  the  onset  of  jamming  is  Et/.tf0  =  15dB  (which  gives 
a  P*<10-5),  and  parametric  in  different  values  of  ■>,  the  jamming  fraction  It  is  clear 
from  the  figure  that  if  the  worst-case  value  of  i  is  always  chosen,  the  dependence  of  Pf 
on  Et/K j ,  rather  than  being  exponential  (as  is  approximately  the  case  for  wideband 
noise),  is  “inverse  linear,”  that  is, 

Pf  oc  (Eb/Kj)~\  worst-case  •>  used.  (1-2) 

Note  that  the  concept  of  worst-case  jamming  through  selection  of  the  jamming  fraction 
■)  is  a  theoretical  concept:  whether  or  not  it  is  practical  to  implement  such  worst-case 
jamming1 2,  it  is  of  interest  to  calculate  the  worst-case  performance  and  ways  for  im¬ 
proving  it. 

A  general  model  of  an  advanced  slow-hopping  FH/CPFSK  communications  system 
is  shown  in  Figure  1-2.  The  advanced  features,  intended  to  improved  the  system  per¬ 
formance  while  under  jamming,  include  error  control  coding,  interleaving.  L  hops  per 
bit  diversity,  and  electronic  counter-countermeasures  (ECCM)  diversity  combining.  In 
[8],  a  “baseline"  version  of  the  system  was  evaluated.  For  this  baseline  system,  at  the 
receiver  the  L  analog  data  sample  values  for  each  bit  are  summed,  and  the  bit  decision 
is  given  by  the  polarity  of  this  sum.  Our  analysis  and  computations  included  a  rigorous 
treatment  of  background  noise,  intersymbol  interference,  and  FM  noise  “clicks."  The 
evaluation,  summarized  in  [15].  determined  that  linear  combining  or  summing  the 
diversity  components  (chips)  of  each  bit  does  not  improve  system  performance  in  worst- 
case  partial-band  noise  jamming.  This  conclusion  holds  for  both  limiter-discriminator 
and  differential  detection  demodulation  techniques. 

1.2  RATIONALE  FOR  FURTHER  STUDIES  OF  ECCM  COMBINING 

The  diversity  sum  method  for  combining  the  L  chips  for  a  given  FH/CPFSK  bit 
can  be  regarded  as  implementing  a  form  of  “soft-decision  metric.”  Note  that  the 
summing  operation  does  not  utilize  any  “side  information”  on  whether  a  particular  chip 
has  been  jammed,  or  how  strongly  it  was  jammed.  Since  there  is  no  mechanism  for 
excluding  or  otherwise  treating  jammed  chips  differently,  when  one  or  more  chips  is 

1See,  for  example,  [20]  and  [29]  for  discussions  of  jammer  implementation. 
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subject  to  jamming,  the  entire  sum  is  corrupted;  this  accounts  for  the  failure  of  the  sum 
metric  to  provide  a  diversity  improvement. 

Now,  for  hopped  MFSK  systems,  it  has  been  shown  [21]  that  if,  instead  of  soft- 
decision  combining  of  the  chips,  we  combine  hard  decisions  on  each  chip,  the  result  is 
that  for  high  Eb/K0  (15dB  or  more)  there  is  a  diversity  gain,  in  the  sense  that  for  a 
particular  J/S  ratio,  more  than  one  hop  per  bit  may  yield  a  lower  bit  error  probability. 
The  improvement  is  due  to  the  hard  decision’s  limiting  a  jammed  chip  to  “one  vote”  in 
the  sum  of  hard  decisions.  Since  no  side  information  is  required,  this  simple  ECCM 
scheme  is  very  attractive  to  consider  for  FH/CPFSK  also,  except  for  the  fact  that 
noncoherent  combining  losses  are  high  for  hard-decision  metrics,  relative  to  soft 
decisions,  in  general. 

In  [8]  it  was  also  shown  that  a  diversity  gain  would  result  from  using  a  “perfect  side 
information”  soft-decision  diversity  combining  scheme  for  FH/CPFSK  which  includes 
only  unjammed  chips  in  the  diversity  sum,  unless  all  chips  happen  to  have  been 
jammed.  Figure  1-3  illustrates  the  theoretical  performance  of  FH/CPFSK  using  such  a 
scheme  in  the  absence  of  any  thermal  noise  (E4/>f0-*cx).  The  error  probability  under 
these  assumptions  is  such  that  for  higher  values  of  Ej/J^j, 

Pe  a  (EJJTj)-1;  (1-3) 

in  the  figure  it  is  evident  that  the  different  curves  for  different  values  of  L  are  related 
such  that  there  is  an  “optimum  diversity”  which  is  a  function  of  E  Jtfj. 

The  implication  of  the  theoretical  results  shown  in  Figure  1-3  is  that  it  is  likely 
that  some  practical  method  exists  for  combining  the  analog  chip  samples  or  soft 
decisions  into  a  metric  that  will  result  in  an  FH/CPFSK  performance  in  worst-case 
partial-band  noise  jamming  that  is  better  than  the  one  obtained  using  the  sum  of  hard 
decisions.  Such  metrics  have  been  found  for  noncoherent  FH/BFSK,  including  a  “self- 
normalizing”  technique  that  does  not  require  side  information  [22,  23]. 

Based  on  all  these  previous  findings,  in  the  work  summarized  in  this  report  we  have 
developed  and  evaluated  soft-decision  schemes  for  combining  FH/CPFSK  diversity 
transmissions  which  produce  a  diversity  improvement  against  worst-case  partial-band 
noise  jamming.  An  analysis  by  Torrieri  [1]  suggested  that  a  concatenation  of  a  convolu¬ 
tional  “outer  code”  with  hard-decision  diversity  combining  (an  “inner  code”)  would 
yield  an  FH/CPFSK  performance  that  improves  as  L,  the  number  of  hops  per  bit, 
increases,  provided  that  soft  decoding  is  used.  Therefore,  in  this  report  we  also  include 
calculations  of  coded  performance  in  combination  with  diversity  combining. 
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1.3  FH/CPFSK  RECEIVER  IMPLEMENTATION  ISSUES 

It  has  been  maintained  [5]  that  a  loss  in  performance  is  experienced,  in  terms  of 
theoretical  error  probability,  by  using  a  differential  detector  instead  of  a  limiter- 
discriminator  for  reception  of  CPFSK,  either  in  Gaussian  noise  or,  when  hopped,  in 
worst-case  partial-band  noise  jamming.  However,  this  conclusion  was  based  on 
comparing  the  performances  of  the  two  receiver  implementations  (parametric  in  the 
modulation  index  and  in  the  receiver  filter’s  bandwidth-time  product)  on  the  basis  of 
the  same  differential  delay  time,  one  bit  (or  chip)  period.  Ekanayake  [4]  has  shown  that 
the  differential  detector’s  performance  can  be  made  very  competitive  with  that  of  the 
limiter-discriminator  by  using  a  differential  delay  value  less  than  a  bit  period.  In  this 
report,  we  do  not  emphasize  the  improvement  in  differential  detector  performance  that 
can  be  obtained  by  selection  of  the  delay,  but  do  develop  the  performance  of  both 
receiver  implementations  when  hard  or  soft  decision  ECCM  combining  is  used  for  de¬ 
modulating  multiple  hops  per  bit  transmissions. 

Another  implementation  issue  concerns  the  phase  transients  or  discontinuities 
experienced  by  an  FH/CPFSK  receiver  during  and  just  following  the  hop  transitions  to 
different  carrier  frequencies.  If  the  signalling  rate  is  increased  slightly  to  allow  for  a 
dwell  time  between  hops,  there  will  be  a  performance  degradation  due  to  the  loss  of 
energy  per  bit,  assuming  the  same  carrier  power  and  average  bit  rate.  Therefore,  it  is 
important  to  calculate  the  amount  of  dwell  time  needed  in  order  to  preclude  significant 
transient  interference  effects  on  the  initial  bits  or  chips  on  each  hop. 

1.4  ORGANIZATION  OF  THE  REPORT 

The  report  is  organized  as  follows:  following  this  introductory  section,  which  in¬ 
cludes  a  summary  of  findings,  the  modelling  and  methodologies  used  in  the  report  are 
summarized  in  Section  2.  Then,  in  Section  3,  we  summarize  our  studies  of  FH/CPFSK 
implementation  issues,  including  hop  transients. 

In  Section  4,  the  major  part  of  our  work  is  presented:  analysis  and  calculations 
showing  the  uncoded  performance  of  FH/CPFSK  against  worst-case  partial-band  noise 
jamming  when  combining  of  the  L  chips  transmitted  per  bit  on  different  hops  takes 
place  according  to  an  “automatic  gain  control”  scheme. 

In  Section  5,  the  performance  of  the  system  using  a  particular  error-control  code  in 
addition  to  ECCM  diversity  combining  is  evaluated  through  calculations. 
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1.5  SUMMARY  OF  FINDINGS 

The  analysis  and  computations  described  in  this  report  support  the  following  con¬ 
clusions  regarding  frequency-hopped  CPFSK  communications  systems. 

Regarding  FH/CPFSK  implementation  issues  (Section  3),  it  is  shown  that  due  to 
the  spectral  attenuation  characteristics  of  the  receiver’s  IF  filter,  the  effects  of  the 
transient  phase  difference  between  the  incoming  hopping  waveform  and  the  receiver's 
synthesizer  are  primarily  confined  to  the  time  interval  (0,  Tc+t»+  r),  referenced  to  the 
input  waveform,  wdth  t=0  at  the  instant  that  the  incoming  signal  begins  to  switch 
frequencies  and  using 

Te  =  CPFSK  modulation  symbol  duration  ss  (IF  bandwidth)-1 

t3  =  synthesizer  switching  time 

r  =  offset  in  receiver  synthesizer  switching  time. 

Therefore,  a  one-symbol  guard  time  is  advisable  before  resuming  the  transmission  of 
data.  Typically,  the  amplitude  shaping  and  other  emissions  control  procedures  used  by 
hopping  radios  dictate  a  larger  off  time,  so  that  these  factors  determine  what  fraction  of 
the  dwell  time  may  be  used  to  transmit  data,  rather  than  switching  transient  effects. 

Regarding  ECCM  combining  techniques  (Section  4),  for  both  limiter-discriminator 
(LD)  and  differential  detection  (DD)  types  of  FH/CPFSK  receiver,  it  is  shown  that 

•  With  no  coding  or  ECCM  diversity  combining,  the  effect  of  worst-case  partial- 
band  noise  jamming  (WCPBNJ)  is  to  require  the  value  of  the  bit-energv-to-jamming- 
noise-spectral-density  ratio  E Jtfj  to  be  43  or  44 dB,  rather  than  about  12  or  15  dB  in 
fullband  jamming  for  the  LD  and  DD,  respectively. 

•  Using  the  optimum  value  of  L,  the  diversity  or  number  of  hops  per  symbol,  it  is 
possible  to  achieve  a  27dB  gain  against  WCPBNJ,  using  either  hard-decision  combining 
or  AGC  combining.  Relative  to  the  DD,  the  LD  receiver  tends  to  suffer  a  loss  in  per¬ 
formance  when  AGC  combining  is  used,  due  to  the  effects  of  FM  noise  clicks,  but 
otherwise  retains  its  usual  advantage. 

•  Coding  alone  can  produce  a  28  or  29 dB  gain  against  the  WCPBNJ. 

•  Coding  plus  optimum  diversity  can  provide  31-33dB  gain.  Soft-decision  decod¬ 
ing  gives  about  1.4  dB  more  gain  than  hard-decision  decoding. 

Several  of  the  analytical  results  are  new.  These  include  the  use  of  a  truncated 
Gaussian  or  a  Gaussian-uniform  mixture  distribution  to  approximate  the  modulo-2ff 
probability  distribution  of  the  differential  phase,  and  expressions  for  its  characteristic 
function  which  include  the  effects  of  FM  noise  clicks  on  the  distribution. 
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2.0  ANALYSIS  MODELS  AND  METHODOLOGY 

Our  studies  concern  the  jammed  performance  of  hopped  binary  FM  communica¬ 
tions,  particularly  under  the  assumption  of  partial-band  noise  jamming  and  the  use  of 
(time)  diversity  and  various  diversity  combining  techniques  to  mitigate  the  jamming. 
In  this  section,  we  present  the  analysis  models  used  to  formulate  measures  of  system  per¬ 
formance,  and  outline  the  methodology  employed  in  calculating  the  performance 
measures. 

2.1  MODELS  FOR  THE  SYSTEM  STUDIED  AND  ITS  PARAMETERS 

In  the  following,  models  for  the  outputs  of  the  several  system  components  are  dis¬ 
cussed  individually. 

2.1.1  Encoding  and  Transmission  Schemes 

Figure  2-1  gives  a  block  diagram  of  the  transmission  scheme  for  the  system 
studied.  Binary  data,  after  error-control  coding  and  possibly  encryption,  are  to  be  trans¬ 
mitted  using  slow-frequency-hopping  digital  FM,  or  CPFSK.  The  coded  symbols  are  to 
be  repeated  on  L  different  hops  in  order  to  increase  the  likelihood  that  some  of  the 
symbols  are  free  of  jamming.  The  figure  suggests  one  of  many  possible  ways  to  accom¬ 
plish  this  objective.  According  to  the  version  shown  in  the  figure,  the  coded  symbols 
are  first  read  into  a  Q-bit  shift  register  (Q-symbol  buffer),  where  Q  is  the  number  of 
symbols  that  can  be  transmitted  in  one  hop  period.  For  example,  if  the  channels 
allotted  to  the  system  support  20  kbps  digital  FM  signalling,  and  the  hop  rate  is 
100  hops/sec,  then  Q  could  be  as  high  as  200. 

When  the  Q  symbols  have  all  been  generated  at  rate  Re  and  stored  in  the  input 
buffer,  they  are  then  transferred  to  a  second  (output)  buffer.  The  transmitter  logic 
reads  this  buffer  L  times  at  the  rate  LRe,  and  this  stream  of  data  “chips”  is  used  to 
frequency-modulate  the  selected  carrier  frequency,  which  is  changed  (hopped)  to  a  new. 
pseudorandomly-selected  value  after  Q  chips  have  been  transmitted.  In  this  manner,  L 
copies  of  the  Q-symbol  sequence  have  been  transmitted  on  L  different,  successive  hops. 
Although  Figure  2-1  suggests  that  the  Q  symbols  are  repeated  in  the  same  order  on 
each  hop,  it  is  of  course  possible  with  a  more  sophisticated  system  to  permute  the 
symbols  or  otherwise  scramble  them  so  that  the  order  of  appearance  of  the  symbols  is 
different  on  each  hop. 
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We  note  that  certain  fundamental  relationships  exist  among  the  digital  rates  at 
various  points  in  the  transmission  logic,  and  among  the  energies  in  the  transmitted 
waveform  which  correspond  to  each  rate.  The  symbols  actually  transmitted  are  keyed 
at  the  rate  Rd  which,  as  we  have  already  noted,  is  a  basic  specification  of  the  communi¬ 
cations  channel  being  used  for  hopping.  The  original  symbol  rate  is  Rc  =  R d/L,  on 
account  of  the  repetitions.  Viewed  another  way,  the  energy  transmitted  per  chip  is  the 
fraction  1/L  of  the  coded  symbol  energy.  If  the  error  control  code  rate  is  r,  then  the 
original  bit  rate  is 

R(,  =  Rcr  =  rRd/L.  (2-1) 

For  example,  for  Rd  =  20kbps,  r=i,  and  L  =  5,  the  equivalent  coded  symbol  and  bit 
rates  are  Rc  =  4kbps  and  Rt  =  2kbps;  between  the  coding  and  the  repetitions,  the  bit 
energy  gets  split  into  ten  pieces  in  this  example. 

Having  briefly  gone  over  the  overall  transmission  scheme,  we  now  discuss  certain  of 
its  aspects  in  greater  detail. 

2. 1.1.1  Error  control  coding 

It  has  been  suggested  in  [1]  that  concatenation  of  a  convolutional  “outer"  forward- 
error-control  code  with  an  “inner”  repetition  code  will  yield  an  FH/CPFSK  performance 
in  worst-case  partial-band  noise  jamming  which  improves  with  the  number  of  hops  per 
bit,  L,  when  soft-decision  outer  decoding  and  hard-decision  inner  decoding  are  used. 
These  terms  will  be  discussed  further  below;  here  we  note  the  effect  that  use  of  L 
hops/bit  has  on  the  input  data  rate.  If  the  maximum  channel  data  rate  Rd  is  20  kbps 
using  CPFSK,  then  the  maximum  coded  symbol  rate  is  Rc  =  20kbps/£  and  the 
maximum  input  data  rate  is  R6  =  20kbps  •  r/L.  For  example,  with  r=  \  and  L  =  4.  the 
system  can  accommodate  a  bit  rate  no  higher  than  2.5kbps.  The  practical  effect  of  this 
limitation  is  that  the  coding  cannot  be  used  for  most  voice  transmissions,  since  the  bit 
rate  needed  to  support  voice  data  is  commonly  regarded  to  be  16kbps,  although  in  some 
situations  an  acceptable  digital  voice  quality  may  be  achieved  with  a  significantly  lower 
data  rate.  Therefore,  when  we  speak  of  coded  results,  including  diversity  transmissions, 
necessarily  we  are  considering  data  other  them  voice  data.  Encryption  (scrambling)  can 
be  used  for  voice  transmissions,  since  the  data  rate  is  not  affected;  such  a  procedure 
provides  a  measure  of  transmission  security  but  does  not  offer  any  improvement  in 
system  performance  against  jamming. 

For  this  study,  it  is  assumed  that  the  error  control  code  used  is  a  convolutional 
code  with  rate  r=\  and  constraint  length  7. 
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2. 1.1. 2  Inter- hop  switching  time 

Another  practical  consideration  which  relates  to  the  transmission  scheme  is  the  fact 
that  the  transmitted  data  rate  during  the  hop  dwell  times  is  not  Rd,  as  depicted  in  our 
simplified  diagrams.  Rather,  Rd  is  the  average  channel  data  rate  since  usually  between 
hops  (that  is,  while  the  synthesizer  frequency  is  changing)  the  transmitter’s  power  is 
smoothly  turned  off  and  then  restored  in  order  to  avoid  generating  spurious  harmonics 
or  splatter.  During  this  interhop  switching  period,  some  buffering  of  the  channel  data 
chips  is  required,  since  they  are  being  produced  during  this  period  when  there  is  no 
output.  It  also  follows  that  the  actual  transmission  of  the  chips  must  be  “bursted”  at  a 
slightly  higher  rate,  Rd  >  Rd. 

In  our  analysis,  for  the  most  part  it  is  permissible  to  overlook  this  implementation 
detail.  However,  in  Section  3  we  shall  be  studying  the  issues  connected  with  the  inter¬ 
hop  switching  time  in  detail,  in  order  to  assess  the  effects  it  has  on  the  system's 
performance. 

2. 1.1. 3  Transmitted  waveform  and  its  parameters 

The  coded  symbol  chip  stream  {dk}  going  into  the  modulator  in  Figure  2-1  is 
assumed  to  be  representable  by  the  bipolar  data  waveform 

d(t )  =  ]rd*p(t  — &T),  dk  =  ±1,  (2-2) 

k 

where  the  pulse  function  p(t)  is  assumed  to  be  rectangular: 

p(t)  =  u(t)-u(t-T)=  rect(t-|T),  (2-3) 

with  u(t)  being  the  unit  step  function.  The  interval  T  is  the  symbol  duration,  so  that 
the  data  rate  is  Rd  =  l/T.  Although  much  attention  has  been  given  to  reducing  the 
transmitted  bandwidth  of  binary  FM  signals  by  employing  non-rectangular  p(t)  pulse 
shapes  and  by  correlative  coding  or  trellis  coding  of  the  data,  here  we  will  assume  the 
p(t)  as  given  and  will  treat  the  dk  symbol  values  as  having  been  generated  inde¬ 
pendently.  The  resulting  modulated  carrier  frequency  is  given  by 

M  =  A  +  /d-d(t),  (2-4) 

and  the  commonly  accepted  measure  of  relative  frequency  deviation  is  the  modulation 
index  h,  defined  by  normalizing  the  deviation  by  the  data  rate: 

h  =  2fd/Rd  =2fdT.  (2-5) 

The  frequency  modulation  of  the  carrier  frequency  fh  by  the  data  waveform  d(t)  results 
in  the  constant  envelope  signal 
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s0(t)  =  const,  x  cos  [wAt  +  6  m  (t)],  (2-6a) 

where  the  data-modulated  phase  6  m  (t)  is 

M t)  =  2irfd  fl  d£  d(0-  (2-6b) 

-oo 

Some  filtering  of  the  transmitted  signal  may  be  done  in  order  to  contain  the  output 
spectrum;  however,  this  would  introduce  amplitude  and  phase  distortion.  Therefore, 
only  the  broadest  type  of  filtering  is  used  on  the  transmitter  output,  such  as  may  be 
needed  to  reject  modulation  products,  etc.  The  preferred  way  to  control  the  bandwidth 
of  the  emitted  spectrum  is  to  smooth  the  data  waveform,  as  in  “tamed  FM”  [2],  in  such 
a  way  that  the  resulting  intersymbol  interference  is  predictable.  In  our  analysis,  any 
effects  of  transmitter  filtering  are  neglected,  and  for  convenience  all  intersymbol 
interference  and  signal  distortion  effects  are  considered  as  resulting  from  the  receiver 
processing. 

2.1.2  Reception  Scheme  and  Received  Waveform 

Figure  2-2  provides  an  overview  of  the  FH/CPFSIv  reception  scheme.  Basically, 
the  processes  involved  are  the  reverse  of  those  for  the  transmission  scheme.  The 
receiver  must  accurately  synthesize  local  oscillator  frequencies  in  order  to  dehop  the 
signal  at  the  proper  times,  and  must  then  be  capable  of  sampling  the  demodulator 
output  at  the  end  of  each  chip  interval.  These  sample,  QL  of  them  for  L  complete 
hops,  are  buffered  so  that  the  receiver  logic  can  in  some  manner  combine  the  L  chips 
belonging  to  a  particular  code  symbol. 

As  the  figure  suggests,  binary  decisions  can  be  made  on  each  chip  as  it  is  received. 
The  resulting  logic  and  buffering  for  such  a  “hard  decision”  procedure  is  simpler  than 
that  required  for  a  “soft  decision”  scheme,  which  involves  A/D  conversion  of  the 
samples  and  storing  the  resulting  QL  multibit  words,  one  per  chip.  With  either  hard  or 
soft  chip  processing,  after  diversity  combining  there  is  the  option  to  perform  binary 
decoding  of  code  symbol  decisions  or  else  soft-decision  decoding  using  code  symbol 
metrics  produced  by  the  combining. 

2. 1.2.1  Received  signal  waveform 

Waveform  selection  in  the  receiver  is  accomplished  using  a  hopping  local  oscillator 
(synthesizer),  assumed  to  be  in  synchronism  with  the  hopping  pattern  of  the  received 
signal  waveform,  and  a  bandpass  filter  centered  at  the  IF  (intermediate  frequency)  of 
the  receiver.  Ideally,  the  signal  portion  of  the  IF  waveform  out  of  this  bandpass  filter  is 
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equivalent  to  that  for  a  system  whose  carrier  frequency  is  fixed,  not  hopping: 

s(t)  =  a(t)  Acos  [^'/ft  +  <£(t)j, 


where  fIF  =  uIF/2n  is  the  If  center  frequency  and  A  is  the  received  signed  amplitude, 
assumed  constant.  The  amplitude  modulation  a(t)  and  the  distortion  of  the  trans¬ 
mitted  phase  modulation  0m(t)  into  the  phase  function  ^(t)  axe  effects  induced  by  the 
IF  filter.  These  quantities  are  given  by 


and 


a(t)  =  u2(t)  +  v2(t) 


<^>(t)  =  tan 


u(t)‘ 

v(t)J’ 


(2-8a) 
(2— 8b) 


where 

u(t)  =  ho(t)*sin0m(t)  (2-9a) 

and 

v(t)  =  h0(  t)  ^  cos  0  m  ( t )  (2-9b) 

are  the  quadrature  and  in-phase  components,  respectively,  of  the  IF  filter  output.  In 
these  expressions,  (*)  denotes  convolution,  and  h0(t)  is  the  equivalent  lowpass  impulse 
response  of  the  IF  filter,  that  is.  that  filter  is  modelled  by 

h/f(t)  =  2h0(t)cosw/ft  (2-  10a) 

or,  in  the  frequency  domain. 

=  H0(/-//f)  +  H  o(/+//f).  (2-10b) 

For  example,  if  the  lowpass  filter  has  the  Gaussian  shape 

H0(/)  =  e-nf7/SB\  (2-11) 

then  the  3dB  bandwidth  B3  is  0.9294  times  the  lowpass  noise  bandwidth  B0  =  j\V Jf. 
For  an  n-pole  Butterworth  filter  with  transfer  function 

'"•wr  =  i  t  (kp  (2"12a) 

the  noise  bandwidth  is  related  to  the  3  dB  bandwidth  by 

B»  - <2-12b) 

=  1.5708B3,  n=  1 

=  1 . 1 1 07 B 3 ,  n=  2 
=  1.0472B3,  n  =  3 
=  1.0262B3,  n  =  4. 
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Now  since  theoretically  the  FM  signal  waveform  has  infinite  bandwidth,  the  IF 
filtering  rejects  not  only  unwanted  signals  but  also  portions  of  the  desired  signal. 
Therefore  the  filtering  introduces  distortion  in  the  form  of  amplitude  modulation  and 
phase  distortion,  as  shown  above.  For  the  common  modelling  assumptions  we  are 
making,  in  effect  the  distortion  caused  by  the  finite  IF  bandwidth  represents  all  the 
distortion  suffered  by  the  waveform  — at  least  all  distortion  due  to  filtering— just  as  the 
noise  at  IF  represents  all  noise  present. 


2. 1.2. 2  Noise  and  jamming 

The  toted  IF  waveform  x(t)  noted  in  Figure  2-2  is  assumed  to  consist  of  the  signal 
term  given  by  (2-7)  plus  stationary,  additive  Gaussian  noise:  x(t)  =  s(t)  +  n(t).  When 
partial-band  noise  jamming  is  present,  covering  the  fraction  7W  of  the  total  hopping 
bandwidth  W,  with  0<7<1,  the  noise  power  in  the  IF  bandwidth  is  modelled  as 
random  on  account  of  the  signal's  hopping  around  in  the  total  bandwidth,  and  has  a 
two-valued  distribution: 

E{n2(t) }  =  <tjv  with  probability  I-7  (2-  12a ) 

=  <r-N  +  cr-j  =  with  probability  7.  ( 2 -  12b ) 

Assuming  both  types  of  input  noise- the  combined  background  and  receiver  input 
noise,  represented  by  cry.  and  the  jamming  noise,  represented  by  a;)  -  have  flat  spectral 
power  density  functions  going  into  the  IF  filter,  the  respective  noise  powers  are 

<■  =  *o  (2-1 3a ) 

and  c7j  =  K  j\\  JF/~,.  (2-13b) 

In  (2-13b),  the  spectral  level  Kj  is  the  average  jammer  noise  spectral  power  density 
referred  to  the  total  bandwidth,  that  is, 

Kj  4  J/W,  (2-13c) 

where  J  is  the  total  jammer  power  at  the  receiver  location. 


The  signal-to-noise  power  ratio  (SNR)  at  the  output  of  the  IF  filter  without  and 
with  jamming  are  time-varying  due  to  the  amplitude  distortion  o(t)  as  a  function  of  the 
data  modulation.  For  a(t)  =  1,  we  define  the  SNR’s 


0  A  _Ai  _  _! _ Erf  _  l.r.E» 

P*  ~  2a?.  ~  \\’jfT  *o~  DL^o 

and 


(2-14) 


.  A  _Al  _  •>  E 4  _  7  r  Efc 
Pj  ~  2a?  “  \V/fT'k;  "  D 'i'Jf/ 


(2-15) 
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In  (2-14)  and  (2-15)  for  convenience  we  use  the  symbol  D  to  denote  the  IF  filter  s 
bandwidth-time  product: 

D  ^  W/fT.  (2-16) 

How  the  phase  of  the  total  waveform  x(t)  is  affected  by  the  noise  can  be  seen  by 
combining  signal  and  noise  terms  at  the  output  of  the  IF  filter  to  obtain 

x(t)  =  a(t)  A  cos  [uijf  t  +  </>(t)]  +  n(t) 

=  [v(t)  +  m(t)]cosw;/-t  -  [u(t)  +  n,(t)]sinw;/-t 

=  R(t)cos|u.;/-t  ■+■  4>(t)],  (2-1  /a) 

in  which  R(t)  and  <I>(t)  are  envelope  and  phase  of  the  total  waveform,  and  we  have 
utilized  the  quadrature  or  Rician  decomposition  of  bandpass  noise: 

n(t)  =  nf(t)coSu./ft  -  ns(t)sin uijpt.  (2-17b) 

We  may  define  a  phase  noise  term  by  writing 

tan <t>(t )  =  Ujt>  ^  =  tan (o(t )  +  »>(t)]:  (2-1S) 

v(t)  +  nc(t) 

the  phase  noise  term  thus  defined  is 


ry(t )  =  tan 


-l 


tan  <!>  -  tan  o 
1  +  tan<t>tano 


=  tan' 


n.v  -  nru 


u*  +  v2  +  n<-v  +  n,.u 


=  tan 


-l 


MO 


a(t )  «j2p  +  Ml) 


(2-19) 


in  which  p  is  the  SNR.  and  and  v,  at  the  same  time  instant  are  the  independent, 
zero-mean,  unit-variance  Gaussian  random  variables  given  by 

av,{i)  =  n,(t)cosO(t)  -  ne(t)sin0(t)  (2-20a) 

<r«/c(t)  =  n.(t)sintf(t)  +  ne(t)cos<£(t).  (2-20b) 

From  observation  of  (2-19)  we  may  state  that  for  low  SNR,  r?  tends  to  be  uniformly 
distributed,  since  for  this  case  r?  is  approximately  the  random  phase  of  the  noise.  For 
high  SNR.  the  distribution  of  r?  tends  toward  a  truncated  or  ‘'aliased''  modulo-27r 
Gaussian  distribution,  since  tan~‘z  s=  x  for  r<  1.  We  also  observe  from  (2-19)  that  this 
additive  phase  noise  is  signal-dependent. 
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2.1.3  Demodulator  Options 

Two  types  of  FM  demodulator  are  considered  in  this  study:  limiter-discriminator 
and  differential  detector,  as  illustrated  in  Figure  2-3. 

Essentially  the  limiter-discriminator  extracts  the  instantaneous  frequency  deviation 
from  the  carrier  (or  IF  center  frequency),  4>(t),  and  is  followed  by  an  integrate-and- 
dump  lowpass  filter.  The  receiver  output  produced  is  given  by  the  differential  phase 
sample  value  at  tk  =  kT  +  t0 

A4>(t*)  =  *(tt)  -  *(tk-T).  (2-21) 

The  differential  detector's  output  is,  when  the  delay  r  equals  T, 
fR(ti)R(t*  —  T)  cos  [w//-tj.  +  4>(tt)]  sin  (t*  —  T)  +  4>(tt  —  T)]^j 

lowpass 

=  |  R ( 1 1 )  R(t k  —  T)  sin  [ A4>( t *. )  +  T] 

=  |R(t*)R(tJfc  —  T) sin  [A<J>(tk )].  uJ/pT  =  2nn.  (2-22 j 

Ekanavake  [4]  indicates  that  differential  detector  performance  in  terms  of  the  proba¬ 
bility  of  correct  chip  decisions  can  be  made  comparable  with  that  of  the  limiter- 
discriminator  by  choosing  r  <  T.  For  r  =  T.  it  has  been  shown  [5]  that  the  differential 
detector  achieves  about  the  same  performance  as  the  limiter-discriminator  receiver 
when  the  modulation  index  h  <  0.5.  but  is  degraded  significantly  for  larger  values  of  h. 


2. 1.3.1  Differential  phase  in  the  absence  of  noise 

The  choice  of  the  sampling  time  tk  is  based  on  sampling  when  the  differential  phase 
in  the  absence  of  noise, 

Ao(t)  =  <£>(t )  —  <£(t  — T),  (2-23) 

has  a  positive  or  negative  peak  value.  In  general,  the  best  sampling  time  for  the  dis¬ 
torted  phase  is  different  than  that  for  the  original,  undistorted  phase  modulation. 

The  undistorted  signal  phase  waveform,  0m(t),  at  the  particular  instant  t,  =  fciT  +  (. 
can  be  written  as 

0m(t])  =  jt  h'^2<ii  +  <4jf ;  (2-24a) 


its  differential  value  is 

Aflm(t,)  =  Mt,)  -  Mtl-T) 

=  ^*,-1  +  ^(<4,- ( 2 -24  b ) 

It  is  evident  that  when  successive  data  bit  values  are  different  (dt  ^  a  maximum 
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(a)  Limiter-discriminator  followed  by  an  integrate-and-dump  filter 


(b)  Differential  detector  for  narrowband  FM 
FIGURE  2-3  DEMODULATOR  OPTIONS 


19 


Analysis  Models 


deflection  of  ±vh  would  be  achieved  the  undistorted  phase  using  a  sampling  time  at  the 
end  of  the  data  symbol  interval  (i.  e.,  c  =  T). 

As  shown  in  Appendix  A,  the  effect  of  the  IF  filter  is  to  produce  the  distorted  and 
delayed  differential  phase  waveform  A<f>(t).  The  delay  causes  the  maximum  deflection 
of  A <f>  to  occur  at  the  sampling  time  offset  e  =  T  +  trf  >  T,  where  td  is  the  delay,  and 
also  causes  the  magnitude  of  the  deflection  to  be  smaller  than  irh.  The  distortion 
induces  a  dependence  of  the  sampled  A<j>  value  on  the  data  sequence,  due  to  intersymbol 
interference  (ISI).  For  the  cases  of  practiced  interest,  only  ISI  from  immediately 
adjacent  chips  is  significant  [6]. 


2. 1.3. 2  Probability  density  function  for  the  limiter-discriminator  receiver  output 

It  is  assumed  that  the  discriminator  output  y(t)  illustrated  in  Figure  2-3  is  the 
derivative  of  the  signal's  instantaneous  phase  with  respect  to  the  center  frequency,  }IF : 

y(t)  =  4>(t) 

d  t  n-if A  a(t)sind(t)  4-m(t) 
dt  A  a(t)  cos<p(t)  -f  n«(t) 

=  ^L{<i(t)A  a(t)[Aa(t)  +  m(t)cos<*>(t)  +  n,(t)sin4>(t)J 
+  A  a(t)  [nc(t)sin<2>(t)  -  n,(t)cos«i>(t)] 

+  nc(t)hi(t)  -  hc(t)ni(t)|.  (2-25) 

Because  of  the  presence  of  the  Gaussian  noise  terms,  from  these  expressions  we 
recognize  the  fact  that  the  discriminator  output  is  a  random  variable  wrhich  can  take 
any  value  on  (  — oc  .  +oo). 

The  receiver  further  integrates  the  discriminator  output  over  a  symbol  interval  to 
obtain  the  receiver  output 

Z(t)  =  f'  dO-(0  (2-26a) 

it-T 

=  $(t)  -  <*>(t-T)  =  A4»(t),  (2-26b) 

where  A<t>  is  the  total  differential  phase.  Note  that  this  differential  phase  can  take 
values  on  (  — oc, +oo),  as  observed  in  [7];  however,  since  the  phases  are  expressible  only 
indirectly  via  the  arctangent  function,  for  example, 


A4»  =  tan 


tan<J>(t)  -  tan$(t-T)  ) 

1  +  tan$(t)tan  <t>(t-T)J' 


(2-27a) 
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(2-27b) 


it  is  customary  to  resolve  the  ambiguity  of  this  function  by  writing 

A4>  =  y  +  2*(N+-  AL), 

where  the  first  term 

y  =  (A#)  modulo  2t  (2-27c) 

is  the  principal  value  of  the  arctangent,  and  N+  and  N _  are  positive  and  negative  “click 
numbers,”  respectively. 

Click  statistics.  For  moderate  and  high  values  of  SNR,  it  has  been  shown  [8,  9] 
that  the  distribution  of  the  difference  of  click  numbers,  N+—  AL,  is  well-represented  by 
that  of  the  random  variable  Nc  whose  absolute  value  is  a  Poisson  random  variable  and 
whose  mean  tends  to  be  of  the  opposite  sign  of  A<6,  the  differential  phase  in  the  absence 
of  noise.  That  is, 

Pr{  \N<  \  =  n]  =  ^-j-exp{-[a  |  },  n=0,  1,  2,  ...  (2-28a) 

where 


Q  =  Nc  =  — 


,-rd<  -Zf  ' 


<^(0  -pa2(() 


(2-28b) 


This  click  average  is  consistent  with  the  fact  that  the  mean  value  of  the  discriminator 
output  is  [10] 

E{  v(t) }  =  y(t)[l-e'pa2(t)];  (2-28c) 

exchanging  the  order  of  integration  and  expectation  yields  the  following  mean  value  for 
A<I>(t): 


E{  z  }  =  P  d^U)-f'  d$y(0e^aU) 

Jl-T  Jt-T 


=  Ay(t)  +  2?r  Nc.  (2-2Sd) 

Distribution  of  principal  value.  The  probability  density  function  for  y  has  been 
found  to  be  [11-13] 

pv,(i)  =  p^(x;  Ay,  r,  A,  U,  V,  W) 


_  1  — r2  — A2 


do 


e _ cos  q 


[l-(rcosx+ Asinxjcosa]2 

U-W(rcos Ay  +  A  sinAy) 


l-g{x,a)+2- 


1-r2- A2 


(2-29a) 


in  which 


g(x,  a)  = 


U -\sino  — W  cos  (x- Ay)  cos  q 
1  —  ( r  cos  x-f  A  sin  x)  cos  a 


(2-29b) 
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U  =  ip[a2(t)  +  a2(t-T)]  ,  (2-29c) 

V  =  §^[a2(t)  —  <*2(t— T)]  ,  (2-29d) 

W  =  J U2- V2  =  P  Ja(t)a(t-T)  (2-29e) 

r  =  <r"2E{nc(t)ne(t-T)}  =  r-2E{n,(t)n.(t-T)},  (2-29f) 

and  A  =  <7-2E{nt(t)ns(t-T)}  = -o-_2E{n.(t)nc(t-T)}.  (2-29g) 


The  domain  over  which  (2-29a)  is  valid  is  usually  taken  to  be  jx—  A<i>|  <  ir.  The 
quantities  U  and  W  can  be  seen  as  the  arithmetic  and  geometric  averages,  respectively, 
of  the  SNR’s  at  times  t  and  t  — T.  The  parameters  r  and  A  are  the  in-phase  and  quadra¬ 
ture  correlation  coefficients  of  the  noise  at  the  output  of  the  IF  filter,  and  are  a  function 
of  that  filter,  assuming  a  flat  noise  spectrum  going  into  the  filter.  For  example,  if  the 
filter  has  a  symmetric  passband  about  the  IF  frequency,  a  result  is  that  A  =  0. 

The  total  differential  phase  for  the  output  of  the  limiter-discriminator  receiver  has 
the  probability  density  function 

Pa*{x)  =  exp { —  |  a  | }  £  *-~J-p*(x±2n7r),  (2-30) 

where  a  is  the  click  number  average  given  in  (2-28b)  and  the  sign  taken  is  the  same  as 
that  of  Ae>.  Figure  2-4  illustrates  the  multi-model  character  of  pa#(x). 

When  the  receiver  is  subject  to  partial-band  jamming  with  probability  ■>,  by  exten¬ 
sion  (2-30)  becomes,  for  AcJ>0. 

->)  =  (1  — £“TP*(*+2nx;  pN) 

n  =  0  71  ’ 

+  7e'ai  £^p,(x+2mr;  pT)  ,  (2-31) 

n  =  0  n' 


using  aQ  and  to  denote  the  absolute  values  of  click  number  averages  under  no 
jamming  and  jamming  conditions,  respectively.  The  SNR  in  the  absence  of  jamming. 
pN,  was  defined  in  (2-14);  the  jammed  SNR,  pr,  is  given  by 


A  Ai  _  A2 
Pt~  24  2(4+t5) 

E».  E> 

-  ps  pJ  1 .  r  . 

Py  +  Pj  D  L  E»' 

*0  + 


(2-32a) 


(2-32b) 
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2. 1.3.3  Probability  distribution  for  the  output  of  a  differential  detector 

The  output  of  the  differential  detector  shown  in  Figure  2-3(b)  is 

z(t)  =  |R(t)R(t- T)sin[4>(t)-<l>(t- T)],  uIFT  =  2nn.  (2-33) 

It  was  shown  in  [8]  that  the  output  of  a  differential  detector  at  the  sampling  time  is 
equivalent  to  the  difference  between  two  independent,  scaled  noncentral  chi-squared 
random  variables,  each  with  two  degrees  of  freedom.  We  denote  this  equivalence  by 


z  ~  C!x2(2;  dj)  —  c2 x2(2;  d2), 
where  the  scale  factors  are 

c12  =  \o7{  -J 1  — r2  T  A) 
and  the  noncentrality  parameters  are 

A  _  2{U— rWcos  Ao  ±  -Jl—  r2  Wsin  A<^} 

12 "  ' 

Therefore,  the  characteristic  function  for  the  random  variable  z  is 
=  C^2(civ;  2.  dj  C^(-c2w,  2,  d2) 


(2-34a) 

(2-34b) 

(2-34c) 


1  +4j/*C]C2+2ji/(c2— Ci) 


fM<udi 
1  1  + 


Ij— c2d2)  — 2t,2C]C2(di-(-d2) 
+4j/2c1c2+2j^(c2-Ci) 


_ . _ 1 _ _  eJ  >C-»2D  \ 

1  +ji^A  +  ^*B  '  |l  +ji/A  +  t/2Bj 


(2-35a) 


where  for  convenience  we  have  defined 


A  4  2(c2-Ci)  =  cr2A 

(2-35b) 

B  4  4c,c2  =  i(l-r2-A2) 

(2-35c) 

C  =  c1d1-c2d2  =  «r2WsinA  4> 

(2-35d) 

and 

D  =  2cic2(d!+d2)  =  5<T’4(U-rWcosA<^-fAWsinA<^). 

(2-35e) 

The  mean  value  of  z  is  easily  obtained  from  (2-34a)  and  is 


E{z}  —  Cj(2-f-di)  —  c2(2  +  d2) 

=  2(cj  c2 )  -f  Cjdj  —  c2d2 

=  <72(A  +  W  sin  At?)  (2-36a) 

=  l  A2,ja(t)  a(t-T)sin  A<t>  forcr2=0.  (2-36b) 
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Since  we  will  be  interested  in  the  probability  of  error  for  various  weighted  combinations 
of  differential  detector  output  samples,  the  characteristic  function  is  a  more  convenient 
description  of  the  distribution  than  the  probability  density  function  (PDF),  for  reasons 
which  will  be  shown  in  Section  4.  For  further  developments  of  the  PDF  for  the  class  of 
cross- correlators  to  which  the  differential  detector  belongs,  the  reader  is  referred  to  [14]. 

2.1.4  Generic  Chip  Combining  Schemes 

Our  analysis  efforts  are  focused  on  ways  to  combine  the  L  chips  belonging  to  a  par¬ 
ticular  data  symbol  so  as  to  give  a  diversity  gain,  if  possible,  against  worst-case  partial- 
band  noise  jamming.  A  generic  formulation  for  the  chip  combining  approach  is  to  base 
the  symbol  decision  on  the  polarity  of  the  decision  statistic  or  metric 

z*  =  £<*,  =  £/(**.)’  (2-37) 

«=i  ?=i 

where  /( zkq)  is  some  function  of  the  gth  chip  belonging  to  the  &th  data  symbol.  For  ex¬ 
ample.  if  hard  decisions  are  made  on  each  chip,  then 

C  ir  q  —  —  1, 

=  0,  z*,<0.  (2-3S) 

This  nonlinear  combining  function  is  known  to  produce  a  diversity  gain.  In  our  pre¬ 
vious  study  [8.15],  we  investigated  the  combining  based  on  /( zkq)=zkv  or  linear 
combining,  and  found  that  this  method  does  not  yield  a  diversity  gain. 

In  general,  in  order  for  a  diversity  gain  to  be  realized,  it  is  necessary  for  the 
receiver  to  give  a  lower  weight  in  the  sum  (2-37)  to  chips  which  have  undergone  jam¬ 
ming  as  the  consequence  of  the  hop  on  which  they  were  transmitted  having  hopped  into 
a  portion  of  the  band  which  is  jammed.  Thus  it  is  required  that  the  function  C,kq  of  the 
chip  decision  variable  zkq  somehow  discriminate  against  jammed  chips.  This  concept 
can  be  expressed  by 

/(zif;  no  jamming)  >  /( zkq\  jamming).  (2-39a) 

In  principle,  this  requirement  can  be  satisfied  by  weighting  ziq  by  the  ratio  of  the  un¬ 
jammed  noise  power,  <r^,  to  the  noise  power  on  its  hop,  giving 

Ct?  =  <^zt,/<r?  =  wfz*f  (2-39b) 

=  zkq,  hop  q  not  jammed  (2-39c) 

=  zt,(T^./(T2r,  hop  q  jammed.  (2-39d) 

Techniques  for  measuring  the  power  on  the  hop  are  discussed  in  [24]. 
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2.2  FORMULATIONS  FOR  THE  PERFORMANCE  MEASURES 

In  this  report  the  performance  of  the  system  is  measured  in  terms  of  the  proba¬ 
bility  of  a  bit  error  for  FH/CPFSK  as  a  function  of  jammer  power  when  the  system  is 
subjected  to  worst-case  partial-band  noise  jamming,  with  and  without  error-control 
coding,  for  different  diversity  combining  schemes  and  numbers  of  hops  per  transmitted 
symbol,  L. 


2.2.1  Uncoded  Bit  Error  Probability 

Conceptually,  the  probability  distribution  for  the  output  of  the  combiner  at  a  given 
time,  z,  is  conditioned  on  the  number  of  hops  which  have  been  jammed.  Let  p(z  j  Z ) 
denote  the  PDF  of  z  conditioned  on  there  having  been  /  out  of  L  hops  jammed.  Then, 
averaged  over  the  possible  jamming  event  (i.e.,  the  possible  values  of  /),  the  probability 
of  an  error  in  deciding  on  the  polarity  of  a  transmitted  binary  symbol  is 

P<  =  Pr{<4=  0}  Pr{z*>  0 1  dt  =  0}  +  Pr{<4  =  1}  Pr{zt<  0  |  <4=1} 

=  \  [Pr{z*> 0 1  <4  =  0}  +  Pr{z*<0|  <4=  1}] 

=  Pr{zt<  0  |  <4=  1} 

f0 

=  )  Pr{l  hops  jammed}  dz  p(z  | /).  (2-40) 


This  formulation  assumes  that  the  two  possible  binary  symbol  values  are  equally  likely, 
and  that  the  conditional  probability  of  error  is  the  same  for  both  <4=0  and  <4=1. 
Assuming  that  a  fraction  7  of  the  available  hopping  channels  are  jammed  and  that  the 
pseudorandomly-selected  hop  frequencies  are  equally  likely,  the  probability  that  l  out  of 
L  hops  are  jammed  is 

Pr{I  hops  jammed}  =  (^)V  (1— (2-41) 

Note  that  this  expression  represents  the  total  number  of  jamming  situations  resulting  in 
l  hops  jammed,  irrespective  of  which  of  the  hops  are  jammed,  or  in  what  order.  Thus 
(2-39)  is  valid  if  the  combiner  output’s  probability  distribution  does  not  depend  on  the 
order  in  which  the  hops  were  jammed  or  not  jammed. 


Since  the  chips  being  combined  axe  independent  by  virtue  of  being  corrupted  by 

noise  in  different  frequency  channels,  their  joint  PDF  is 

l 


p(ztl,  zt2,  ...,  zkL )  =  If  p?(zt?).  (2-42) 

?=i 

For  convenience,  we  may  arbitrarily  designate  the  first  l  samples  as  occurring  on 
jammed  hops. 
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For  an  unjammed  hop  the  noise  power  on  the  hop  is  <r2=o-2r,  and  for  jamming  it  is 
cr7  =  a7T=  where  we  adopt  the  notation 


,S'4:  =  i4  =  i+  . 

°N  7E  JHj 


Using  this  notation  it  is  convenient  also  to  write 

P«(z*4)  =  P (z*,;  *«) 

=  P(zi9! 

=  p(zt?;  k*n), 

This  notation  also  allows  us  to  rewrite  (2-42)  as 


hop  q  unjammed; 
hop  q  jammed. 


P(Z*1-  Z*2>  Zti)  =  P(Zi9)  fj)  II  P?(Z*5’  KCT/</)- 

unjammed  jammed 

L-l  hops  /  hops 


(2-43) 


(2-44a) 
(2 -44b) 

(2-45) 


The  characteristic  function  for  the  combiner  output  is,  using  a  similar  approach  to  the 
notation, 

1  r*  » 


C.-(i/;  l)  =  E  exp{jj/£  /(zt9)} 

«=i 


=  |C<(";  «»£.))'  |  <?<(*-;  (2-46) 

in  which  C((i/;  a2)  is  the  conditional  characteristic  function  for  C*9  =  /( zt?)  when  the  hop 
noise  power  is  a'. 

If  the  mathematical  formulation  for  C((»/;  a-2)  is  simple,  the  preferred  method  for 
obtaining  the  probability  of  error  is  first  to  obtain 

p(z|f)  =  H  di/  <T‘"C.{k  /)  (2-47) 

J  -00 

for  substitution  in  (2-45). 


2.2.2  Coded  Bit  Error  Probability 

The  decoding  schemes  we  will  be  considering  in  this  report  include  both  hard- 
decision  decoding  and  soft-decision  decoding.  By  “hard”  we  mean  that  the  inputs  to 
the  decoder  are  binary  symbol  decisions.  “Soft”  refers  to  a  non-binary  input  to  the 
decoder.  For  the  selected  constraint  length  seven,  rate  one-half  convolutional  code,  the 
error  probability  for  the  decoded  bits  is  bounded  by  [25] 
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Pj  =  P b(L;  7,  ECCM  combining  scheme,  decoding  input) 

<  36Qio  +  211 Q12+  1404 Q14+  11633Q16+  •••  (2-48) 

in  which  the  quantities  Q2m  are  the  probabilities  of  decoder  error  in  comparing  the 
correct  path  segment  to  a  path  segment  that  differs  in  2m  symbols.  If  the  decoder  is 
presented  with  hard  decisions  {dk}  on  the  transmitted  code  symbols  {dk},— with 
average  error  probability  P«(£)  =  Pe(I;  7,  ^Ek)-the  Q2m  are  given  exactly  by 

=  g(2f)[P.(L)]!”-*[l-P.(i)]‘+K2m  ){[1-P.(£)1P.(£)}”  (2-49.) 

-  (JmJkl  to  (2-49b) 

where  the  form  given  by  (2-49b)  is  derived  in  Section  4.2.1.  For  example, 

Q2  =  P»(I),  (2-50a) 

Q4  =  3[Pe(I)]2 -2[Pe(E)]3,  (2-50b) 

Q6  =  10  [P<(L)]3  — 15  [Pe(L)]4  +  6  [Pt(L)]5,  etc.  (2-50c) 

Using  A  =  Pe(L),  the  polynomials  used  in  (2-48)  are 

Q,0  =  126  A5  —  420  A6+  540  A7  -  315  A8  +  70  A9  (2-51a) 

Q12  =  462A6-  1980  A7 +  3465  A8 -3080  A9  +  1386A10-252AU  (2-51b) 

Qh  =  1716A7-  9009  A8  +  20020A9-  24024  A10 -i-  16380  A11 

-  6006  A12+  924  A13  (2-51c) 

Q16  =  6435  A8 -40040  A9  +  IOSIO8A10-  163800AU+  150150A12 

-  83160  A13  +  25740  A14  -  3432  A15.  (2-51d) 

Instead  of  presenting  the  hard  symbol  decisions  {dt}  to  the  decoder,  the  inputs 

may  be  a  soft-decision  metric  such  as  the  number  of  the  L  chip  decisions  {dkt\  g=l,  2, 
... L }  which  are  positive  (or  conceptually,  their  unthresholded  sum).  In  that  situation,  it 
is  shown  in  [25]  that  the  value  of  the  Q2m  used  in  (2-48)  are  bounded  by 

0.2m  <  ${4P«(1;  7,  £E*)[1-P«(1;  7,  ^Et)]}m£,  (2-52) 

in  which  it  is  emphasized  that  each  chip  represents  the  fraction  r/Z,  =  1/2X  of  the 
energy  of  an  original  encoded  bit. 
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2.3  SUMMARY  OF  CALCULATION  METHODOLOGY 

In  general,  we  have  found  it  more  convenient  to  use  the  characteristic  function 
method  in  studying  sums  and  weighted  sums  of  L  receiver  samples.  In  the  following  we 
summarize  how  the  probability  of  bit  error  may  be  calculated  from  the  characteristic 
function. 

2.3.1  Characteristic  Function  Method  for  Obtaining  the  Bit  Error  Probability 

The  cumulative  probability  distribution  P(Z)  =  Pr{z<Z}  for  a  random  variable  z 
may  be  written  in  terms  of  its  characteristic  function  C,( v)  as  follows  [16]: 

Pr{z  <  Z}  =  1-1  P#Im{CJ(t/)e-jt'Z}.  (2-53) 

0 

Therefore,  the  bit  error  probability  for  a  binary  system  with  decision  statistic  z  is 

Pt  =  Pr { z  < 0 1  <4=1}  =  (2-54) 

J  0 

When  z  is  a  weighted  sum,  that  is. 

L 

z  =  ]Tw,z?.  (2-55) 

?=i 

the  probability  becomes 

p<  =  tt  (II I  c'*»(J/W«)l)sinf]Car8[C,*»(*'w»)])-  (2-56) 

j0  ?=i  '  \=i  ' 

2.3.2  Characteristic  Function  with  FM  Clicks  Included 

In  Section  2.1,  a  general  expression  for  the  probability  density  function  of  the  differ¬ 
ential  phase  sample  A <t>  at  the  output  of  the  limiter-discriminator  receiver  was  shown  to 
be  a  weighted  sum  of  shifted,  modulo-2^  PDF’s,  p^(-).  The  corresponding  characteristic 
function  is 

CWtd  t)  =  [  dzp^z;  7)  ei"r 

J  —  oc 

=  (l-7)e_Q°  [  dz  e^rp^(z+2n7r;  pN) 

n  =  0  U’  J 

+  7e  f  dz  J1'1  pv,(z+2riff;  pT) 

n  =  0  U-  J  -oc 
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=  (l-7)e-U°£ 


-aovS'aJ 


n\ 


pn) 


, _  -Qi  v'  Q?  -iv2xn  n  ,  > 

+  7 e  2^  -j  e  w(^,  Pr)- 

n  =  0  Ti‘ 


(2-57) 


The  individual  characteristic  functions  C^( u;  p)  may  be  factored  out  of  the  sums,  and 
the  sums  completed,  to  yield  the  closed  form 

CA*( v,  7)  =  (l-7)exp{-o0+Q0  e~JI/2,r}  C^(v;  pN) 


+  7  exp{-aj  +Qj  e~j*/2,r}  Q(  1/;  pT) 

z -1  ,  -a0(l-cos2)ri/)-ja0sin2ffi/  ~  . 

C 1  —  7 )  e  W(i/;pA) 

+  7e-«.0-c~2^)-j<.,sin2„c,  ; 


2.3.3  Characteristic  Function  with  Adaptive  Weights 


(2-58a) 


(2-58b) 


If  a  weight  w0  is  given  to  a  A4>  sample  from  an  unjammed  hop,  and  to  one  from 
a  jammed  hop.  the  characteristic  function  for  A$  becomes 


r  1  .  \  n  x  -a0(l-cos27ru0J/)-jQ0sin27ru0i/ .  . 

7)  =  (l-7)e  cu.{w0u;  PA.) 

-Qi(l-cos2xu',i/)-ja,sin27ru)1i/  ,  , 

+  7e  1  1  1  C^u;  pT). 


(2-59) 
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3.0  FH/CPFSK  IMPLEMENTATION  STUDIES 

In  this  section  we  consider  certain  aspects  of  the  implementation  of  frequency- 
hopped  CPFSK  systems.  First,  in  Section  3.1,  we  discuss  various  issues  affecting  the 
implementation  of  current  frequency-hopped  systems.  Then,  in  Section  3.2,  we  use 
models  of  synthesizer  switching  transient  waveforms  to  analyze  their  effect  on  the 
emitted  signal  spectrum  and  the  de-hopped  receiver  IF  signal  during  the  switching 
times  between  hops. 

3.1  DISCUSSION  OF  IMPLEMENTATION  ISSUES 

The  implementation  issues  that  are  the  focus  of  the  following  discussions  include 
factors  affecting  the  selection  of  a  hopping  rate,  emissions  control  measures  and  hop 
transient  effects,  and  synchronization  requirements. 

3.1.1  Factors  Affecting  the  Selection  of  a  Hopping  Rate 

Certain  fundamental  limitations  on  frequency  hopping  rate  follow  from  system  re¬ 
quirements.  Suppose  that  the  maximum  signalling  rate  which  can  be  supported  by  a 
hopping  slot  or  channel,  in  order  to  avoid  adjacent  channel  interference  (for  example),  is 
Re  Also,  let  Rd  denote  the  rate  of  the  data  symbols  to  be  sent  over  the  hopping 
channel,  and  let 

T,,  =  R^1  =  hop  duration  or  dwell  time  (3- la) 

Q  =  RdT^  =  number  of  data  symbols  per  hop  (3- lb) 

To  =  portion  of  T^  needed  for  overhead  operations.  (3-lc) 

By  “overhead  operations”  we  refer  to  the  various  practical  operations  which  must  be 
implemented  to  support  the  transmission  of  a  frequency  hop;  these  include 

•  retuning  the  transmitter  and  receiver  synthesizers, 

•  controlling  the  fall  and  rise  of  the  transmitted  power  between  hops  to 

minimize  spectral  splatter, 

•  initializing  the  state  of  the  receiver  filters,  and 

•  acquiring  synchronization  on  the  hop,  if  necessary. 

During  the  portion  of  the  hop  interval  in  which  the  Q  symbols  of  user  data  are 
being  transmitted,  the  transmitted  symbol  rate  R,  must  be  higher  than  Rd  because  of 
the  overhead  operations.  Thus,  the  actual  symbol  rate  on  the  hop  is 
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R,  =  -  G.  =  M.  <  Re.  (3-2a) 

l^—lo  1  h—  io 

Solving  for  Rh  gives  the  inequality 

R»  <  (3-2b) 

where  RCT0  is  the  number  of  symbols  which  could  have  been  transmitted  during  the 
overhead  time  at  the  maximum  hop  channel  rate. 

For  example,  if  Rc  =  20kbps,  Rd=  16kbps,  and  To=0.2ms,  the  maximum  hopping 
rate  would  be  1000  hops/sec;  if  To=  2.0ms,  the  maximum  hopping  rate  would  be  100 
hops/sec.  Typically,  the  “off  the  air”  times  of  tactical  VHF  hopping  radios  is  on  the 
order  to  10-20%  of  TA  [26,  29]. 

If  the  data  rate  Rd  is  low  compared  to  the  channel  rate  Rc,  the  upper  limit  to  RA 
approaches  To1,  and  it  is  possible  to  have  the  situation  where  the  hopping  rate  is 
greater  than  the  data  rate,  commonly  referred  to  as  fast  hopping.1  A  fundamental  lim¬ 
itation  which  is  not  expressed  by  the  inequality  in  (3-2b)  is  that  the  hopping  rate  must 
be  less  than  Rc.  since  for  fast  hopping  the  hopping  rate  is  dominant  in  determining  the 
hop  bandwidth. 

It  has  been  said  [28]  that  the  SINCGARS  hop  rate  of  about  100  hops  per  second 
“was  determined  by  factors  such  as  signal  density,  assurance  of  communications  in  a 
mutual  interference  environment,  ability  to  enter  the  net  when  out  of  synchronization, 
technology  and  cost  limitations  of  developing  a  suitable  synthesizer,  and  reduced  mean 
time  between  failures. v 

On  the  other  hand,  and  Italian  manufacturer  claims  [27]  that  the  use  of  a  direct 
sequence  spread  spectrum  modulation  on  the  hops,  rather  than  the  commonly  used 
digital  FM  or  CPFSK,  would  permit  higher  hop  burst  rates  (R>)  without  incurring 
mutual  interference,  and  also  would  enhance  the  synchronization  performance,  with  the 
result  that  a  hopping  rate  four  times  higher  than  that  of  SINCGARS  is  practical.  In 
fact,  the  burst  rate  during  the  hop  dwell  times  in  the  JTIDS  and  EPLRS  systems  is 
5 MHz,  greater  than  the  hop  channel  spacing  of  3  MHz,  and  both  systems  employ  error 
control  coding2  and  a  form  of  direct  sequence  spreading-in  that  they  utilize  MSK  (digi¬ 
tal  FM  with  a  modulation  index  of  /i  =  0.5,  which  is  very  similar  to  the  bi-phase 
modulation  most  often  associated  with  direct  sequence  systems. 

1Some  analysts  categorize  hopping  speed  on  an  absolute  basis,  calling  hopping  faster  than  500  hops/sec 
fast  hopping. 

2For  Army  use,  neither  JTIDS  or  EPLRS  will  be  used  for  voice  communications.  This  allows  the 
system  designer  much  more  freedom  in  selecting  modulation  parameters. 

32 


Implementation  Studies 


3.1.2  Emissions  Control  and  Hop  Transient  Effects 

As  alluded  to  previously,  at  the  points  in  time  when  the  hopping  signal  switches 
from  one  carrier  frequency  to  another,  there  is  the  possibility  of  spectral  splatter  if  the 
amplitude  of  the  waveform  is  not  shaped  properly.  It  is  well  known  that  “sharp  edges” 
on  signal  pulses  give  them  higher  spectral  sidelobes  — and  therefore  higher  potential  for 
adjacent  channel  interference— than  do  “smooth  edges.”  For  this  reason,  both  EPLRS 
[32,  p.  7-4]  and  JTIDS  [33,  p.  23]  implementations,  for  example,  control  the  rise  and  fall 
times  of  the  hops.  In  addition,  the  JTIDS  transmissions  are  subject  to  strict  controls  on 
the  amount  of  power  emitted  during  the  hop  dwell  time  as  well  as  the  duty  cycle  of 
transmissions.  In  [29],  it  is  stated  that,  in  order  to  reduce  the  amount  of  splatter  into 
nearby  nets,  it  is  necessary  that  neither  the  rise  time  or  the  fall  time  of  the  hop  pulse  be 
much  less  than  a  bursted  symbol  time,  TC=R71,  thus  requiring  that 

min  (rise  time,  fall  time)  >  aTc,  (3-3 J 

where  the  factor  a  is  said  to  be  typically  between  the  values  1  and  4. 

A  separate  concern  related  to  emission  control  is  the  emission  of  RF  energy  in 
adjacent  and  other  hop  channels  due  to  the  characteristics  of  the  hopping  waveform 
during  the  transient  period  between  hops.  In  effect,  the  waveform  is  the  result  of  a 
frequency  modulation  of  a  carrier  at  the  center  of  the  hopping  band  by  a  multi-level 
pulse,  one  level  per  possible  hop  center  frequency,  and  the  spectrum  of  such  a  modula¬ 
tion  can  be  very  wide,  containing  many  spurious  harmonics  related  to  the  frequency 
switching.  For  this  reason.  VHF  hopping  radios  using  an  indirect  synthesizer,  such  as 
the  British  JAGUAR  [34,  pp.  6-10,  7-7],  employ  a  “hop  mute”  or  power-down  control 
function  during  the  time  that  the  synthesizer  circuitry  indicates  an  “out  of  tune" 
condition.  If  a  direct  frequency  synthesizer  is  used,  it  is  possible  to  adequately  suppress 
the  spurious  harmonics  while  maintaining  a  very  small  switching  time  [29],  in  which 
case  the  controlling  concept  in  waveshaping  is  the  prevention  of  splatter. 

Even  with  the  use  of  waveshaping  to  reduce  splatter,  there  may  be  undesirable 
effects  at  the  desired  receiver,  such  as  [29]  ringing  in  the  IF  filter  as  it  is  started  up  with 
the  new  hop  input,  as  well  as  phase  discontinuities  arising  from  Doppler  shifts,  fre¬ 
quency-selective  fading,  indeterminacy  of  the  IF  phase  prior  to  reception  of  the  pulse 
(assuming  the  signal  power  is  reduced  between  pulses),  and  imperfect  synchronization  in 
the  receiver. 

An  analysis  of  phase  transients  due  to  hopping  was  presented  in  [30]  for  the  case  of 
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a  MFSK  system  that  hops  once  per  symbol.  The  transients  were  shown  to  give  rise  to  a 
spreading  of  signal  energy  so  that  there  is  a  loss  of  energy  in  the  transmitted  symbol’s 
sub-band  and  undesired  energy  in  the  other  sub-bands.  For  hopping  systems  such  as 
VHF  tactical  radios  which  have  many  symbols  per  hop,  there  is  a  concern  that  the 
phase  transient  will  degrade  the  first  symbols  on  the  hop -or,  conversely,  that  if  there  is 
a  delay  in  transmitting  symbols  until  the  transient  settles  out,  then  an  effective  loss  of 
usable  signal  energy  is  incurred,  lowering  the  performance  of  the  system  for  a  given 
limit  to  transmitter  power  [1,  29]. 

3.1.3  Synchronization  Considerations 

Generally,  hopping  radio  systems  use  preambles,  special  sequences  of  frequencies 
and/or  data,  in  order  to  establish  synchronization  between  hopping  radios.  This  tech¬ 
nique  assumes  that  a  coarse  synchronization  has  already  been  achieved,  so  that  the 
receiver  knows  what  hop  frequencies  are  upcoming  but  not  the  exact  time  of  arrival, 
due  to  the  dependence  of  the  time  of  arrival  upon  the  relative  positions  of  transmitter 
and  receiver,  which  may  change  from  message  to  message. 

The  timing  references  in  hopping  radios  are  designed  to  Lave  an  accuracy  which  at 
least  guarantees  that  hopping  by  the  transmitter  and  receiver,  once  synchronized,  can 
maintain  synchronism  over  the  duration  of  a  typical  multi-hop  message.  Nevertheless, 
there  is  always  some  residual  difference  between  the  timing  of  the  incoming  signal  and 
that  of  the  receiver,  usually  neglected  in  analyses,  but  sometimes  modelled  for  the 
purpose  of  particular  studies  (e.g.,  [31]).  As  mentioned  above,  a  difference  in  switching 
times  between  the  incoming  signal  and  the  receiver’s  synthesizer  can  produce  significant 
phase  discontinuities.  *"*■ 

In  the  following  subsections,  some  of  the  hopping  transient  considerations  we  have 
discussed  are  given  further  attention  in  the  form  of  detailed  analysis. 

3.2  ANALYSIS  OF  HOP  TRANSIENT  EFFECTS 

In  [30]  it  has  been  shown  that  the  phase  transient  at  the  output  of  an  indirect, 
phase-lock-loop  based  synthesizer  can  be  modelled  by  the  waveform 

6,(t)  =  ^e-<“"'TJ=sin(w„jw‘2t),  (<1,  (3-4) 

fn  Jl^C2 

where  w„  =  2 nfn  is  a  loop  resonant  radial  frequency,  C  is  a  damping  factor,  and  A/ is  the 
desired  step  change  in  frequency.  Thus,  for  example,  we  observe  that  in  general  there  is 
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some  ringing  and  some  overshoot  in  the  phase  output.  In  this  section,  we  utilize  a 
simplified  model  for  the  phase  transient  in  order  to  analyze  the  spectrum  of  the  hopping 
transmissions  and  the  effects  of  phase  discontinuities  at  the  receiver  IF  when 
transmitter  and  receiver  are  not  in  perfect  synchronization. 

3.2.1  Model  for  Synthesizer  Transient  Waveforms 

For  the  purpose  of  studying  the  effects  of  the  transients  occurring  between  hops,  we 
neglect  the  data  modulation  and  any  hop  amplitude  control,  and  look  very  closely  at 
the  interval  between  two  hops,  taking  t  =  0  to  be  the  moment  at  which  the  frequency 
switching  begins.  Thus  we  write  for  the  received  signal  waveform 


S0(t)  =  Ar  COS  [0r(t)j  =  ArCOsfJ  dt^r(t)] 

(3-5a) 

where  the  instantaneous  frequency  or  phase  derivative  is  given  by 

MO/2*  =  fx. 

t  <  0 

(3-5b) 

=  fx  +  *tsign(/2-/j), 

0  <  t  <  t. 

(3— 5c) 

t  >t,. 

(3-5d) 

This  formulation  expresses  the  concept  that,  prior  to  the  end  of  the  hop,  the  carrier  has 
been  settled  at  fx  and  that  beginning  at  t  =  0  this  carrier  is  linearly  frequency-modulated 
for  an  interval  equal  to  t*.  the  switching  time,  after  which  the  carrier  is  /2.  It  is 
assumed  in  (3-5)  that  either  the  frequency  rate  of  change  6  is  given,  in  which  case  it 
turns  out  that  the  switching  time  is  the  frequency-dependent  quantity 

t,  =  |/2-/i|/6,  6  given,  (3-6a) 

or  the  switching  time  t..  is  given,  in  which  case  the  frequency  rate  of  change  must  be 
frequency-dependent : 

6  =  I/2-/1IA.,  t,  given.  (3-6b) 

3. 2. 1.1  Received  and  local  oscillator  transient  phase  functions 

Whether  6  is  specified  or  t,  is  specified,  using  F=f2- fx  the  transient  phase  function 
is  given  by  (assuming  0r(O)  =  0  for  convenience) 


$  r  ( t )  =  [ldu0r(u) 

Jo 

=  2*fxt, 

t  <0, 

(3— 7a) 

=  2;r/]t  +  irFt7 /t». 

0  <  t  <  t,, 

(3-7b) 

=  2 nfxt  -1-  2nF(l— 1>)  +  irFt,, 

t  >  t,. 

(3-7c) 
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Assuming  that  downconversion  is  employed  by  the  receiver,  if  the  transmitter  and 
the  receiver  are  in  perfect  synchronism  the  receiver  synthesizer  produces  the  local 
oscillator  frequency  fio  =  ffiF  +  fiF ’  w^ere  Irp  and  fJF  are  the  selected  (hopping)  channel 
carrier  frequency  and  the  fixed  IF  center  frequency,  respectively.3  However,  allowing 
the  receiver  to  be  late  by  the  time  interval  r,  the  receiver  synthesizer  waveform  is 
assumed  to  have  the  instantaneous  frequency  or  phase  derivative 

^LoW/2*  =  0r(t-r)/2jr  +  fJp.  (3-8) 

This  assumption  leads  to  the  local  oscillator  transient  phase  function  (neglecting  any 
constant  of  integration  which  may  differ  from  that  for  the  received  phase  function) 

=  +  6r(u—  r)] 

=  2llf]Ft  +  0r(t  — t)  +  2  Tr/jT,  (3-9) 

since  6r(  —  r)  =  —2 7r/jT. 

3. 2. 1.2  Mixer  output  transient  phase  difference 

The  mixer  output  then  is,  prior  to  the  IF  bandpass  filter, 

const,  x  cos  [0LO(t)-0r(t)]  =  Acos[2rr///rt  +  A0(t)]  (3-10a) 

where  the  transient  phase  difference  A0  is 

A0(t)  =  flr(t-r)  -  flr(t)  +2rr fxr.  (3-10b) 

The  form  that  this  function  takes  depends  on  the  relation  between  t»,  the  synthesizer 
switching  time,  and  r,  the  receiver  synchronization  offset.  For  indirect  synthesizers 
-and  indeed  most  synthesizers -it  can  be  assumed  that  |r|  <t*.  In  that  case  and  in  the 
case  of  r  >0,  the  transient  phase  difference  is  found  to  be 


o' 

II 

<1 

t<0 

(3-lla) 

=  -*Ft2/t„ 

0  <  t  <  r 

(3-1 lb) 

=  — jrF(2rt  — r2)/t,, 

T  <  t  <  t. 

(3-llc) 

=  —  2nFr  +  nF(t  —  r  —  t,)2/t», 

t,  <t<t,  +  r 

(3- 1 Id ) 

=  —  2irFY, 

t  >  t#  +  T. 

(3—1 le) 

Note  that  if  r= 0,  that  is,  there  is  no  synchronization  error,  then  there  is 

no  transient 

phase  difference  and  the  combination  of  hopping  and  de-hopping  is  indeed  “transparent’’ 


3Using  this  method  it  is  necessary  to  use  additional  pre-filtering  at  RF.  such  as  done  by  the  Jaguar 
radio  [34],  to  prevent  images  of  non-selected  channels  from  appearing  at  the  mixer  output. 
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as  assumed  by  most  simple  analyses  of  hopping  systems.  Also,  we  observe  that  (3-11) 
describes  AS  as  a  piecewise  continuous  function  with  a  succession  of  linear,  quadratic, 
linear,  quadratic,  and  linear  line  segments,  as  illustrated  in  Figure  3-1. 

For  the  case  that  r  >  t,  >  0,  the  transient  phase  difference  takes  the  functional  form 


A0(t)  =  0,  t<0  (3-12a) 

=  -*Ft  2/t,,  0<t<t.  (3-12b) 

=  -2irF(t-t.)-*Ft.,  t,  <t<r  (3-12c) 

=  -2itFt  +  irF(t  —  r  —  tj)2/t,,  r  <  t  <  r +  t«  (3-12d) 

=  -2 nFr,  t  >  tj  +  r.  (3-12e) 


As  illustrated  in  Figure  3-1,  this  functional  form  has  the  same  structure  as  the  one  for 
the  case  of  r  <  t,. 

A  more  convenient  model  of  the  transient  phase  difference  for  analytical  purposes  is 
easily  found  to  be 

A0(t)  =  0,  t  <0  (3-13a) 

=  —  it  Ft  +  jr  Ft  COS  ^  )  •  0<t<t,  +  r  (3-13b) 

=  -2 ttFt.  t  >  t»  +  r.  (3-13c) 

Note  that  this  shifted  cosine  function,  like  (3-11)  and  (3-12),  also  has  quadratic  seg¬ 
ments  at  endpoints  and  a  linear  segment  in  the  middle,  to  an  excellent  approximation, 
but  is  continuous  on  the  interval  (0,  t.«  +  r)  rather  than  piecewise  continuous. 

3.2.2  Effect  of  Receiver  Filtering  on  Transients 

The  harmonic  content  of  the  transient  phase  difference  can  be  discovered  by  con¬ 
sidering  expansions  of  the  mixer  output  signal  waveform  prior  to  the  IF  bandpass  filter: 

s(t)  =  A  cos  [27r//ft  +  Afl(t)].  (3-14) 

If  the  functional  form  for  A0(t)  given  by  (3-11)  or  (3-12)  is  used,  then  the  analysis 
involves  expansion  of  terms  such  as 

cos  (at2)  =  j  dzcos(ax2) 
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where  C(z)  is  the  Fresnel  cosine  integral  [17,  eq.  7.3.1].  The  mathematics  involved  is 
similar  to  that  for  analysis  of  “chirp  radar”  systems  [35].  Note  that  the  harmonic 
structure  is  not  yet  apparent  in  (3-15),  but  requires  further  development.  However,  if 
instead  of  the  models  for  the  transient  phase  difference  given  by  (3-11)  and  (3-12)  we 
use  the  shifted-cosine  one  given  by  (3-13),  we  have  terms  to  evaluate  such  as 

cos(acosit)  =  J0(a)  +  2^  (-1)"  J2n(a)cos(2n£>t)  (3-16a) 

oc  n  =  1 

=  Y^{-\)nen3in{a)cos{2nbt)  (3-16b) 

,  n=0 

and 

sin(acosk)  =  2^  (-1)"  J2n+1(a)cos[(2n+l)6t],  (3-16b) 

n  =  0 

in  which  the  J«(i)  axe  the  Bessel  functions  of  integer  order  [17,  ch.  9]  and  („  is  Euler’s 
constant,  which  takes  the  values  c0  =  l  and  f„  =  2  for  n>0.  Note  that  using  (3-13)  for 
the  transient  phase  difference  will  therefore  give  the  harmonics  directly. 

Applying  (3-16)  to  the  expression  for  the  mixer  output  gives  the  expression 

s(t)  =  Acos(27r/;/-t),  t<0  (3-17a) 

=  Acos(2,/,ft-,FT)£(-l)"<„J,„('FOcos(2aii) 

n  =  0 

+2  A  sin  (2»/,ft  -  »Fr)  y;  ( - 1 )"  r)  cos 

0  <  t  <  ta  +  r  (3-17b) 

=  A  cos  (27r//ft— 2?r  Ft),  t  >  t»  +  r.  (3-17c) 

Since  it  is  true  that  Jo(0)  =  1  and  Jn(0)  =0  for  n>0,  we  see  that  this  somewhat  complex 
expression  reduces  to  simply  s(t)  =  Acos(2jr//ft),  that  is,  no  transient  phase  difference 
at  all,  for  the  case  of  perfect  synchronism  (r  =  0). 

From  (3-17)  we  see  immediately  that  during  the  interval  (0,  t,  +  r)  the  mixer  out¬ 
put’s  energy  is  spread  about  the  IF  center  frequency  to  the  harmonics  of 

he  =  |(t.+r)-1.  (3-18a) 

The  distribution  of  energy  among  these  harmonics  is  dependent  on  *Ft,  the  argument 
of  the  Bessel  functions.  Assuming  a  synch  timing  offset  of  about  1%  of  the  symbol 
period,  Tc,  gives  a  value  for  r  of  about  0.5/isec.  For  a  typical  VHF  hopping  radio,  the 
size  of  F  can  be  from  25  kHz  to  over  50  MHz,  but  more  likely  would  be  from  25  kHz  to 
around  256  x  25  kHz  =  6.4  MHz  because  the  radios  often  are  set  up  to  hop  in  a  sub-band, 

rather  than  over  the  entire  30-88 MHz  VHF  band.  Thus  the  value  of  nFr  ranges  from 

around  to  more  than  3.2  s-.  For  small  values  of  the  argument,  the  Bessel  functions 
are  well  approximated  by 
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J„(*Fr)  «  ,Ft<1. 

(3- 18b) 

For  large  values  of  the  argument,  we  have  [17,  eq.  9.2.1] 

Jj.(*fr)  «  (4^  ^cos(,F,-5) 

(3-18c) 

and 

J2n+i(’rF’r)  a  (  J  J^sin(irFr  =) . 

(3-18d) 

From  (3-18c)  and  (3-18d),  we  observe  that  the  harmonics  for  large  values  of  vFr  are  of 
approximately  equal  magnitude,  which  agrees  with  the  flat  spectra  for  chirp  waveforms 
shown  in  [35]. 


It  is  clear  that  the  transient  phase  which  will  appear  at  the  input  to  the  de¬ 
modulator  will  be  considerably  different  from  that  at  the  mixer  output,  since  the  IF 
bandpass  filter  will  attenuate  these  harmonics  to  different  degrees.  For  deriving  the 
effect  of  the  IF  bandpass  filter  on  the  transient  phase  difference  when  it  is  assumed  that 
a  Gaussian-shaped  filter  is  used,  a  detailed  transient  analysis  is  presented  in  Appendix 
C.  The  implications  of  that  analysis  are  that  the  IF  filter  will  indeed  reject  the 
harmonics  of  fA6.  Therefore  in  anticipation  of  the  rejection  of  the  harmonics,  we  make 
use  the  unit  step  function  U(t)  to  reformulate  (3-17)  as  follows: 

s(t)  =  Acos(21r/;ft)-U(-t)  +  A  J 0(nFr) cos (2vf[Ft-nFr)  •  [U(t)-U(t-t»-r)] 

+  Acos(2nfIFt-2nFr)  -  (3-19) 

The  corresponding  in-phase  and  quadrature  components  are 

I(t)  =  A  -U(-t)  +  AJ0(irFr)cos(irFr).[U(t)-U(t-tJ-r)] 

+  Acos(2?rFr)  •  U(t  —  tj  —  t)  (3-20a) 

and 

Q(t)  =  AJ0(7rFr)sin(7rFr).[U(t)-U(t-t.-r)] 

+  Asin(2irFr).U(t-t.-r).  (3-20b) 

To  characterize  the  IF  filter's  response  to  these  inputs,  we  use  the  notation  S(t)  to 
denote  the  equivalent  lowpass  filter’s  response  to  the  unit  step  function,  that  is. 

S(t)  =  [  dzh0(t-x)  (3-21a) 

J0 

=  [  dzh0(x)  =  PG[2fiW/f(t-t<<)],  (3-21b) 

J0 

where  PG(-)  denotes  the  Gaussian  cumulative  probability  distribution  function  and  td  is 
the  filter  delay.  The  quadrature  components  at  the  output  of  the  IF  filter  then  are 
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v(t )  =  A-S(-t)  +  AJ0(5rFr)cos(7rFr).[S(t)-S(t-t.-r)] 

-f  Acos(2?r.Fr)  •  S(t  — 1«  — r)  (3-22a) 

and 

u(t)  =  A  J0(?r.fY)  sin  (fffY)  •  [S(t )  —  S(t  —  t,  —  r)] 

+  Asin(2r  Fr)-S(t-t.-r).  (3-22b) 

Assuming  that  the  filter  delay  is  approximately  one  symbol  duration,  that  is,  about 
50/isec,  and  that  the  combined  switching  time  and  synch  offset  t»  +  r  are  much  less  than 
this  delay,  say  about  10/isec,  the  filtered  transient  amplitude  and  phase  functions  a,(t) 
find  4>«(t)  at  the  output  of  the  IF  filter,  given  by 

a,(t)  =  [u2(t)+v2(t)]1/2  (3-23a) 

and 

0,(t)  =  tan~‘[u(t)/v(t)],  (3-23b) 

are  as  sketched  in  Figure  3-2  for  the  assumption  of  a  large  value  of  irFr.  The  cases  con¬ 
sidered  in  that  figure  are 

ttFt  =  m,  m± i,  m±\,  (3-24) 

where  m  is  an  integer.  It  is  evident  that  the  transition  event  takes  place  over  the  time 
interval  defined  by 

|t-td  — |(t,  +  r)|  <  3<r„  =  *  Tc,  (3-25) 

assuming  that  W/fTc  =  1,  and  that  in  general  the  signal  is  subject  to  an  attenuation  or 
energy  loss  during  this  interval  when  the  phase  is  switching.  For  this  reason,  it  is  to  be 
expected  that  the  system  noise  at  the  IF  filter  output  will  largely  determine  the  value  of 
phase  during  this  transient  interval.  The  critical  case  is  seen  to  be  it  Ft  equal  to  an 
integer  plus  since  for  this  case  the  signal  amplitude  decreases  momentarily  to  zero 
and  the  (modulo-2?r)  phase  changes  rapidly  by  ±ir,  even  in  the  absence  of  noise.  It  is 
highly  likely  that  an  FM  noise  click  or  rapid  2?r  phase  shift  will  occur  in  this  case  also. 

The  conclusion  to  be  drawn  from  this  analysis  is  that  the  effects  of  the  transient 
phase  difference  between  the  incoming  hopping  waveform  and  the  receiver’s  synthesizer 
are  primarily  confined  to  the  interval  (0,  Te+t.  +  r),  referenced  to  the  input  of  the 
receiver  with  t  =  0  at  the  instant  that  the  incoming  signal  begins  to  switch  frequencies. 
Therefore,  a  one-symbol  guard  time  is  advisable  before  resuming  the  transmission  of 
data.  Typically,  as  has  been  mentioned,  the  amplitude  shaping  and  other  procedures 
designed  to  control  emissions  dictate  a  larger  off  time,  and  so  are  the  controlling  factor 
in  the  determination  of  what  fraction  of  the  dwell  time  may  be  used  to  transmit  data. 
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Key:  A:  Ft  =  integer 

B:  Ft  =  integer  +  i 

C:  Ft  =  integer  +  * 

D:  Ft  =  integer  -  i 

E:  Ft  =  integer  -  \ 


FIGURE  3-2  IF  SIGNAL  AMPLITUDE  AND  PHASE  TRANSIENT  WAVEFORMS 
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4.0  INVESTIGATION  OF  ECCM  COMBINING  TECHNIQUES 

In  this  section  we  summarize  investigations  of  the  effectiveness  of  various  ad  hoc 
diversity  combining  techniques  for  improving  the  uncoded  performance  of  FH/CPFSK 
in  worst-case  partial-band  noise  jamming.  Techniques  used  in  conjunction  with  the 
limiter-discriminator  receiver  axe  discussed  in  Section  4.2,  and  those  pertaining  to  the 
differential  detector  are  treated  in  Section  4.3. 

4.1  BACKGROUND:  NONCOHERENT  COMBINING  LOSS 

When  the  FH/CPFSK  signal  is  accompanied  by  stationary  Gaussian  noise  only 
(either  background  noise  or  continuously  present  noise  jamming),  the  use  of  multihop 
per  symbol  diversity  is  not  called  for.  except  perhaps  to  reduce  the  interceptability  of 
the  signal  by  spreading  the  signal  energy  in  frequency.  As  illustrated  in  Figure  4-1,  for 
noise-only  conditions  the  unjammed  probability  of  bit  error  increases  with  L ,  the  num¬ 
ber  of  times  the  bit  is  split  into  chips  which  are  transmitted  on  different  hops,  with  the 
bit  recovered  by  examining  the  sign  of  the  sum  of  the  demodulated  differential  phase 
samples.  This  loss  in  performance  as  L  increases  is  commonly  termed  “noncoherent 
combining  loss"  (NCL). 

4.1.1  Diversity  Sum  Error  without  FM  Noise  Clicks 

An  understanding  of  the  mechanisms  which  give  rise  to  NCL  for  FH/CPFSK  sys¬ 
tems  can  be  gained  by  considering  a  simplified  analysis  of  differential  phase  combining. 
If  FM  noise  clicks  are  ignored,  the  differential  phase  samples  for  L  hops/symbol 

zi<t  =  A4>(ti?),  g=l,  2,  ...,  L  (4-1) 

for  the  A:th  symbol  are  approximately  Gaussian  random  variables,  with  identical  means 
=  A<t>  and  variances  al  =  Llp.  That  is,  such  is  the  case  when  the  symbol  energy  is 
kept  constant  as  L  changes,  so  that  the  effective  SNR  when  the  samples  are  taken  is 
pj L.  The  mean  value  is  the  (distorted)  differential  phase  in  the  absence  of  noise. 

Using  then  the  approximation  that  the  samples  are  Gaussian  random  variables,  the 
sum  of  the  L  differential  phases  is  approximately  Gaussian  with  mean  Lpt  =  LA<j>  and 
variance  Lcrl  =  X2/p,  denoted  by 

L 

z*  = 

»=> 

*=  G(LAo,  L2 / p)  —  L-G(A4>,  1  /p).  (4-2) 
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In  (4-2)  we  use  G(p,  a2)  to  represent  a  Gaussian  random  variable  with  mean  p  and 
variance  <r2,  and  in  the  last  line  of  (4-2)  we  make  use  of  the  fact  that  multiplication  of 
such  a  random  variable  by  a  constant  K  gives  K-G(/i,  <r2)  =  G(Kp,  K2<r2). 

Now,  since  the  symbol  decision  is  made  of  the  basis  of  the  sign  of  zk,  (4-2)  implies 
that  in  the  absence  of  FM  clicks,  the  symbol  error  probability  for  any  value  of  L  is 
approximately  given  by 

P.(L)  =  Pr{zt< 0 1  dk=  1} 

«  Pr{T-G(Ad>,  1/p)  <  0} 

=  Pr{ G(A<£,  1/p)  <0},  (4-3) 

which  does  not  depend  on  the  value  of  L.  That  is,  there  is  no  combining  loss. 

4.1.2  Effects  of  FM  Noise  Clicks 

However,  when  the  clicks  are  included  in  the  analysis,  by  raising  the  characteristic 
function  for  one  of  the  A4>  samples1  to  the  Lth.  power  we  can  deduce  that  the  sum  of 
the  differential  phase  samples  is  conditionally  Gaussian,  given  the  value  of  Ne.  the 
number  of  clicks.  The  average  number  of  clicks  in  the  sum  is  L  times 
aL=e~p/L\Ad>\ /2ir,  the  average  number  of  dicks  for  one  hop,  and  if  Nc=n ,  the 
Gaussian  variable's  mean  and  variance  are 

P n  =  LAo  -  sgn(A<?)-27r n,  a\  —  L2!p.  (4-4) 

Thus  with  clicks  included  in  the  analysis,  the  error  probability  is 

P t(L)  *  £  1 h2jJle-LoL  pr(G(A^-^H,  1/p)  <0}  (4-5a) 

ZZb  n-  " 

=  £  Qc[(i*-2p)^l,  (4-5b) 

n  =  0  11  ’  u 

which  clearly  does  depend  on  the  value  of  L.  In  (4-5b),  Qc(  •)  is  the  Gaussian  comple¬ 
mentary  cumulative  distribution  function. 

To  provide  a  numerical  example,  let  A <p  be  positive.  Then  we  have  for  L  =  1  and 
for  L  —  2  the  following  results: 

P«(l)  =  e"01Qc(A^^)  +  aie'aiQG[(A<i»-2T)^]  +  ... 

»:  Qie_t>1  -♦  qj  for  high  SNR; 


1The  characteristic  function  shown  in  (2-52a)  may  t>e  used,  with  7  =  0  and  q0  =  q£. 
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and  Pe(2)  =  e  2°2QG(A<5^p)  +  2a2e  °2Qg[(A<&  —  7r) ^p]  +  ••• 

«  2a2e'2a2  -  2o2  for  high  SNR.  (4-6b) 

Thus  for  high  SNR  the  ratio  of  P«(2)  to  P«(l)  is  about  equal  to 

Pe(2)/P«(l)  «  2a2/ai  =  2er/2.  (4-6c) 

Because  of  the  amplitude  distortion  a(t),  the  SNR  is  less  than  Ek\>f0;  a  typical  value  is 
p  =  0.7Et\jro.  For  E6\Xo  =  10dB,  then,  (4-6c)  gives  a  ratio  of  66.2.  In  Figure  4-1,  the 
ratio  of  error  probabilities  for  L  =  1  and  L  =  2  when  E6\X0  =  lOdB  is  0.0133/0.0002 
=  66.5.  This  agreement  confirms  that  the  so-called  noncoherent  combining  loss  for 
limiter-discriminator  detection  of  CPFSK  is  attributable  to  the  FM  noise  clicks. 

When  the  FH/CPFSK  system  is  subject  to  full-band  noise  jamming,  the  NCL 
shows  up  as  different  asymptotes  (minimum  error)  for  fixed  Eb\N0  as  jamming  power  is 
decreased.  For  example,  in  Figure  4-2  we  show  the  effect  of  full-band  jamming  on  the 
error  probability  for  an  L  hop/bit  FH/CPFSK  system  with  linear  (sum)  combining  of 
the  differential  phase  chips  when  Et\R0  is  selected  to  give  a  1(T5  Pc  for  L=  1. 


4.2.  HARD-DECISION  COMBINING  RECEIVER 


For  the  hard-decision  (HD)  receiver  for  FH/CPFSK.  a  polarity  (±1)  decision  is 
made  as  each  chip  is  received. -equivalent  to  estimating  the  original  chip  sequence  dkl. 
by  d^-and  the  bit  or  symbol  decision  after  all  L  chips  have  been  received  is  according 
to  the  rule 


(4-7) 


4.2.1  Analysis  of  the  HD  Receiver 

A  correct  symbol  decision  is  made  if  a  majority  of  the  chip  decisions  are  correct2: 
in  the  case  of  a  tie,  which  is  only  possible  if  L  is  even-valued,  a  random  binary  decision 
is  made.  Therefore,  the  probability  of  error  for  a  hard  decision  combining  receiver  is 

(L- l)/2 

P <(L;  RAr,  R ,,  y)  =  £  (£)  (l-p,)*pf-*,  L  odd  (4-Sa) 

k—0  XK/ 

(L- 2)/2 

=  J2  (%)  (i-pi)*pf'"‘  +  s(^2)^1-p,^p,^/2,  L  even  ^4_sb^ 

where 

p,  SP,(1;JR1,..  JRj.  ,)  (4-8c) 

and  where  we  have  used  the  notation  RA.=  E6/X0  and  ’RJ  =  Ehj H j. 

2Note  that  this  rule  is  equivalent  to  using  the  sign  of  the  median  of  the  L  chips,  for  L  odd. 
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By  reversing  the  order  of  these  sums,  the  functional  form  for  P t(L)  can  also  be 
written 

2*71  +  1  - 

P ,{L)  =  T  for  L=2m+1  (odd)  (4-9a) 

i  =  m  +  l'  *  ' 

for  L=2m+2  (even).  (4-9b) 


By  substituting 


/ 2tti-\-2\  /2tti4"1\  ,  /2m+l 

V  k  )  ~  \  k  )  +  l  Jfc-1 

into  the  the  equation  for  L  even,  and  noting  that 

(2mk+1)  =  0  for  k>  2m  +  1, 


) 


(4-10a) 


(4-10b) 


it  is  straightforward  to  establish  that  the  functional  form  of  Pe  as  a  function  of  p,  for  L 
even  is  the  same  as  for  the  next  lower  odd  value  of  L.  For  example,  functionally. 


P«(2:  Pi)  =  Pe(l:  Pi)  =  Pi.  (4-lla) 

P«(4;  Pi)  =  Pe(3;  p,)  =  3p?  —  2pf .  etc.  (4-llb) 


A  general  expression  can  be  found  by  noting  that  [17,  §6.6] 


P((2m+1;  pj 


(?™+*I!  [Pldi  (l-x)mim 

m!  mi}.  J  o 

(2m  +  l)!  ^  (m\  (-1)"  m  +  n  +  1 

ml  ml  ^{nJm+n+l^ 


(4-12) 


However,  it  should  be  noted  that  pj  itself  is  a  function  of  L,  as  stated  in  its  definition 
above.  The  practical  application  of  this  observation  about  the  functional  form  of  Pe  for 
the  HD  receiver  is  that  in  programming  any  calculations  of  Pe(L),  we  may  consider  only 
odd  values  of  L.  The  error  probability  for  even  values  of  L  can  be  obtained  after 
computing  P«(£-l)  simply  by  replacing  the  ratios  R^  and  R^  with  RA,  x(L-l)/L  and 
RJx(T-l)/T  and  then  reusing  the  subroutine  for  Pe(L-l)  to  compute  P«(T). 

In  [8]  it  was  shown  that  a  precise  expression  for  p!  is 

Pi  =  ((1—  ■) )  +  7  F{KT/L)  )isi  patterns  (4— 13a) 

where  for  a  particular  ISI  pattern  (indexed  by  the  parameters  Ac>,  U,  V,  and  W)  and 
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using  a  limiter-discriminator  receiver. 


F(  R)  =  1  —  e~“  +  e 


a  ,  -o  rsm 


in  Ail’'7  ,  eXPft 


U— V  sinz+Wcoszl 
1  +  rcos  cos  i  J 
1-l-rcos  A^cosx 


...  fU  —  Vsinz- 

+  e-  WsinA^r 

4?r  J_,/2  1-n 


U— Vsinz-W  cos  A^coszi 


-rcosz 

rcosx 


(4-13b) 


An  accurate  approximation  for  modulation  and  receiver  parameters  h=0.7 ,  W,fT  =  1.0, 
and  a  Gaussian  IF  filter,  is  the  following: 


(F(  R)> 


jexp{  -.974717R}, 


R  <  0.5 


.346664exp{  -  .705883  R  +  .  115918  R-1}  R>0.5. 


(4-14) 


Figure  4-3  shows  a  comparison  of  this  approximation  with  the  exact  unjammed 
error  probability.  Note  that  a  linear  scale  is  used  for  the  SNR. 


4.2.2  Results  Using  the  HD  Receiver  with  Diversity 

Using  the  approximation  given  above  for  the  chip  error  probability  as  a  function  of 
SNR,  -we  calculated  the  uncoded  performance  of  hard- decision  combining  of  X-hop/bit 
FH/CPFSK  in  worst-case  partial-band  noise  jamming.  The  worst-case  result  was  found 
by  maximizing  P e(L;  •>)  with  respect  to  *>.  the  partial-band  jamming  fraction,  for  given 
values  of  RjV  and  R^. 

In  Figures  4-4  and  4-5.  for  E6/Jf0=20dB,  we  show  the  worst-case  P,  for  L= 1  to  7 
and  for  L-l  to  12,  respectively.  As  noted  above,  characteristically  for  the  hard-decision 
receiver,  Pe(T=2m)  >  Pe(T=2m— 1),  that  is.  the  error  is  higher  for  L= 2  than  for  L—  1. 
higher  for  L— 4  than  for  L= 3,  etc.  This  behavior  is  attributable  to  the  necessity  of 
making  a  random  decision  in  the  case  of  a  tie  when  L  is  even-valued.  Except  for  this 
predictable  variation,  the  trend  is  that  the  value  of  Ry  required  for  a  10'5  uncoded  bit 
error  probability  decreases  as  L  increases,  indicating  a  diversity  gain. 

For  clarity,  the  error  curves  for  only  the  odd  values  of  L  axe  plotted  in  Figure  4-6. 
Here  the  facts  that  there  is  a  diversity  gain  and  an  optimum  value  of  diversity  are  more 
clearly  observable.  Not  that  the  error  curves  are  “inverse  linear”  for  the  given  scale, 
with  slopes  equal  to  -\{L+ 1),  revealing  that  the  worst-case  error  is  inversely 
proportional  to  the  i(L  +  l)th  power  of  the  signal  enerv-to-jamming  noise  density,  R^. 
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FIGURE  4-3  EXACT  AND  APPROXIMATE  CPFSK  ERROR  PROBABILITY 
COMPARISON  FOR  LIMITER-DISCRIMINATOR  DETECTION 


~-5  5  15  25  35  45 

E t/Xj  (dB) 

FIGURE  4-4  HARD-DECISION  DIVERSITY  PERFORMANCE  FOR  L=  1  TO 
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E JKj  (dB) 

FIGURE  4-5  HARD-DECISION  DIVERSITY  PERFORMANCE  FOR  L=1  TO  12 


52 


E JXj  (dB) 

FIGURE  4-6  HARD-DECISION  DIVERSITY  PERFORMANCE  FOR  ODD  L 
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If  Eb/Jf0  is  reduced  (the  thermal  and  system  noise  is  increased  relative  to  the  signal 
energy),  the  noncoherent  combining  losses  (NCL)  will  affect  the  diversity  gains  against 
the  jamming.  In  Figure  4-7,  for  Ei>/-N‘0  =  15dB,  we  observe  that  the  maximum  value  of 
optimum  diversity  is  L—l.  for  a  narrow  range  of  Rj  around  the  value  Rj«20dB,  and  in 
Figure  4-8,  for  Ej/>f0  =  13dB,  the  highest  value  of  diversity  for  which  improvement  is 
noted  is  L— 3.  It  is  apparent  also  from  Figure  4-8  that,  as  the  thermal  noise  increases, 
the  NCL  for  hard-decision  combining  is  significantly  greater  for  even  values  of  L  than  it 
is  for  odd  values. 


4.3  LIMITER-DISCRIMINATOR  TECHNIQUES 

In  this  subsection  we  summarize  our  analysis  and  present  numerical  results  con¬ 
nected  with  studies  of  possible  ECCM  weighting  schemes  to  enhance  the  performance  of 
a  limiter-discriminator-based  FH/CPFSK  receiver  against  partial-band  noise  jamming. 


4.3.1  Approximations  to  the  Differential  Phase  PDF  and  CDF 

Analytically  it  is  quite  convenient  to  utilize  the  probability  of  error  expression  in¬ 
volving  the  characteristic  function  (CHF)  of  the  diversity  sum  decision  variable,  z: 


Pe 


i  i 

5  TT 


(4-15) 


The  convenience  associated  with  using  this  form  arises  from  the  fact  that  the  CHF  for 
the  sum  of  L  independent,  identically  distributed  random  variables  is  simply  that  for 
one  variable,  raised  to  the  Xth  power: 

Cz{i')  =  [Clq{v)]L  (4-1 6a ) 

i 

when  z  =  £  z?.  (4-16b) 

i=i 

A  general  expression  for  Ctg{ u)  was  given  in  Section  2.3,  including  FM  noise  clicks 
and  partial-band  noise  jamming.  To  use  that  expression,  it  is  necessary  to  formulate 
the  CHF  corresponding  to  p^(x),  the  probability  density  function  (PDF)  for  a  single  2 
interval  of  possible  A4>  sample  values.  Thus  we  require  an  expression  for 


£*(•')=[  dxe5"1  p„(z).  (4-17) 

J  -00 

Unfortunately,  an  exact  analytical  expression  for  this  CHF  cannot  be  found  in  closed 
form,  since  p^(z)  itself,  given  previously  in  (2-29),  is  in  integral  form.  Therefore  we 
have  investigated  several  approximate  expressions. 
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o 


-5  5  25  25  35  45 

E»/.n -j  (dB) 

FIGURE  4-8  HARD-DECISION  DIVERSITY  PERFORMANCE  FOR  Et/x0  =  13dB 
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4. 3. 1.1  Expansion  of  the  CHF  in  terms  of  cumulants 

The  CHF  for  a  random  variable  x  can  be  expressed  in  terms  of  its  moments  or  in 
terms  of  its  distribution’s  semi-invariants  or  cumulants.  The  general  expression  is 


Cr{u) 


r  =  0 


or 

r! 


mr 


(4-18) 

(4-19) 


in  which  the  quantities  { mr }  are  the  (noncentral)  moments  and  the  {  «„  }  are  the  semi¬ 
invariants  or  cumulants  [17].  For  example,  in  terms  of  the  moments,  some  of  the 


cumulants  can  be  shown  to  be 

K!=m,  (4-20a) 

Kj=uJ  =  m2-m]  (4-20b) 

*3=  m3-3m1m2  -  mf  (4-20c) 

k4=  m4  —  4  k^'3—  3  k\  —  6  k\k2—  nt  ■  (4 — 20d ) 


For  the  Gaussian  distribution,  all  the  cumulants  higher  than  and  k2  (the  mean  and 
variance)  are  zero. 

Since  for  high  SNR  the  PDF  for  the  modulo-27r  differential  phase  is  approximately 
Gaussian  [13.  IS],  we  reasoned  that  the  CHF  C^.(u)  could  be  well-approximated  by  a 
truncated  version  of  (4-19), 

C^{v)  %  ^expj  Y  “p  Kn  j,  (4-21a) 


using  exact  cumulant  values  obtained  from  calculations  of  the  moments  about  Ac>: 

mr  =  E{(v'-  A«i)r}  =  |  dx  xrp„,(x+A<£),  (4-21b) 

using  the  expression  given  in  (2-29)  for  the  PDF  of  \p. 

For  zero  SNR  (noise  only),  p^x)  is  the  uniform  PDF,  with  zero-valued  odd 
moments  and  even  moments  given  by 

m2r=2rT I  (4-22a) 

which  give  rise  to  zero- valued  odd  cumulants  («i=  k3=  k5=  0)  and 

k2=  i,r2=  3.29  (4-22b) 

*<  =  -£v4=  -13.0  (4-22c) 

«6  =  —  If  jr«=  —  135S.  (4-22d) 


It  is  evident  from  (4-22)  that  for  sufficiently  low  SNR,  the  higher  cumulants  of  the 
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differential  phase  PDF  do  not  diminish  quickly  enough  to  permit  representation  of  the 
CHF  with  a  reasonably  small  number  of  cumulants.  Therefore,  this  approach  to  ap¬ 
proximating  the  CHF  was  abandoned. 

4. 3. 1.2  CHF  based  on  approximation  of  the  PDF 

For  high  SNR,  it  is  well  known  [18]  that  the  PDF  of  the  modulo- 2  tt  differential 
phase  i/>  is  well  approximated  by  a  Gaussian  PDF.  In  Figure  4-9,  we  show  a  direct 
comparison  of  the  exact  PDF  pv,(z)  with  the  PDF  for  a  Gaussian  random  variable  with 
the  same  mean  and  variance.  As  seen  in  the  figure,  the  approximation  is  excellent  for 
high  SNR  but  for  low  SNR  begins  to  depart  significantly  from  the  exact  PDF,  especially 
for  ip  values  near  to  A4>  ±  tt.  Moreover,  the  region  over  which  the  Gaussian  PDF  has 
significant  value  increases  with  the  value  of  the  variance,  whereas  the  exact  PDF  has 
the  finite  domain  |V’  —  A4>|  <  n;  therefore,  use  of  the  Gaussian  characteristic  function 

CG{u)  =  ei'“'-*V/2  (4-23) 

in  the  error  probability  calculation  is  not  expected  to  approximate  the  the  system 
performance  faithfully  for  low  SNR. 

For  this  reason,  we  have  developed  an  approximation  to  the  differentia]  phase 
characteristic  function  starting  with  a  “truncated  Gaussian”  PDF  that  we  define  by 

(  Ke  /2°°  for  |  a: |  <  7r 

prc(x)  i  |  (4-24a) 

l  0  elsewhere. 

The  properties  of  this  distribution  are  discussed  in  Appendix  B.  For  example,  the  nor¬ 
malization  factor  K  in  (4-24a)  is  given  by 

K_1  =  cr0f2n  erf  ^7r/<r0j2j,  (4-24b) 

and  it  can  be  shown  that,  as  the  parameter  <r0  -»  oo,  the  PDF  pTG(x)  becomes  a  uniform 
PDF  on  the  interval  (-tt,  tt).  Also,  the  variance  of  the  distribution  <t2tg  is  related  to  o\ 
by  the  expression 

r  _  2/2  2 

<r\G  =  (Tq  1  -  42*  e  *  ',°y/(T0erf  .  (4-24c) 

The  characteristic  function  for  the  truncated  Gaussian  distribution  may  be  approxi¬ 
mated  by  the  “mixture”  characteristic  function 

CM{ v)  =  fsinc(i/)  +  (l-t)e  <’0>'  /2,  (4-25a) 
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where  the  mixture  parameter  e  is  SNR-dependent  and  is  given  by 

2  2 

t  =  4*  e  '  ‘Jujjr042erf  ^/cr042^J  ;  (4-25b) 

also,  sinc(i')  =  sin(7u/)/5rt/  is  the  CHF  for  the  uniform  distribution  over  the  interval  (  —  ?r, 
*).  The  first  and  second  terms  of  CM(v)  correspond  to  low-SNR  (uniform)  and  high- 
SNR  (Gaussian)  versions  of  the  differential  phase  CHF,  that  is,  w  may  be 
interpreted  as 

C„M  =  +  [1  (4-25C) 

with  f  decreasing  as  the  SNR,  p,  increases. 

In  effect,  the  use  of  (4-25)  approximates  the  differential  phase  PDF  p^(x)  by  a 
mixture  of  of  PDF's.  Figure  4-10  illustrates  the  accuracy  of  this  PDF  approximation 
when  <tq  is  taken  to  be  the  inverse  function  of  the  truncated  Gaussian  variance  a2  given 
by  (4-24c).  Evidently,  for  low  SNR  this  value  of  crl  is  too  high,  judging  from  the  curves 
in  Figure  4-10  when  Ei,/^o=0dB. 


Taking  the  nfxture  approximation  concept  a  step  further,  we  found  by  trial  and 
error  that  the  differential  phase  PDF  is  well  approximated  by  the  simple,  heuristic  ap¬ 
proximation  based  on  (4-2oa)  with  c  and  o20  given  by 


f 


,-2SNR 


(4-26a) 


and 


2  <t2-V(/3 

■  r-< 


(4-26b) 


The  formula  for  a20  is  based  on  matching  the  actual  variance,  a2 .  with  that  of  the 
approximation,  (1  -c)al  +  x2( /"S.  In  Figure  4-11.  we  observe  that  this  simplified  mixture 
approximation  gives  a  better  fit  to  the  actual  PDF. 

In  Figure  4-12  we  compare  the  unjammed  CPFSK  probability  of  error  for  L=  1,  as 
obtained  using  exact  and  approximate  PDF’s,  averaged  over  the  intersymbol  inter¬ 
ference  (ISI)  patterns.  The  click  average  used  by  the  approximations  was  estimated  by 

q,  =4^e  for  pattern  t.  (4-27) 

Lt r 

For  high  Ej/Xq,  when  e  <C  1  and  also  a  <  1,  the  approximations  all  converge  to  Gaussian 
approximations  and  give  the  same  result,  as  shown,  which  underestimates  the  error  as  if 
0.15dB  were  added  to  Efc/JT0.  For  Eb/K0<  5.5dB,  all  of  the  approximations  over¬ 
estimate  the  error  probability;  it  is  somewhat  surprising  that  the  use  of  a  purely 
Gaussian  approximation,  taking  care  to  match  only  the  first  two  cumulants  of  the  exact 
differentia!  phase  distribution,  actually  results  in  the  best  fit  to  the  exact  Pf  curve. 


Gaussian 


FIGURE  4  10  COMPARISON  OF  EXACT  AND  ANALYTICAL  MIXTURE  APPROXIMATION  PDF’S 


Gaussian  IF  filter  Pattern:  111 


FIGURE  4  11  COMPARISON  OF  EXACT  AND  SIMPLIFIED  MIXTURE  APPROXIMATION  PDF’S 
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4.3.2  Adaptive  Gain  Control  (AGC)  Scheme 

For  the  case  of  frequency-hopped  noncoherent  BFSK  in  partial-band  noise  jam¬ 
ming,  it  has  been  shown  that  an  “AGC”  or  adaptive  gain  control  ECCM  scheme  is 
successful  in  mitigating  the  effects  of  the  jamming  [21,  22].  Under  this  scheme,  the 
different  chips  constituting  a  transmitted  symbol  axe  weighted  in  inverse  proportion  to 
the  noise  power  on  the  hops  on  which  they  axe  individually  received.  Therefore, 
analysis  and  computations  were  performed  to  evaluate  the  jammed  performance  of 
FH/CPFSK  when,  as  discussed  in  Section  2.1.4,  the  differential  phase  samples  {zt,} 
from  different  hops  are  weighted  by  the  factors  {w,}: 

L 

Z  k  =  Zkr  (4-2Sa) 

9 

Instead  of  the  weights  being  the  inverses  of  noise  power  (<r2)  as  previously  discussed,  in 
this  study  they  are  taken  to  be  the  inverses  of  noise  standard  deviation  ( o )  according  to 
the  relation 

w,  =  1 joq  (4-28b) 

=  1/(7  V.  hop  not  jammed  (4-2Sc) 

=  1/c rT,  hop  jammed.  (4-2Sd) 

The  use  of  o  instead  of  <r2  is  reasonable  in  view  of  the  fact  that  the  variance  of  the 
phase  noise  is  proportional  to  that  of  the  IF  noise:  that  is,  when  corrupted  by  noise,  the 
differential  phase  samples  are  random  variables  with  scale  factors  {<r,}.  Normalization 
by  <7,  then  has  the  effect  of  equalizing  the  noise  power  received  on  each  hop.  thereby 
reducing  the  tendency  of  the  stronger  noise  on  the  jammed  hops  from  dominating  in  the 
sum  of  weighted  differential  phase  samples. 

4.3.2. 1  Analysis 

Without  loss  of  generality,  we  can  formulate  the  error  probability  for  the  L 
hops/symbol  FH/CPFSK  system  when  an  AGC  scheme  is  used  by  assigning  the  weights 
w0=  1  to  unjammed  hops  and  uj,  to  jammed  hops,  where 

w\  =  ff  yv  /  °  x  =  J  p  t  /  P  s  =  I/''!"*-  (4-29) 

In  Section  2.2.1.  an  analysis  of  possible  jamming  events  made  use  of  a  conditional 
characteristic  function  (CHF)  for  the  sum  of  the  samples,  conditioned  upon  the  number 
of  hops  which  have  been  jammed.  An  unconditional  CHF  for  a  weighted  differential 
phase  sample  (>kq=vcqZkq  can  be  written 

C^u:  ■>)  =  (1  ~y)C((o:  <r2v)  +  7  C((  v;  k<t2v). 


r 


G4 


(4-30) 
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The  advantage  of  the  unconditional  form  is  that  the  CHF  for  the  weighted  sum  is 
simply  the  Xth  power  of  the  single-sample  CHF  (4-30),  giving  for  the  zth  ISI  data 
pattern  the  error  probability  expression  [a  special  case  of  (2-56)] 


’.<(£;  r)  =  I  CM-;  7)l‘sm(Laig!C>i  7)] ), 


(4-31) 


which  must  be  averaged  over  the  variation  in  parameters  due  to  the  ISI  patterns,  as 
discussed  in  Appendix  A. 

Now,  making  use  of  the  FM  click  noise  analysis  of  Section  2.3,  we  use  (2-59)  to 

write  (suppressing  the  ISI-pat tern-dependent  notation) 

rt  ■  \  n  \  -Qo(l-cos2wi/)-ja0sin2jri/  , 

Cc(t/;  *>)  =  (l-t)e  C^{u; 

-Q1(l-cos2?ru)1i/)-ja1sin2iru1i/  -  ,  . 

+  *>e  C^{wxv;  pT)  (4-32) 

in  which  the  pattern-dependent  CHF  for  the  modulo-2?r  differential  phase  without  clicks 
is  denoted  by  p)  when  the  SNR  is  p,  and  the  click  averages  q,  and  a0  are  given  by 

(2-28)  under  jammed  and  unjammed  conditions,  respectively.  For  convenience,  we 
separate  the  signal  and  noise  components  of  the  differential  phase: 

C  =  Ac  +  A  r).  (4-33a) 

Then.  Cv(i/;  p)  can  be  written 

p)  =  v\  p),  (4-33b) 

and  we  use  the  mixture  approximation  (4-25),  calibrated  by  taking  <tq  to  be  the  inverse 
of  the  variance  o\-=-o\n  and  explicitly  accounting  for  the  possibility  of  a  nonzero  phase 
noise  mean  value  pJr)  by  writing 


p)  =  e  r,{e{p)-Cv{i')  +  [1  -f(p)]-C,G(«r0i/)}. 


(4-33c) 


The  required  mean  and  variance  may  be  calculated  using  the  exact  modulo-2^  differ¬ 
ential  phase  probability  density  function  p^x)  given  as  (2-29): 


and 


VAr,  =  PAr,(p)  =  [  d XX  p^(xj  A<£=0) 

J  —  ir 

^  Art  =  ffAr)(p)  —  f  dx(x  ~PAn)  P^(I>  A($  =  0). 

J  —w 


(4-34a) 

(4-35a) 
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The  analytical  forms  of  the  magnitude  and  phase  of  the  CHF  needed  to  perform 
the  calculation  indicated  in  (4-31)  are  somewhat  tedious  but  are  straightforward  to 
derive.  Again  suppressing  the  notation  necessary  to  remind  us  of  the  ISI  pattern  depen¬ 
dence,  the  real  and  imaginary  parts  are  given  by 

Re{Q(i/;  7)}  =  (l-7)e  a°^  C°s2,rt/^cos [v(A<^  +  /io)-Qosin2?rj/] 

*  taCj/W  +  U-'o)  CG{ff00v)} 

+  7e"’a^1_c°s“,rUll/^cos[u;1i/(A^-i-/ii)-a1sin27rt(JiJ/] 

*  [<iC(/(wii/)  +  (l-cx)CG(ie1<r01i/)]  (4-36a) 

and 

Im{Q(i/;  7)}  =  (l-7)e~Q°*1_COS“,ri/^sin[i/(A^-f/io)-oosin25rv] 

x  [f0CV(v)  +  (l-co)  CG(croot')} 


+  *>e 


-a1(l-cos2xUjt/) 


sin  [iyi«/(A^  +  A,i)  — oisin2iru;1^] 


x  [qC7t.(tw1v)  +  (l-<1)  CG{v)xa0ii/)\,  (4-36b) 

in  which  we  have  used  the  subscripts  0  and  1  in  a  consistent  way  to  denote  values  of  the 
parameters  without  and  with  jamming,  respectively. 

The  total  error  probability  then  is.  after  averaging  over  the  three  ISI  data  patterns. 

P.(£;  i)  =  J 1)  +  P Ml;  l)  +  2P,3(£;  7)  ].  (4-37) 


4. 3. 2. 2  Numerical  results 


In  the  computations  of  FH/CPFSK  error  probability  for  L  hops  per  symbol  and 
partial-band  noise  jamming,  the  mixture  approximation  (4-25)  based  on  the  truncated 
Gaussian  distribution  was  utilized.  Inadvertently  instead  of  the  mixture  parameter  e(p) 
given  by  (4-25b),  we  in  effect  used 


I  _ 


ill 

l-( 


(4-3S) 


but  this  departure  from  the  analytical  procedures  described  previously  did  not  have  a 
great  influence  on  the  result,  since  computations  based  on  the  simplified  mixture 
parameters  given  by  (4-26)  did  not  produce  different  results  in  the  regions  of  the 
computed  curves  which  were  of  interest.3  Files  of  pre-computed  first  and  second 
cumulants  were  used  as  well  as  pre-computed  ISI  parameters  based  on  the  assumption 
that  /»  =  0.7,  W/fT  =  1.0,  and  that  a  Gaussian-shaped  IF  filter  is  used  in  the  receiver. 


3Further  comments  on  this  point  are  made  below. 
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In  Figure  4-13,  the  FH/CPFSK  uncoded  bit  error  for  L=  1  is  plotted  vs.  E b/Nj  for 
the  case  of  E4/X0  =  20dB  and  for  various  values  of  the  partial-band  jamming  fraction,  7. 
It  is  clear  from  the  figure  that  the  worst-case  value  of  7  (in  terms  of  maximizing  the 
error)  is  a  function  of  E b/K j.  The  upper  envelope  on  the  overlapping  curves  is  easily 
determined  to  be,  approximately, 


max  P« ( 1 ;  7)  =  P«(l;  7  =  1), 
_  0.238 

“  E»/V 


E b/Kj  <  ldB 
E JKj  >  ldB. 


(4-39a) 

(4-39b) 


This  dependence  is  the  well  known  “inverse  linear”  characteristic  of  the  error  perform¬ 
ance  of  communications  systems  in  worst-case  partial-band  noise  jamming.  It  also  can 
be  determined  from  the  figure  that  the  worst-case  value  of  7  is  approximately 

7wc  =  e^/^IST  ’  Ej/Jf j  >  ldB.  (4-39c) 


Although  not  shown  in  Figure  4-13,  there  is  a  different  asymptotic  value  for  Pt  for  each 
parametric  curve  to  the  left,  that  is,  as  E j/J'fj  decreases  (jammer  power  increases).  The 
value  of  the  asymptote  is  approximately  given  by  assuming  that  the  error  probability  is 
0.5  when  jammed,  and  negligible  (for  the  Ek/Xo=20dB  case  presented)  when  not 
jammed,  giving  the  value  7/2  for  the  asymptote. 

A  set  of  L=  1  curves  similar  to  those  shown  in  Figure  4-13  was  generated  for  the 
simplified  mixture  approximation  to  the  PDF  given  by  (4-26).  The  appearance  of  the 
curves  cannot  be  distinguished  from  those  in  Figure  4-13;  examination  of  the  computer 
printout  reveals  that  the  error  curves  computed  using  the  simplified  mixture  model  give 
a  value  for  the  error  that  is  generally  lower,  but  only  to  the  extent  that  the  second 
significant  figure  is  affected. 

Keeping  all  other  parameters  the  same,  in  Figure  4-14  we  show  the  effect  of 
increasing  the  value  of  the  diversity,  L,  from  1  to  2.  The  worst-case  Pe(2;  7)  can  be 
seen  to  have  a  generally  more  rapid  decrease  with  Eb/ifj  than  was  the  case  for  L=  1.  In 
fact,  for  values  of  E b/Kj  between  about  7  and  20  dB,  the  slope  of  the  upper  envelope  of 
the  parametric  curves  in  Figure  4-14  is  very  close  to  being  -2,  indicating  a  dependence 
upon  the  inverse  of  {Eb/X 7)2.  For  E b/}(j  greater  than  20dB,  for  the  range  of  Pe  values 
shown  the  slope  of  the  worst-case  error  is  approximately  -1,  indicating  a  dependence 
upon  (Et/J'fj)-1. 


For  the  curves  corresponding  to  7  values  of  0.02,  0.01,  and  0.005  in  Figure  4-14,  an 
unusual  dependence  on  E b/ X  j  can  be  observed.  The  curves  imply  that  there  is  a  peak 
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E»/*,  (dB) 

FIGURE  4-13  PERFORMANCE  OF  FH/CPFSK  IN  PARTIAL-BAND  JAMMING 
FOR  L=  1  HOP/BIT,  L/D  RECEIVER,  AND  AGC  WEIGHTING 
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E >/Xj  (dB) 

FIGURE  4-14  PERFORMANCE  OF  FH/CPFSK  IN  PARTIAL-BAND  JAMMING 
FOR  L= 2  HOPS/BIT,  L/D  RECEIVER,  AND  AGC  WEIGHTING 
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value  of  the  error  for  some  finite  value  of  Eb/Kj,  whereas  we  would  expect  the  error  to 
increase  uniformly  as  the  value  of  Eb/Nj  is  decreased  (jammer  power  increased).  This 
behavior  was  checked  out  thoroughly,  and  was  found  to  be  an  artifact  of  the  computa¬ 
tional  method,  which  uses  the  trapezoidal  rule  for  numerical  integration  to  calculate  the 
error  probability  according  to  (4-31).  The  spacing  between  the  samples  of  the  inte¬ 
grand  was  decreased  for  the  purposes  of  checking  the  result  for  fixed  7  in  Figure  (4-14), 
and  it  was  found  that  with  a  small  enough  spacing  the  error  probability  as  calculated  by 
this  integration  rule  does  indeed  behave  as  expected.  However,  the  value  of  the  result 
for  the  portion  of  the  curve  which  touches  the  upper  envelope  (worst-case  error)  was 
found  to  be  unaffected  to  any  significant  degree  by  the  decrease  in  spacing,  so  the  whole 
set  of  curves  was  not  recomputed  in  the  interest  of  saving  computational  time. 

Figures  4-15,  4-16,  and  4-17  show  the  effect  of  raising  the  value  of  L  to  3,  4,  and  6, 
respectively.  In  Figure  4-15,  for  L  =  3  the  slope  of  the  worst-case  error  probability  (the 
upper  envelope  of  the  parametric  curves  for  fixed  values  of  7),  is  seen  to  be  approx¬ 
imately  -3  for  E j/Xj  in  the  range  of  10  to  20  dB.  Similarly,  for  L= 4,  in  Figure  4-16 
the  slope  is  about  —4  for  the  same  general  range  of  Eb/Nj  values.  However,  in  Figure 
4-17,  the  slope  is  not  -6.  as  might  be  expected;  it  is  closer  to  being  -4.  This  behavior 
can  be  explained  by  the  fact  that,  even  though  Ek/X0  is  very  high  (20 dB),  as  the 
transmitted  symbols  (bits  if  no  coding  is  used)  axe  subdivided  increasingly  finer  by 
making  L  higher  the  system  becomes  subject  to  significant  noncoherent  combining 
losses. 

A  composite  of  the  worst-case  error  performances  of  FH/CPFSIv  for  different 
values  of  L  using  AGC  weighting  is  presented  in  Figure  4-18.  In  that  figure  we  observe 
that  for  low  values  of  E b/Nj  (high  values  of  jammer  power),  the  tendency  is  for  the 
error  to  increase  as  L  is  increased,  but  for  high  values  of  Eb/Kj  (low  values  of  jammer 
power),  the  tends  to  be  a  diversity  improvement  as  L  increases.  This  behavior  is 
reasonable  in  view  of  the  fact  that  for  strong  jamming,  the  worst-case  value  of  7,  the 
jamming  fraction,  tends  to  be  7  =  1,  corresponding  to  fullband  jamming,  and  we  know 
that  noncoherent  combining  losses  cause  the  error  to  increase  with  L  for  a  stationary 
Gaussian  noise  channel. 

For  reference,  in  Figure  4-18  there  also  is  plotted  a  curve  which  corresponds  to 
fullband  jamming  and  no  system  noise.  Against  this  reference,  we  can  observe  the 
effect  of  choosing  the  value  of  L  which  gives  the  lowest  worst-case  error  probability -the 
lower  envelope  of  the  curves  parametric  in  L.  From  this  point  of  view,  we  see  that  if  it 
were  possible  to  use  of  the  optimum  value  of  diversity  (which  depends  on  Eb/Xj)  the 
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E»/^ j  (dB) 


FIGURE  4-15  PERFORMANCE  OF  FH/CPFSK  IN  PARTIAL-BAND  JAMMING 


FOR  L=3  HOPS/BIT,  L/D  RECEIVER,  AND  AGC  WEIGHTING 
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FIGURE  4-16  PERFORMANCE  OF  FH/CPFSK  IN  PARTIAL-BAND  JAMMING 
FOR  1=4  HOPS/BIT.  L/D  RECEIVER,  AND  AGC  WEIGHTING 
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-5  5  15  25  35  45 

E*/ X j  (dB) 

FIGURE  4-17  PERFORMANCE  OF  FH/CPFSK  IN  PARTIAL-BAND  JAMMING 
FOR  L= 6  HOPS/BIT.  L/D  RECEIVER,  AND  AGC  WEIGHTING 
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FH/CPFSK  with  diversity  S~- 

Worst-case  partial-band  noise  jamming--— 
£  =  #  hope  per  bit,  E,/XD  =  20dB  §— 
AGC  chip  weighting,  L/D  detection  =- 
h=0.~,  W',rT  =  1.0,  Gaussian  IF  filter  — 


I  laiumaigBB 


s.w  t;';3 

B1U1 


il 


Reference:  fullbana  jamming 
with  no  system  noise 

(E,/Xp  =  xl 


(dB) 

FIGURE  4-18  DIVERSITY  PERFORMANCE  OF  FH/CPFSK  IN 

PARTIAL-BAND  NOISE  JAMMING  WITH  AGC  WEIGHTING 
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worst-case  error  performance  of  the  system  could  be  improved  greatly  from  the  situation 
holding  for  no  diversity.  Instead  of  requiring  about  44 dB  to  produce  a  1CT5  error 
probability. -a  loss  of  32 dB  in  performance  from  the  ideal  (fullband  jamming)  case, 
which  requires  about  12dB-with  optimum  diversity,  the  loss  due  to  the  worst-case 
jamming  can  be  held  to  less  than  10  dB  for  some  value  of  L. 

4.4  DIFFERENTIAL  DETECTION  RECEIVER  TECHNIQUES 

The  differential  detector  implementation  of  the  receiver  for  FH/CPFSK,  shown  in 
Figure  2-3  and  discussed  in  Section  2.1.3,  is  somewhat  easier  to  analyze  than  the 
limiter-discriminator  receiver.  Even  though  the  probability  density  function  (PDF)  for 
the  differential  detector's  output  is  not  available  in  closed  form,  its  characteristic 
function  (CHF )  is. 


4.4.1  Characteristic  Function  for  the  Differential  Detector  in  Jamming 

It  has  been  shown  [S]  that  the  output  of  the  differential  detector,  conditioned  on 
the  CPFSK  intersymbol  interference  parameters,  is  equivalent  to  the  difference  of  two 
scaled,  independent  noncentral  chi-squared  random  variables.  Let  z,,  denote  the  5th 
chip  on  a  particular  hop  of  an  L-hop/bit  FH/CPFSK  transmission ,  given  the  ith  inter- 
svmbol  interference  (ISI)  pattern,  where  i  may  be  1,  2.  or  3  (see  Appendix  A).  We  may 
write  the  CHF  for  the  jammed  receiver  as 


= 


+ 


1  c 

4  2 


(")  + 


\c^{uy 


(4-40a) 


where  the  pattern-dependent  CHF.  assuming  an  unspecified  ECCM  weighting  scheme 
using  weights  w0  and  w)  for  unjammed  and  jammed  hops,  respectively,  is  given  by 

=  ( 1  —  -> )  CAW'  A<?,.  j-U.o.  |\v,0,  trl) 

+  1  Cj(Wji/;  A0,,  jUjj,  jWjj,  <r 7- )  (4~40b) 

From  (2-35),  the  conditional  CHF  for  an  IF  filter  that  is  symmetric  about  the  center 
frequency  ( A  =  0 )  is 


n,  .  T.  2x  1  fji/<r2WsinA4+i<TV(U-rW 

c'("' 40' L'  "• ' 1  =  ui.vd-^rl - uVvii-^) 


cos  A*)' 


Without  loss  of  generality,  we  may  replace  z  with  the  normalized  variable. 


7.  = 


4  7. 


2  \  1 

wnr  si- r 


(4-41 ; 


(4-42) 
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which  yields  (after  leaving  off  the  prime) 
n  !  \  1  1  \  1  f 4^0 -Si'2,!?  ( 

c.,,(v)  =  <1—>)^?exP  (^7717  J 


+  1 


;exp 


f  4 j  v  a  (  ioj  j  w0 )  -  8  k* 2  /?  ( w:  /  w0 ) ' 

\  -  .  n  /  /  \  -> 


1+4ki/(w1/u;0)  [  1+4ki/  (Wj/wq) 


with 

and 


a  =  Wl0sin  A6i/ L-ll-r2 


(4-43a) 

(4-43b) 

(4-43c) 


,?  =  (Ui0  — r\Vt0cos  A<^,)/Ljl  — r2 , 

in  which  we  have  used  the  fact  that  the  jammed  and  unjammed  parameters  are  related 
by  the  parameter  *,  defined  in  (2-42).  in  the  following  way: 


A  ?  /  2 

K  —  (T  j/  <7  \ 

h,i  =  L  iC,/k  and  \\  ,i  —  W  ,o/K- 


and 


(4-43d) 

(4-43eJ 


4.4.2  AGC  Weighting  Scheme 

4.4.2. 1  Analysis 

For  the  AGC  weighting  scheme,  we  use  the  weights  w<j  —  \  and  K/;=er2./<T^-=  k_1  for 
the  differential  detector.  This  choice  is  logical  in  view  of  the  fact  that  the  differential 
detector  output  is  proportional  to  noise  power  (<r2).  Substitution  of  these  weights  into 
(4-43)  gives  the  CHF 


1 

1  +  4e2 


( 1  - 1 )  exp 


—  Sr2  J 
1  +  4r2 


exp 


4jra  i 
l+4i-2J 


+  7  exp 


—  Sr2/?/ 

1  +4r2  J 


exp 


4jro/K  )  \ 

l+4r2  J  J 


(4-44) 


The  imaginary  part  needed  for  the  error  probability  calculation  using  the  methodology 
described  in  Section  2.3  is 


I  m{[C,,/r)]i}  = 

1  C,„ 

(j')\L  sin  |Targ[C/fl( 

01} 

(4-45a) 

where 

|C;,,(l')!2  = 

{Re[ 

C„»j}2+  {Im[C. 

.»I)J 

(4-45b 1 

and 

[ 

[  Re  C„(r)j J 

arg '  C.-tlU')\  = 

tan"! 

( 4 -45c ) 

7G 


Combining  Techniques 


The  required  real  and  imaginary  parts  of  the  CHF  are  given  by 


Re[C;,i(t/)]  = 


_  1 


-8^^  ( 


4isa 


i+4„5  eitp\r^?/c“\rr4i?j 


+  7  exp 


{=££}“•{£&} 


and  "  ’  '  ^ 

lm[a„(,)]  =  (d-7) 


+  7  exp 


—  8  u"2  0  j  k\ 


sin 


\Al 'aj  f 


1  +4e2  J  [l+4i/\  , 

Using  numerical  integration  the  error  probability  for  the  t'th  pattern  becomes 
P„(i:  7)  = 

and  the  total  error  probability  is.  after  averaging  over  the  three  ISI  patterns, 
P <(L:  7)  =  UPel(l:  7 )  +  PAL:  7)  +  2P<3(I;  7)  ]• 


(4-46a) 


(4-46b) 


(4-47) 


(4-4S) 


4. 4. 2. 2  Numerical  results 

Calculations  of  the  differential  detector  error  probability  (4-4S)  were  made  using 
numerical  integration  of  (4-47).  In  Figure  4-19.  the  case  of  L=1  (no  diversity)  is 
presented.  In  this  figure.  Pf  is  plotted  as  a  function  of  Eb/Nj  for  fixed  values  of  7.  the 
partial-band  jamming  fraction,  when  E6/No=20  dB.  The  worst-case  jamming  result, 
the  upper  envelope  of  the  family  of  constant-7  curves,  is  seen  to  be  ‘‘inverse  linear"  in 
form,  that  is.  proportional  to  (E b/Kj)~]: 

maxPefl:  7)  =  0.275  x  (E b/Xj)~l  Eb/ br  j>2dB  (4-49a) 

=  P«(l:  1)  E*/Jr,<2dB.  (4-49b) 

This  result  is  in  agreement  with  the  comparable  result  for  equal  chip  weighting  pre¬ 
sented  in  [8],  which  was  computed  using  a  different  program.  For  L=  1.  of  course,  there 
is  no  effect  produced  by  using  any  kind  of  weighting,  since  the  bit  decision  is  based  on 
only  one  chip. 

Similar  results  for  L  =  1  are  shown  in  Figures  4-20  and  4-21  for  E6/N0  =  15dB  and 
13dB,  respectively.  As  the  value  of  Eb/H0  is  reduced  from  20dB,  the  asymptotic  error 
probability  for  high  E t/X;  rises  as  shown  in  these  two  figures.  The  fact  that  the 
parameters  of  the  worst-case  jamming  are  dependent  on  the  value  of  E6/X0  (that  is,  on 
the  amount  of  system  noise  present)  is  often  overlooked  in  jamming  analyses.  These 
figures  illustrate  that  the  results  can  be  quite  sensitive  to  the  assumed  value  of  Eb/S0. 
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Results  for  X= 2  for  E4/>fo=20dB  and  15dB  axe  presented  in  Figures  4-22  and 
4-23,  respectively.  In  Figure  4-22,  for  E b/tfj  greater  than  about  5dB,  the  upper 
envelope  of  the  constant-7  curves  has  slope  equal  to  —2  on  the  scale  shown,  indicating 
that  the  worst-case  probability  of  error  for  X=2  using  the  AGC  weighting  is 
proportioned  to  (Eb/Kj)~2.  This  result  in  itself  establishes  that  the  AGC  technique 
works  in  mitigating  the  effects  of  the  partial-band  jamming.  Note,  however,  from  what 
can  be  observed  in  Figure  4-23,  that  as  the  amount  of  thermal  or  system  noise  in¬ 
creases,  the  dependence  of  the  worst-case  error  on  Eb/K j  becomes  more  complex,  since 
there  is  a  “jamming-noise-dominant  region”  to  the  left  and  a  “thermal-noise-dominant 
region”  to  the  right. 

Similar  results  are  shown  for  X=3,  4,  5,  and  6  in  Figures  4-24  to  4-27.  The  fact 
that  the  increase  in  diversity  produces  a  deterioration  in  performance  for  low  values  of 
Eb/Kj  can  be  observed  by  noting  the  value  of  the  error  probability  at  Eb/K,=  — 5dB  in 
these  figures;  the  steady  increase  in  the  error  probability  is  understood  to  be  the  effect 
of  noncoherent  combining  losses,  since  for  strong  jamming  the  worst-case  value  of  is 
7  =  1  (fullband  jamming),  corresponding  to  a  stationary  Gaussian  noise  channel. 

The  fact  that  the  increase  in  diversity  tends  to  produce  an  improvement  in  perform¬ 
ance  for  high  values  of  Efc / -N" 0  and  E b/Kj  can  be  observed  in  Figures  4-19,  4-22,  and 
4-24  to  4-27  by  noting  the  worst-case  (maximum)  value  of  Ek/Xj  required  to  produce 
an  error  probability  of  10”5.  As  X  increases  to  X= 6,  this  value  is  seen  to  decrease.  How¬ 
ever,  another  trend  observable  from  these  figures  is  that  the  range  of  E k/tfj  values  for 
which  7  =  1  yields  the  worst-case  error  probability  tends  to  increase  as  the  order  of 
diversity  X  is  increased.  The  worst-case  value  of  7  (among  those  plotted)  at  Pe  =  10~5 
for  Ej/>fo=20dB  is  seen  from  these  figures  to  be  7=0.0001,  0.01,  0.1,  0.2,  0.5,  and  1.0 
for  X  =  1,  2,  3,  4,  5,  and  6,  respectively.  Thus  for  X>6  the  worst-cr.ce  value  of  EJtfj 
needed  to  produce  a  10” 5  error  will  begin  to  increase  with  X. 

A  composite  of  the  worst-case  error  probability  performance  for  FH/CPFSK  in 
partial-band  noise  jamming  using  AGC  weighting  and  differential  detection  is  given  in 
Figure  4-28.  We  observe  first  that  for  low  values  of  E b/}fj  (strong  jamming)  the  error 
probability  increases  as  X  increases,  due  to  the  fact  that  for  low  Eb/Kj  the  worst-case 
jamming  is  fullband  jamming  (7  =  1),  and  noncoherent  combining  losses  are  proportional 
to  X  on  the  fullband  jamming  (Gaussian)  channel.  Next  we  observe  that  the  value  of 
E b/Kj  required  to  produce  Pt(X)  =  10-5  tends  to  decrease  as  X  is  increased  to  the  value 
of  6,  as  previously  noted,  indicating  a  “diversity  gain”  is  in  effect.  Examination  of  the 
“under-envelope"  of  the  fixed-X  curves  revealed  that  the  best  value  of  X-in  the  sense  of 
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FIGURE  4-22  PERFORMANCE  OF  FH/CPFSK  IN  PARTIAL-BAND  JAMMING 

FOR  1  =  2  .  DIFFERENTIAL  DETECTION,  AND  AGC  WEIGHTING 
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FIGURE  4-23  PERFORMANCE  OF  FH/CPFSK  IN  PARTIAL-BAND  JAMMING 
FOR  L  =  2  .  DIFFERENTIAL  DETECTION.  AND  E»/Jf0=15dB 
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FIGURE  4-24  PERFORMANCE  OF  FH/CPFSK  IN  PARTIAL-BAND  JAMMING 
FOR  L= 3.  DIFFERENTIAL  DETECTION,  AND  AGC  WEIGHTING 
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E t/Jfj  (dB) 

FIGURE  4-25  PERFORMANCE  OF  FH/CPFSK  IN  PARTIAL-BAND  JAMMING 
FOR  1=4.  DIFFERENTIAL  DETECTION.  AND  AGC  WEIGHTING 
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FIGURE  4-2G  PERFORMANCE  OF  FH/CPFSK  IN  PARTIAL-BAND  JAMMING 
FOR  L-b.  DIFFERENTIAL  DETECTION.  AND  AGC  WEIGHTING 
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E»/ }f j  (dB) 

FIGURE  4-28  DIVERSITY  PERFORMANCE  OF  FH/CPFSK  IN  PARTIAL-BAND 
JAMMING:  DIFFERENTIAL  DETECTION,  AGC  WEIGHTING 
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giving  the  least  value  of  worst-case  Pe-is  a  function  of  E6/ JM 'j  (or,  looking  at  it  another 
way,  a  function  of  P<).  If  the  best  value  of  L  is  always  used,  the  performance  of  the 
system  against  worst-case  partial-band  noise  jamming  is  improved  from  the  very  poor 
performance  of  no  diversity  {L= 1)  to  within  3dB  of  the  system’s  performance  in  full- 
band  (non-worst-case)  jamming.  Thus  the  use  of  optimum  diversity  is  capable  of 
limiting  to  3dB  the  effects  of  the  jammer’s  efforts  to  inflict  the  maximum  possible 
disruption  of  the  system. 

It  is  interesting  to  compare  the  set  of  curves  in  Figure  4-28  with  those  for  the 
limiter-discriminator  receiver  shown  previously  in  Figure  4-18.  In  doing  so  we  first 
observe  from  the  “reference”  curves  in  both  figures  that  in  fullband  jamming  the 
limiter-discriminator  FH/CPFSK  receiver  performs  better  than  the  one  using  differ¬ 
ential  detection,  a  well  known  result.  For  example,  using  the  values  of  h= 0.7  and 
W jp T  —  1 . 0  the  limiter-discriminator  requires  about  11.6dB  to  produce  Pe  in  the 
Gaussian  channel,  compared  to  14.7dB  for  the  differential  detector.  Yet,  in  worst-case 
partial-band  noise  jamming  with  both  receivers  using  AGC  weighting  and  optimum 
diversity,  we  observe  that  the  differential  detector  outperforms  the  limiter-discriminator 
by  requiring  a  smaller  value  of  E b/tfj  to  achieve  the  same  error  probability.  For 
example,  the  optimum  diversity  for  both  receivers  to  achieve  Pe  =  l(T5  is  L= 6;  the 
differential  detector  requires  Ej/J^ySslTdB  to  accomplish  this  level  of  performance, 
while  the  limiter-discriminator  requires  Ej/Jf  Jw23dB. 

The  clue  to  understanding  the  better  performance  of  the  differential  detector  when 
diversity  is  employed  for  FH/CPFSK  is  to  compare  the  variation  in  E b/Nj  required  to 
achieve  P*=0.10  in  both  Figures  4-18  and  4-28  as  the  value  of  L  is  increased  from  L= 1. 
The  curves  in  Figure  4-18  are  spaced  further  apart  than  those  in  Figure  4-28;  this 
demonstrates  that  the  limiter-discriminator  receiver,  although  it  enjoys  a  “head  start” 
or  advantage  over  the  differential  detector  when  L=  1,  is  subject  to  higher  noncoherent 
combining  losses  (associated  with  FM  click  noise  as  discussed  in  Section  4.1)  than  is  the 
differential  detector.  This  point  is  also  brought  out  by  the  comparison  shown  in  Figure 
4-29,  in  which  we  have  given  both  types  of  receiver  the  same  performance  (Pt  =  10-5) 
for  L— 1  and  no  jamming.  The  asymptotic  performances  of  the  receivers  for  high  E  JN  j 
values  is  the  performance  in  noise-only  for  the  stated  values  of  Et/K0.  We  observe  in 
Figure  4-29  that  the  larger  noncoherent  combining  losses  prevent  the  use  of  L= 2  from 
improving  the  system  performance  when  a  limiter-discriminator  is  employed,  for  this 
level  of  system  performance  without  jamming;  however,  the  differential  detector,  which 
has  smaller  noncoherent  combining  losses,  yields  improvement  for  L  as  high  as  1=3. 
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brought  out  by  the  comparison  shown  in  Figure  4-29,  in  which  we  have  given  both 
types  of  receiver  the  same  performance  (Pt  =  10-5)  for  L=  1  and  no  jamming.  The 
asymptotic  performances  of  the  receivers  for  high  E b/tfj  values  is  the  performance  in 
noise-only  for  the  stated  values  of  Eb/}f0.  We  observe  in  Figure  4-29  that  the  larger 
noncoherent  combining  losses  prevent  the  use  of  L= 2  from  improving  the  system 
performance  when  a  limiter-discriminator  is  employed,  for  this  level  of  system 
performance  without  jamming;  however,  the  differential  detector,  which  has  smaller 
noncoherent  combining  losses,  yields  improvement  for  L  as  high  as  L= 3. 

The  somewhat  poor  performance  of  the  limiter-discriminator  for  L  >  1  can  perhaps 
be  improved  by  adjusting  the  receiver  parameters  h  and  W/fT.  In  [5],  for  example,  it  is 
said  that  the  performance  of  FH/CPFSK  in  worst-case  partial-band  noise  jamming  is 
optimized  for  /i=0.6  and  W/fT=0.75.  Variation  in  these  parameters  to  find  the  best 
combination  under  the  jamming  conditions  is  outside  the  scope  of  the  work  summarized 
in  this  report,  however. 

Another  interesting  way  to  look  at  what  we  observe  in  this  comparison  is  that  the 
differential  detector,  whose  output  is  proportional  to  sinA4>,  performs  further  processing 
on  the  limiter-discriminator  output,  A4>:  the  sine  function  both  limits  the  value  of  the 
receiver  output,  since  |sinx|<l,  and  nullifies  the  effects  of  any  2v  phase  shifts  (clicks) 
superimposed  on  A4>  due  to  the  phase  ambiguity  property  of  the  sine  function: 
sin(x±2*-)  =  sinx. 


I 

I 


4. 5. 1.2  Comparison  of  limiter-discriminator  results 

We  now  compare  the  uncoded  receiver  performances  using  diversity  combining  and 
limiter-discriminator  detection  shown  in  Figure  4-18,  for  AGC  soft-decision  combining, 
and  in  Figure  4-6.  for  hard-decision  combining.  Both  figures  assume  Eb/N0  =  20dB.  It 
is  apparent  at  once  from  comparing  these  two  sets  of  curves  that,  for  limiter- 
discriminator  detection  of  FH/CPFSK,  the  hard-decision  receiver  performs  better  than 
the  soft-decision  receiver.  For  example,  using  the  best  value  of  L ,  a  10-5  error 
probability  is  achieved  for  Ej/X^ndB  for  hard-decision  combining  and  for 
Et/ifj  =  23dB  for  soft-decision  combining.  This  is  a  surprising  result,  since  soft-decision 
procedures  generally  outperform  hard-decision  ones. 

The  reason  for  this  surprising  result  is  that  when  diversity  is  used,  the  bit  or 
channel  symbol  energy  is  split  into  L  chips,  with  the  consequence  that  the  value  of  the 
chip  energv-to-noise-density  ratio  is  decreased  from  E J}f0  by  the  factor  \/L.  This 
decrease  in  signal  energy  per  observation  interval  gives  rise  to  a  disproportionately 
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FH/CPFSK  with  diversity 
Worst-case  partial-band  noise  jamming 
£  =  #  hops  per  bit,  AGC  chip  weighting 
-  Differential  detection  (Et/>f0  =  14.7dB)  - 

---  -  Limiter-discrimator  (Ek/X0  =  11.6dB)  . 

(These  Ej/X0  values  give  10-5  for  £  =  1  and 

,  no  jamming)  i 

A=0.7,  W'/fT=  1.0,  Gaussian  IF  filter 


E ,/*,  (dB) 

FIGURE  4-29  COMPARISON  OF  THE  DIVERSITY  PERFORMANCES  OF 
FH/CPFSK  IN  PARTIAL-BAND  JAMMING 
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higher  risk  of  FM  noise  clicks,  especially  on  jammed  hops.  Now,  if  soft-decision 
combining  is  used,  the  sum  of  the  differential  phase  samples  can  be  “contaminated"  by 
the  occurrence  of  a  click  on  any  or  all  of  the  hops  from  which  the  samples  axe  taken. 
For  example,  if  L= 2,  supposing  that  one  hop  is  jammed  and  one  is  not,  the  value  of  the 
sum  of  the  samples  is 

w1-A$1  +  w0-A$2  =  A<)>(w1+  w0)  —  2ir(w1Nj  +  w0N2)  +  wlAr]l  +  w0Atj2 

=  A^(1  +  ^)_2,r(^  +  Ar2)  +  7^  +  A,?2  (4~50) 

in  which  A<j>  is  the  nominal  (noiseless)  differential  phase  (assumed  positive),  the  Ar/’s 
are  the  differential  phase  noise  samples,  *  is  the  ratio  of  jammed  noise  power  to 
unjammed  noise  power,  and  the  Ns  are  click  numbers  jamming— whose  averages  are  qj 
and  qq  with  and  without  jamming,  respectively.  Typically,  A<i>  is  between  0.5tt  and  rr. 
and  a0  is  very  small.  Thus  while  the  AGC  weighting  takes  care  of  the  phase  noise,  in 
effect  equalizing  the  noise  for  jammed  and  unjammed  hops,  it  does  not  remove  the 
clicks,  if  and  when  they  occur,  although  their  contributions  on  jammed  hops  are 
attenuated. 

Now,  when  hard  decisions  on  each  hop  are  used,  any  false  decision  that  is  due  to 
jamming  or  to  the  occurrence  of  clicks  (which  are  disproportionately  more  likely  on 
jammed  hops)  is  automatically  “contained”  to  one  “vote”  in  the  overall  T-hop  decision. 

4. 5. 1.3  General  comparison  without  coding 

Comparing  now  the  limiter-discriminator  FH/CPFSK  receiver  with  either  hard  or 
soft  decisions  and  the  differential  detector  with  soft  decisions,  it  would  appear  from  the 
computed  results  that  the  choice  of  receiver  among  these  three  to  mitigate  the  effects  of 
■worst-case  partial-band  noise  jamming  will  depend  on  the  value  of  EJXo-  If  there  is 
relatively  little  thermal  noise  (high  Ej/X0),  the  ECCM  receiver  of  choice  would  be  the 
combination  of  hard  decisions  and  the  limiter-discriminator.  If  the  thermal  noise  is 
definitely  not  negligible  (Eb/K0<  20dB),  the  differential  detector  with  AGC-weighted 
soft  decisions  should  be  used. 

4.5.2  Coded  Performance  Results  and  Comparisons 

A  number  of  computations  were  made  of  the  performance  of  a  coded  FH/CPFSK 
system,  assuming  that  an  “outer  code”  is  used  in  conjunction  with  an  “inner  code” 
consisting  of  one  ot  the  combinations  of  receiver  design  and  L- hop/code  symbol  ECCM 
combining  techniques.  As  discussed  in  Section  2,  the  error-correcting  outer  code 
selected  for  the  study  is  a  rate-i,  constraint  length  7,  convolutional  code,  with  either 
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hard  or  soft  decoding.  The  comparisons  to  be  made  are  among  the  following  cases: 

A.  Limiter-discriminator  detection  with 

1.  Hard-decision  diversity  combining,  hard  outer  decoding 

2.  Soft-decision  AGO  diversity  combining,  hard  outer  decoding 

3.  Hard-decision  diversity  combining,  soft  outer  decoding. 

B.  Differential  detection  with 

1.  Hard-decision  diversity  combining,  hard  outer  decoding 

2.  Soft-decision  AGC  diversity  combining,  hard  outer  decoding 

3.  Hard-decision  diversity  combining,  soft  outer  decoding. 

The  methodology  for  calculating  upper  bounds  on  the  bit  error  probability  when 
coding  is  used  -was  summarized  in  Section  2.2.2  (pp.  27,  28).  Since  the  rate  of  the 
selected  code  is  r=i,  note  that  numerical  values  for  these  bounds  as  a  function  of  Eb/Kj 
can  be  thought  of  as  being  developed  in  two  steps:  (1)  Find  the  value  of  the  uncoded 
error  probability  for  3dB  less  than  E b/Xj  AND  3dB  less  for  Eb/X0:  and  (2)  map  the 
value  of  the  uncoded  error  probability  into  the  value  for  the  coded  error  probability 
bound. 

For  example,  equations  (2-48)  through  (2-51)  imply  that  10-7,  10~6,  10-5,  10~\ 
10-3,  and  10-2  hard-decision  decoding  error  probability  bounds  result  from  code  symbol 
error  probabilities  of  0.007177,  0.01116,  0.01711,  0.02572,  0.03773,  and  0.05409, 
respectively.  Thus  if  we  are  interested  in  estimating  the  region  of  Eb/Nj  for  which  the 
decoded  bit  error  probability  is  between  10” '  and  10-2,  then  we  may  look  at  the  curves 
for  the  uncoded  error  probability  and  find  the  range  of  Eb/Xj  for  which  that  quantity  is 
between  0.0072  and  0.054;  to  these  E b/Xj  values  we  then  add  3dB.  Applying  this 
technique  to  Figure  4-6  (page  53),  we  find  that  for  Ej,/jfo=20dB  the  uncoded  error 
probability  is  in  the  range  of  0.007  to  0.05  for,  approximately,  5.5dB  <  Eb/Xj  <  14.5. 
implying  that  for  hard-decision  decoding  the  coded  error  probability  will  be  between 
10-7  and  10-2  for  Eb/Xj  between  8.5 dB  and  17.5 dB.  Some  approximation  is  involved 
because  the  curves  in  Figure  4-6  are  valid  for  E4/.N"o  =  20dB,  not  17dB. 

While  this  technique  is  useful  for  estimating  coded  results,  its  primary  usefulness  is 
for  checking  the  results  for  reasonableness.  In  all  of  the  coded  results  presented  below, 
the  error  probabilities  were  recalculated  as  functions  of  Eb/X  j  for  fixed  values  of  Eb/X0 
and  of  7.  the  partial-band  jamming  fraction,  and  from  these  parametric  curves  the 
worst-case  performance  was  determined. 
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4.5.2. 1  Results  for  hard-decision  diversity  combining  and  hard  decoding 

In  Figure  4-30,  the  worst-case  partial-band  noise  jamming  performance  of 
FH/CPFSK  is  presented  as  a  function  of  E b/K  j  for  Eb/Xo=20dB,  for  different  values  of 
L,  the  number  of  hops  per  code  symbol,  and  for  both  limiter-discriminator  and  differ¬ 
ential  detection,  assuming  hard-decision  diversity  combining  and  hard  decoding.  Inter¬ 
symbol  interference  has  been  accounted  for  in  the  calculations,  using  the  parameter 
values  h= 0.7  and  W/jrT=1.0,  and  assuming  a  Gaussian  IF  filter. 

We  first  observe,  as  a  check  on  these  results,  that  the  range  of  Eb/)f  j  for  which  the 
coded  bit  error  probability  falls  between  10-7  and  10~2  is  in  agreement  with  that 
predicted  from  the  relevant  uncoded  results  shown  previously  in  Figure  4-6,  for  the 
limiter-discriminator  receiver.  (Uncoded  hard-decision  combining  performance  results 
were  not  computed  for  the  differential  detector.) 

The  coded  performance  of  the  receiver  using  limiter-discriminator  detection  is  seen 
to  be  uniformly  better  than  that  for  the  receiver  using  differential  detection.  While  for 
the  unjammed  Gaussian  channel  the  limiter-discriminator  is  about  3.1  dB  better  com¬ 
pared  to  the  differential  detector  (using  the  value  of  Eb/}C0  required  to  obtain  a  10-5 
error  probability  as  a  basis  for  comparison),  we  see  from  Figure  4-30  that  this  advan¬ 
tage  is  reduced  to  about  1.3 dB  in  worst-case  partial-band  noise  jamming  when  coding 
arid  diversity  are  employed.4  This  same  margin  is  observed  for  no  coding  and  no 
diversity  in  worst-case  jamming,  when  the  BER  curve-fit  approximation  (4-14)  is  used 
to  calculate  the  jammed  error  probability  for  the  limiter-discriminator  receiver. 

For  both  receiver  types,  the  optimum  value  of  L  (“optimum  diversity”)  is  seen  to 
depend  on  the  value  of  Eb/Xj.  For  achieving  Pt  =  10-5,  that  value  is  L= 3,  in  contrast 
to  the  much  higher  value  (<  11)  without  coding.  The  reason  why  a  lower  value  of  diver¬ 
sity  is  optimum  with  coding  than  without  coding,  for  achieving  the  same  error  proba¬ 
bility,  is  that  with  coding  the  symbol  error  probability  at  the  output  of  the  diversity 
combining  portion  of  the  receiver  is  much  higher  (between  .007  and  .05,  as  we  have 
seen)  and  is  in  the  region  of  the  curves  parametric  in  L  for  which  lower  values  of  L  tend 
to  be  optimum;  quite  often  in  this  region  the  worst-case  value  of  y,  the  partial-band 
jamming  fraction  is  close  to  unity,  therefore  indicating  a  nearly  Gaussian  channel  under 
the  noise  jamming,  so  that  an  increase  in  L  from  L=  1  tends  incurs  a  noncoherent  com¬ 
bining  loss  in  performance  for  a  lower  value  of  L  than  without  coding. 


4 A  summary  comparison  of  E b/N j  requirements  to  achieve  Pt  =  ]0-5  with  and  without  diversity  and 
coding  will  be  given  below. 
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FIGURE  4  30  CODED  PERFORMANCE  WITH  HARD-DECISION  COMBINING  AND  HARD  DECODING 
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4. 5. 2. 2  Results  for  soft-decision  diversity  combining  and  hard  decoding 

The  worst-case  partial-band  noise  jamming  coded  performance  of  limiter-discrim¬ 
inator  and  differential  detection  FH/CPFSK  receivers  is  presented  in  Figure  4-31  for 
the  case  of  A GC  soft-decision  diversity  combining  and  hard  decoding  when  Eb/K0  = 
20 dB.  It  might  have  been  anticipated  from  the  relative  performances  of  these  two 
receivers  with  AGC  combining  but  without  coding  that  the  differential  detector  would 
display  a  better  performance  with  coding.  However,  the  results  in  Figure  4-31  indicate 
that  the  limiter-discriminator  consistently  achieves  a  better  performance,  although  for 
L>  1  its  advantage  is  small,  on  the  order  of  ldB  or  less  for  a  given  value  of  L. 

For  L=  1,  the  advantage  of  the  limiter-discriminator  over  the  differential  detector 
at  the  10~5  bit  error  probability  level  shown  in  Figure  4-31  is  about  3.7dB.  It  should 
be  noted  that  the  same  quantity  is  about  1.4  dB  in  Figure  4-30,  or  2.3  db  less,  even 
though  in  principle  the  combining  and  decoding  portions  of  the  respective  receivers  are 
identical  when  L=  1.  This  difference  is  apparently  due  in  part  to  the  use  of  two  differ¬ 
ent  approaches  to  calculating  the  probability  of  error  at  the  output  of  the  combiner 
(input  to  the  decoder).  For  hard-decision  combining,  as  noted  in  connection  with 
Figure  4-30,  an  approximation  (curve  fit)  was  used  since  it  was  sufficient  for  those 
calculations  to  utilize  a  BER  curve  for  L=  1.  For  AGC  soft-decision  combining,  the 
probability  of  coded  symbol  error  was  computed  using  a  mixture  approximation  for  the 
characteristic  function  of  the  limiter-discriminator  output  samples.  Although  as  shown 
in  Figure  4-3  (page  50)  the  curve-fit  approximation  is  very  good  for  the  Gaussian 
channel  — the  two  computational  methods  agree— evidently  for  worst-case  jamming  they 
disagree  slightly,  with  the  difference  being  magnified  when  coding  is  applied. 

For  the  optimum  values  of  AGC  diversity  combining,  the  limiter-discriminator 
receiver  achieves  a  10~5  bit  error  probability  for  about  a  1.1  dB  smaller  value  of  E fc/ JsT ^ 
than  does  the  differential  detection  receiver. 

4. 5.2. 3  Results  for  hard-decision  diversity  combining  and  soft  decoding 

Figure  4-32  gives  the  coded  bit  error  probability  for  FH/CPFSK  in  worst-case 
partial-band  noise  jamming  for  both  limiter-discriminator  and  differential  detection 
types  of  receiver,  assuming  hard-decision  diversity  combining,  soft-decision  decoding, 
and  Et/Jfo=20dB.  As  in  Figures  4-30  and  4-32,  the  limiter-discriminator  is  seen  to 
achieve  a  better  performance  consistently  for  the  range  of  Pe  values  shown.  At  the  10"5 
level  of  performance,  when  optimum  diversity  is  used,  the  limiter-discriminator’s  advan¬ 
tage  is  about  1.5dB. 
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FIGURE  4  31  CODED  PERFORMANCE  WITH  AGC  SOFT-DECISION  COMBINING  AND  HARD  DECODING 


FIGURE  4  32  CODED  PERFORMANCE  WITH  HARD  DECISION  COMBINING  AND  SOFT  DECODING 
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In  comparing  coded  results  for  L—\  shown  in  Figures  4-30,  for  hard  decoding,  and 
4-32.  for  soft  decoding,  we  note  that  there  are  differences  in  the  values  of  EJKj 
required  to  achieve  a  particular  value  of  bit  error  probability,  with  these  values  being 
higher  for  soft -decoding.  This  draws  attention  to  the  fact  that  both  performances  are 
actually  upper  bounds  to  the  probability  of  bit  error,  not  the  actual  bit  error 
probability,  and  the  two  bounding  techniques  give  rise  to  different  numbers  for  the 
same  conditions.  When,  for  example,  L= 3  is  considered,  it  is  readily  apparent  that  the 
soft  decoding  produces  about  a  1.4  dB  gain  over  the  hard  decoding  by  making  use  of  the 
additional  information  made  available  to  it.  This  number  is  very  close  to  the  1.5  dB 
gain  predicted  in  [1]. 

4. 5. 2. 4  Summary  comparison  of  coded  performance  results 

Supposing  that  optimum  diversity  is  always  used,  a  kind  of  lower  bound  (with 
respect  to  the  order  of  diversity)  to  the  error  probability  upper  bounds  shown  in  Figures 
4-30  to  4-32  can  be  constructed  as  shown  in  Figure  4-33.  It  is  apparent  from  this 
figure  that  the  variation  in  the  value  of  E JKj  required  to  obtain  a  bit  error  of  10-5  is 
relatively  slight,  being  within  a  range  of  about  3dB  for  different  combinations  of 
receiver  type,  diversity  combining  technique,  and  decoding  technique. 

A  consistent  pattern  (in  addition  to  the  advantage  of  using  limiter-discriminator 
detection  instead  of  differential  detection)  evident  in  Figure  4-33  is  the  following 
ranking  of  the  combinations  of  diversity  combining  and  decoding  techniques  that  were 
studied:  generally  for  Pt  values  of  10'5  or  less,  best  performance  is  obtained  using  hard- 
decision  combining  and  soft  decoding.  AGC  soft-decision  combining  plus  hard  decoding 
is  second  best,  and  the  combination  of  hard-decision  combining  and  hard  decoding  is 
third  best.  Some  form  of  soft  processing  therefore  is  advantageous. 

It  is  interesting  to  note  that  for  10-,<<  Pe<  10~2,  the  ranking  of  soft  decoding  and 
soft  combining  is  reversed.  The  reason  for  this  behavior  seems  to  be  the  fact  that  for 
these  values  of  the  error  probability,  the  optimum  diversity  is  L  =  1  for  soft  combining 
and  L=2  (or  3)  for  soft  decoding.  However,  as  noted  before  the  relative  behavior  of  the 
two  decoding  error  probability  bounds  is  uncertain  for  L  —  1,  and  for  these  higher  values 
of  P«  the  tightness  of  the  bounds  can  be  expected  to  diminish. 

A  summary  of  E requirements  to  achieve  P,=  1CT5  when  Et/.N'0  =  20dB  is 
given  in  Table  4-1.  For  reference,  the  first  case  shown  in  the  table  is  the  worst-case 
FH/CPFSK  performance  for  no  coding  and  no  diversity,  noting  the  differences  obtained 
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TABLE  4-1 

REQUIREMENTS  FOR  10"5  Pe  IN  WORST-CASE  JAMMING  WHEN  Et/Xo=20dB 


FH/CPFSK:  h  =  0.7,  W/fT  =  1.0,  Gaussian  IF  filter 


Limiter-discriminator 

Differential  detector 

Case 

E*/-MdB) 

Gain 

E»/jr.,(dB) 

Gain 

Reference:  no  coding 
no  diversity 

43.01 

— 

44.43 

_ 

44. 12 

... 

No  coding,  optimum  diversity 

hard  combining 

<16.7'  (1=11) 

>26.3 

not  computed 

AGC  combining 

<22.S2  (6) 

>21.3 

17.0  (6) 

27.4 

Coding  only  (Z,  =  l ) 

hard  decoding 

13.S1 

29.2 

15.2 

29.2 

11.52 

32.6 

soft  decoding 

15.01 

2S.0 

16.4 

2S.0 

Coding,  optimum  diversity 

hard  comb.,  hard  decoding 

12.31  (3) 

30.7 

13.7  (3) 

30.7 

AGC  comb.,  hard  decoding 

11.02  (2) 

33.1 

12.1  (2) 

32.3 

hard  comb.,  soft  decoding 

10.91  (3) 

32.1 

12.4  (3) 

32.0 

Notes: 

1.  All  hard-decision  combining  results  for  the  limiter-discriminator  were  calculated 
using  the  curve-fit  approximation  to  the  bit  error  probability. 

2.  All  AGC  soft-decision  combining  results  for  the  limiter-discriminator  were  calcu¬ 
lated  using  the  numerical  integration/characteristic  function  technique. 

3.  All  differential  detector  results  were  calculated  using  the  numerical  integration/ 
characteristic  function  technique. 
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for  the  limiter-discriminator  for  the  two  computational  methods  used.  The  second  case 
shown  in  the  table  is  that  of  optimum  diversity  without  coding,  for  which  it  is  evident 
that  diversity  alone  can  produce  a  gain  in  performance  a s  high  as  26  or  27dB,  in  terms 
of  reducing  the  required  E t/tfj. 

The  third  case  shown  in  Table  4-1  is  that  of  coding  alone,  without  the  use  of 
diversity  (i.e.,  L  =  1).  Here  we  observe  that  coding  alone  can  produce  a  gain  of  28dB  or 
higher.  It  is  somewhat  surprising  that  the  hard-decoding  gains  are  consistently  greater 
(by  about  ldB)  than  the  soft-decoding  gains  in  the  cases  studied.  But  on  the  other 
hand,  for  coding  alone  there  is  no  reason  to  expect  a  different  performance  with  soft 
decoding,  since  L—  1.  The  ldB  difference  therefore  is  due  to.  the  difference  in  the 
probability  of  error  upper  bounds  used. 

The  fourth  and  final  case  shown  in  Table  4-1  is  the  combination  of  optimum 
diversity,  combining  method,  and  decoding  technique.  Total  coding  plus  diversity  gains 
of  about  31-33  dB  are  seen  to  be  possible.  We  observe  also  that  the  addition  of 
diversity  to  coding  provides  an  incremental  gain  in  performance  of  about  2-4 dB, 
whereas  if  coding  is  added  to  diversity  an  improvement  of  about  5dB  is  experienced. 
Thus  it  would  appear  that  the  more  influential  factor  in  the  improvement  in  perform¬ 
ance  is  coding,  especially  when  it  is  considered  how  much  the  performance  is  improved 
using  coding  alone. 
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APPENDIX  A 

INTERSYMBOL  INTERFERENCE  ANALYSIS 


We  have  observed  that  the  IF  bandpass  filter  has  the  effect  of  distorting  the  data- 
modulated  phase  waveform  0m(t),  resulting  in  the  distorted  phase  waveform 


<j!>(t)  =  tan- 


f  ho(t)*sin0m(t)1 
|ho(t)*cos0m(t)  J’ 


(A- 1 ) 


as  well  as  in  an  amplitude  modulation,  a(t).  For  most  cases  of  practical  interest,  it  is 
sufficient  to  consider  intersvmbol  interference  (ISI)  effects  (i.e.,  the  overlapping  of  filter 
responses  from  different  transmitted  symbol  or  chip  intervals)  due  to  immediately 
adjacent  symbols  or  chips  [6].  Therefore,  in  what  follows  we  consider  chip  patterns  on 
the  FH/CPFSK  hops  which  are  the  periodic  extensions  of  the  patterns 

111.  000  (all  one’s  or  zeros); 

010,  101  (alternating  one’s  and  zeros);  (A-2) 

and  0110,  1001,  1100,  0Q11. 


The  “present  chip”  in  each  of  these  sequences  is  indicated  by  the  underlining.  The 
patterns  in  (A-2)  were  selected  because  they  generate  the  eight  possible  3-chip  patterns 
in  a  very  simple  manner  and  can  be  analyzed  easily. 


Using  the  steady-state  filtering  approach  of  Tjhung  and  Wittke  [19]  and  Pawula 
[11,  12],  we  recognize  that  the  periodic  extension  of  the  in-phase  and  quadrature  signal 
waveforms  for  the  ith  pattern  yields  the  Fourier  series 

PC 

sin^m(t;  pattern  i)  =  ^  a*,  cos{2xkfp\  +  ?ki)  (A-3a) 

and  i=1 

oo 

cos0m(t;  pattern  i)  =  y]  0ki  cos  (2irkfPit  +  nki),  (A-3b) 

k  =  1 


where  fPi  is  the  appropriate  fundamental  frequency.  Then  the  responses  of  the  IF  fil¬ 
ter’s  lowpass  equivalent  h0(t)  to  these  series  axe 

u(t;  pattern  i)  =  ho(t)*sin0m(t;  pattern  t) 

=  ^  \H0(kfpi)\akicos[2nkfPit  +  ipki  -  B(fc/P.)]  (A-4a) 

k  =  1 

and  v(t:  pattern  t)  =  ho(t)*cos0m(t;  pattern  i) 

=  £  \tto(kfr,)\0kiCos{2nkfpi  +  Tikl-B{kfpt)].  (A -4b) 

i  =  l 
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where  |H0(/)j  is  the  magnitude  of  the  lowpass  filter's  frequency  transfer  function  and 
B (/)  is  its  phase  delay  as  a  function  of  frequency.  The  analysis  is  considerably  simpli¬ 
fied  if  we  introduce  time  offsets  tu  and  U  such  that  sin  [0m(t— tu)]  and  cos[0  m  (t— tt.)]  , 
respectively,  are  even  functions.  Also,  we  assume  that  the  filter  has  a  linear  phase  delay 
characteristic,  that  is, 

B  (f)  =  2*ftd.  (A-5) 

Together,  these  simplifying  factors  allow  us  to  rewrite  (A-4a)  as 

A  _°c 

u(t;  pattern  i)  =  ^  |H0(fc/Pj)| aticos[2xfc/Pjft  +  tu+ td)]  (A-6a) 

i  =  l 

and  v(t;  pattern  i)  =  |H0(fc/Pi)|  Ski  cos  [2nkfPi(t  +  U+  td)] .  (A-6b) 

k=\ 

Patterns  111.  000.  For  the  patterns  111  and  000,  the  undistorted  signal  is  a  pure 


sinusoid,  with 

0m(t;  Ul)  =  -Mt;  000) 

=  irht/T  =  2irfdt.  (A-7) 

Thus  for  these  patterns  sin0m(t)  and  cos0m(t)  are  single-term  Fourier  series,  and  their 
distorted  versions  are 

u(t;  111)  =  -u(t:  000)  =  a0sin  [x/i(t  -  td)/T]  (A-Sa) 

and  v(t:  Ul)  =  u(t;  0Q0)  =  &0cos[nh(t  -  td)/T],  (A-Sb) 

using  the  parameter 

ao  —  |H0(/d)|-  ( A  -  9 ) 

These  quantities  in  turn  give  the  ISI  parameters 

A*(t;  111)  =  —  A0(t;  000)  =  *h  (A-lOa) 

a2(t;  111)  =  a2(t;  0Q0)  =  ag  (A-lOb) 

and  a2(t-T;  111)  =  a2(t-T;  0Q0)  =  a^.  (A-lOc) 


Note  that  for  these  patterns  the  parameter  values  are  independent  of  the  sampling  time. 

Patterns  010,  1Q1.  For  the  patterns  010  and  1Q1,  the  original  frequency  modulation 
is  a  ±fd  squarewave  with  period  2T,  so  that  0m(t)  is  an  even,  bipolar  triangular  wave 
with  amplitude  \rh  and  period  2T.  For  pattern  010,  the  triangular  wave's  first  positive 
peak  occurs  at  t  =  T,  using  the  convention  that  the  baseband  data  waveform  is  such 
that  d(t)  =  dk  for  (k—  l)T<t<kT.  We  note  that  pattern  1Q1  is  simply  the  negative  of 
pattern  010.  The  functions  sin0m(t;  010)  and  cos0m(t;  010)  therefore  are  even,  with 
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periods  2T  and  T,  respectively.  Expansion  of  these  functions  in  Fourier  series  yields 


sin0m(t;  010)  =  -sin0m(t;  1Q1) 


4/ico„  /i_t\  vs'c°s[(2A:  l)irt/T] 
_  _^cos  (Jirfc)  ^  -—p—j- 


(A-lla) 


cos0m(t;  010)  =  cos0m(t+T;  1Q1) 


2  (\  i  o cos  (2fcirt/T) 

=  asm  _• 


(A-llb) 


Assuming  the  filter  rejects  harmonics  with  />  1/T,  we  find  that  the  signal  quadrature 
components  for  these  patterns  are 


u(t;  010)  =  — u(t;  101) 

=  -#cos(i^)|H0(^)|  j-L,  cos^t-t^/T] 
=  -c1cos[5r(t-t<i)/T! 


(A-12a) 

(A-12b) 


and 

v(t;  010)  =  v(t;  1Q1) 

=  Asia  (J»k)[l-J4ij|H„(4,)|coS[2*(t-tJ)/Tl 

(A~12c) 

=  c2  -  c3cos[2ir(t-t<j)/T]. 

(A-12d) 

These  quantities  at  the  sampling  time  t  =  T  +  td  give  the  ISI  parameters 

Atf(trf+T;  010)  =  — AcKt,,  +  T;  1Q1) 

=  2 tan  ‘(cjzhs) 

(A-13a) 

a2(td  + T;  010)  =  a2(td  +  T;  1Q1)  =  c2  +  (c2  —  C3)2 

(A-13b) 

and 

a2(trfi  010)  =  a2(td;  1Q1)  =  c2  -f  (c2— c3)2. 

(A-13c) 

in  which  the  constants  are  defined  as 

C,  S  cos(H)|H„{^)| 

(A-14a) 

c2  =  -^sin(*7r/i)  =  sine  (\h) 

(A-14b) 

and 

Cj-i(^sin(i''“)|Ho(+)l' 

(A-14c) 

Distorted  and  undistorted  phase-related  waveforms  for  the  010  pattern  axe  illustra- 


105 


ISI  Analysis 


ted  in  Figure  A-l. 

Patterns  OHO,  1Q01.  For  patterns  01.10  and  1Q01,  the  original  frequency  modula¬ 
tion  is  a  ±/d  squarewave  with  period  4T,  so  that  0m(t)  for  these  patterns  is  an  even, 
bipolar  triangular  wave  with  amplitude  vh  and  period  4T.  For  pattern  0110,  the 
triangular  wave’s  first  positive  peak  occurs  at  t  =  2T.  Pattern  1Q01  is  the  negative  of 
pattern  0110.  The  functions  sin0m(t;  0110)  and  cos0m(t;  0110)  therefore  are  even,  with 
periods  4T  and  2T,  respectively.  Their  Fourier  series  expansions  are  most  easily  written 
down  as  those  for  pattern  010,  but  with  h  replaced  by  2 h  and  with  t  replaced  by  *t: 


sin0m(t;  0110)  =  -sin0m(t;  1Q01) 


-^hros(irh)  V  C°S  [C2fe—  1  )?rt/2T] 
*  (  )  hi  (2fc-l)2-4/i2 


and 


cos0m(t;  0110)  =  cos0m(t;  1Q01) 


(A-15a) 


sin  (irk) 
irh 


1-2 

1=1 


cos  (kirt/T) 

e~h2 


(A-15b) 


Assuming  the  filter  passes  only  harmonics  of  these  Fourier  series  up  to  /=  1/T,  we  find 
that  the  signal  quadrature  components  for  these  patterns  are 


u(t;  0U0)  =  — u(t;  1001) 

=  -^cos(»k)|i-L,|H„(^r)|cos[.(t-ti)/2T] 

+  |H„($,)|cos  [3»(t-t4)/2T]|  ( A-16a) 

=  -c4cos[ir(t-trf)/2T]  -  c5cos[3ir(t-td)/2T]  (A-16b) 


and  v(t;  0110)  =  v(t;  1Q01) 

-  ^ilH4)lcos[2,(t-tJ)/T)l 
=  c6  -  c7cos[jr(t-td)/T]  -  c8cos[2jr(t-td)/T]. 

At  the  sampling  time  t  =  td  +  T,  these  functions  give  rise  to  the  ISI  parameters 


(A-16c) 

(A-16d) 


I 
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A^(td  +  T;  0110)  =  -A^(td  +  T;  lflOl) 

= 

a2(td  +  T;  01.10)  =  a2(td  +  T;  1Q01)  =  (c6  +  c7  —  cg)2 

(A-17a) 

(A-17b) 

(A-17c) 

and 

a2(td;  0U0)  =  a2(td;  lflOl) 

(A-17d) 

=  (c4  +  c5)2  +  (c6— c7— c8)2 

(A-17d) 

in  which  the  constants  are  defined  as 

c<  -  T(i8-4V)c°s(ir'‘)|Ho(A')i 

(A-18a) 

Cs  -  7r(9?_/4^2)cos(?rM  Ih°(4^)| 

(A-18b) 

c6  4  sin('k>  =  sine (h) 

(A-18c) 

c;-iSc«'Ih»<2t)I 

(A-18d) 

and 

(A-lSe) 

Figure  A-2  illustrates  the  various  phase-related  quantities  for  pattern  0110. 

Patterns  1.100,  0011.  Since  pattern  1100  and  its  negative,  pattern  001 1.  are  time- 
shifted  versions  of  patterns  0110  and  1Q01,  we  can  immediately  write  that 


u(t;  1100)  =  — u(t;  0Q11)  (A-19a) 

=  u(t  +  T;  0110) 

=  c4  sin  [7r(t  —  td)/2T]  —  c5sin  [3?r(t — td)/2T)  (A-19b) 

and 

v(t;  1100)  =  v(t;  Oflll)  (A-19c) 

=  v(t  +  T;  0110) 

=  c6  +  C7  cos  [ir(t— td)/T]  -  cgcos[2jr(t-td)/T].  (A-19d) 

The  corresponding  ISI  parameters  at  the  sampling  time  axe 

A*(td  +  T;  1100)  =  — A^(td+T;  Oflll)  (A-20a) 

=  (A -20b) 

a2(td  +  T;  1100)  =  a2(td+T;  Oflll)  (A-20c) 

=  (c4+c5)2  +  (c6— c7— c8)2  (A-20d) 


I 
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and  a2(td;  0110)  =  a2(td;  1Q01)  (A-20d) 

=  (c6+c7-c8)2  (A-20d) 

Example  values  of  the  ISI  pattern  waveform  parameters  are  given  in  Table  A-l.  A 
composite  plot  of  all  the  differential  phase  waveforms  A<t>(t)  for  the  eight  patterns  is 
presented  in  Figure  A-3  for  example  values  of  the  input  parameters  h  and  W/f.T. 


TABLE  A-l 

ISI  PATTERN  PARAMETER  VALUES 

Conditions:  h  —  2/dT  =  0.7,  D  =  W/fT  =  1.0,  Gaussian  IF  filter 
CNR  =  carrier- to-noise  power  ratio 


Data  Patterns 

A<3(rad.) 

U/CNR 

V/CNR 

W/CNR 

111,  000 

±2.1991 

0.6806 

0.0000 

0.6806 

Oil,  100 

±1.710S 

0.7784 

-0.0997 

0.7719 

010,  101 

±1.2239 

0.8696 

0.0000 

0.8696 

110,  001 

±1.7108 

0.7784 

0.0997 

0.7719 

110 


2.  00 


3.  00 

bit  time  (t/T) 


4.  00 


5.  00 


111 - 

Oil 

110 - 

010 

101 - 

“IIFFERE 


6.  00 


111, 
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APPENDIX  B 

PROPERTIES  OF  THE  TRUNCATED  GAUSSIAN  DISTRIBUTION 
B.l  DEFINITION 

Let  a  truncated  Gaussian  probability  density  function  (PDF)  be  defined  as 


Prci1)  — 


Kt-«, 


lx  <  a 


otherwise. 


(B-l) 


In  the  text,  this  PDF  is  used  as  a  point  of  departure  for  development  of  an  approxima¬ 
tion  to  the  modulo-2^  differential  phase  PDF  pv(x).  with  a~i r. 


B.2  NORMALIZATION  FACTOR 

The  integral  of  the  PDF  (B-l)  must  equal  unity.  Therefore  we  have 


■a  ra 

1  =  dzprG(x)  =  2K  dxe 
—  a  J  0 


-*72* 


/•umar  2  a 

=  2K<r042  due"'  ,  umax  =  a/<70>[2 

o 

=  K<T0iJ5rerf(a  /*o42). 


(B-2a  ] 


Note  that 


I\  1  =  (70-[2lr  erf  (a/  <t0-^2  ). 


lim  K  = 

2  a 


(B-2b) 


(B-2c) 


B.2  VARIANCE 


Since  the  mean  value  is  zero,  the  variance  of  the  distribution  whose  PDF  is  given 
by  (B-l)  is  calculated  as 

c\G  =  2K  H  dxxV*272^ 

J  0 


ffl2/2<7o 

=  2>[2K<7q  du-fue-* 

J  o 

=  2'J2K<7q t(§;  a2/2<To) 

a  P(j:  a2/2 *1) 

<T°  erf(  a/cr0^ 2)  ’ 


(B-3a) 


(B-3b) 
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in  which  i{u\  x)  is  the  incomplete  gamma  function  and  [17] 
p/  .  A 

P(r-  THT 

=  ivrf) ,F,|li  "+1;  *> 


(B-4a) 

(B-4b) 


is  the  cumulative  probability  distribution  function  (CDF)  for  a  chi-squared  random 
variable  with  2u  degrees  of  freedom  [17,  chapter  26].  In  (B-4b),  the  ^  function  is  the 
confluent  hypergeometric  function.  Now,  since  [17] 

P(f:  *)  =  P(J;  ^  (B~5a) 

=  erf(  Jx )  -  (B-Sb) 

it  follows  that 

P(f:  nyiaD^eiHa/cr^)  -  ^e~a7/2(Tl  (B-5c) 

Substituting  this  result  in  (B-3b)  gives 


JL  J2  e-a/2<T° 


i  uo  _ 

Tc  ~  °  ■ 


An  empirical  inversion  of  this  formula  which  we  have  developed  for  a=r  is 


tTo  —  a  TQ' 

=  exp  | 


[  1  —  Inc  2tg  j  ln{  t 2  /  3 )] c 


< ri-r-  <  1 


°  TG>  P 


(B-6) 

(B-7a) 
( B  -  7b ) 


An  alternate  expression  to  (B-6j  is  obtained  by  substituting  (B-4b)  into  (B-3)  twice, 
making  use  of  the  fact  that  erf(-Jx)  =  P(|:  x)  in  the  denominator,  which  results  in 


0 

—  <r 

>F,(1;  1:  0 

C  =  — , 

“  3 

iF,(l;  f;  C )’ 

2  trV 

_  a? 

l+lc+i-U'+l 

■Mc3+-  ■ 

~  3 

l  +  K+HC2+§ 

■Hc3+-~ 

_  a? 
~  3 

{1 

~  i5<l  +  5TsC2  + 

16  ,-3  1184  (■*  , 

4725'  779625'  ‘r 

«  £ 
3 

{1 

~  +  5T5^2} • 

Solving  for  <rl  gives  the  approximation 


2  ^  2a‘( i  i  10/ 1  . 


(B-Sai 


(B-8b) 


(B-9) 
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3er2 

this  approximation  is  satisfactory  for  (?2rc>2  when  g=tt,  or  when  — For 
smaller  values  of  <?\G~  we  have  the  empirical  relations  given  in  (B-7). 

B.3  CHARACTERISTIC  FUNCTION 

The  characteristic  function  for  the  truncated  Gaussian  distribution  is  found  from 


CTG{ v)  =  k|“  dieil/re_r/2<7° 

erff  -  erf/  -  -A _ lifgA 

T2  )  ertj  g04~2  42  J 


2erf(a/<r042) 


Re 


{erff  _fi _ Ji^sY 

r  VoNr2  >T2  Jj 

erf(a/<r0-j2) 


(B-lOa) 

(B-lOb) 

(B-lOc) 


This  exact  form  is  not  suitable  for  computation.  Alternate  forms  for  the  characteristic 
function  can  be  developed,  however.  The  first  one  makes  use  of  the  symmetry  of  the 
integral  in  (B-lOa)  to  write 

CTG(u)  =  2I\ 


‘dze^ 


cos  ux 


(B-10a) 


=  2Kvr4K:  rv., 

^  (2n)! 


n  =  0 


iz  z?"  e-r!,2al. 


using  a  power  series  for  the  cosine  function.  Next,  making  use  of  the  fact  that  (2n)! 
equals  22rT(n-f  * )  n'/TH .  and  transforming  the  integral  by  y=x7/2crl  results  in 


craM  - 


■q2/2(t  o 


_  y  prr  xy'  (  —  \is2<rl)n  y(n+};  a? /2al) 

-  K*»'I2’.L  s - r(B+j) 


/*  e-y 


n  =  0 


=  P(»+J;  o?/2 »3). 


n  =  0 


(B-lOb) 

(B-lOc) 


Another  development  from  (B-lOa)  uses  integration  by  parts  to  obtain 


Ctg(u)  ~  2I\ 


sin  i/x~z  /2(To 


fl  +  2K  fQ 
0  ^<t2Jo 


,  -X72G  . 

die  smu 
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=  2Kae~a‘/2a'°s-^-^ 


au 


+  L  p(n+3.  a2/2 ai), 

n=0  n 


(B-lla) 


after  undergoing  a  development  similar  to  that  used  to  obtain  (B-lOc).  After  noting 
that 

2  2i  +  l 

r  (B-llb) 


P(n+f;  x2)  =  P(|;  x2)  -e"r  g  <  P(f;  z2), 


and  substituting  this  bound  into  (B-lla),  we  find  the  bound 
CTG(v)  <  2K ae 


-a2/2^o  • 


sm  at/ 
ai/ 


oc  ( _ 1 , .  2  2  \  n 


n  =  0 


n: 


=  2Kae"fl2/2ff“§ig^  +  K<r0^P(§;  a2 /2^)e'^ °i/2 


2  _2  , 


-v‘vl/2 


where 


=  t  +  (1-0-e 


f  6  O  J|e-,,/'-‘,J(erf(^or2)]-'. 


(B-llc) 

(B-lld) 


This  second,  mixture  expression  for  the  truncated  Gaussian  CHF  is  preferred  because, 
using  the  bound  as  an  approximation,  it  is  readily  discerned  from  it  that  the  CHF  has 
the  form 

CTG(v)  =s  {  ■  Cr(v‘,  a)  +  (1  -  0  •  o'o),  (B-12) 


in  which  Cv(u\  a)  is  the  CHF  for  a  random  variable  uniformly  distributed  on  (-a.  a) 
and  CG(i'\  al)  is  the  CHF  for  a  zero-mean  Gaussian  random  variable  with  variance  <72. 
Further,  the  mixture  parameter  £  is  seen  to  approach  the  value  {  =  1  as  co-^^c.  that  is. 
as  the  SNR-+0,  and  approaches  zero  as  SNR-»oc. 
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APPENDIX  C 

FUNCTIONS  RELATED  TO  THE  GAUSSIAN  FILTER 

In  the  analysis  of  CPFSK  systems,  it  is  con  ton  to  assume  an  IF  bandpass  filter 
with  noise  bandwidth  W,f  and  a  “Gaussian-shaped  spectral  transfer  function: 

H  lr(f)  =  H„(/-/,f)  +  H  „(/+/, r)  (C-la) 

where 

H  (C-lb) 


C.l  FILTER  IMPULSE  RESPONSE 

The  equations  above  imply  that  the  impulse  response  of  the  bandpass  filter  is 
modelled  by 

hJF(t)  =  2h0(t)cos(2ir/;ft).  (C-2a) 

with  the  “equivalent  lowpass  filter”  impulse  response  being 


MO  =  P 

J  -  oc 


=  WIFf2e~2^'jF{t-idr 


(C-2b) 


where 


I  -(t-td)2/  2c\ 
'  ^ 

a,  £  (2W ]F^)~\ 


(C-2c) 


(C-2d) 


This  lowpass  filter  response  is  nonrealizable  in  that  it  is  nonzero  for  t  <  0.  However,  if 
the  filter  phase  delay  td  is  sufficiently  long,  the  integral  over  positive  values  of  the  argu¬ 
ment  is  very  close  to  unity,  indicating  that  the  filter  response  for  t  <0  can  be  neglected: 


where 


[  dtho(t)  =  QG(-2W;ft1,'f5r)  ssl,  W;rtd  >  1 

Jo 

=  A1>V/’ 


(C-3a) 


(C-3b) 


=  5  l1  ~  erfl 


(C-3c) 


is  the  complementary  cumulative  probability  distribution  function  for  a  Gaussian  ran¬ 
dom  variable  and  erf(-)  is  the  error  function. 
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C.2  RESPONSE  TO  NARROWBAND  WAVEFORMS 

The  response  y(t)  of  the  bandpass  filter  to  a  general  narrowband  waveform  x(t) 
with  in-phase  and  quadrate  components  I(t)  and  Q(t), 

x(t)  =  I(t)  cos u jpt  -  Q(t)sin uJFt  (C-4) 

using  uJF  =2 *fIF,  can  be  developed  as  follows: 

ft 

y(t)  =  du[I(u)cosw;/-u  —  Q(u)sinw;fu]  •2h0(t-u)cosw/f(t— u) 

J  —  OC 
ft 

=  j  du  I(u)h0(t— u)  |cosw/ft  +  cos uIF (2u— t)j 

du  Q(u)h0(t-u)  {sin  w/f.t  +  sinw/f  (2u-t)| 

•t  rt 

du  I(u)h0(t-u)  -sinuijFt  du  Q(u) h0(t-u)  (C-5a) 

—  OC  J  —  OC 


cos~7/.t 


=  u(t)cos^/ft  - v(t)sinu>;/-t  (C-5b) 

with  u(t)  and  v(t)  being  the  in-phase  and  quadrature  components  of  the  filter  output, 
respectively. 

The  same  result  can  be  obtained  using  the  ‘‘complex  envelope"  notation  to  write 


and 


w’ith  the  understanding  that  the  effect  of  the  bandpass  filter  in  terms  of  its  equivalent 
lowpass  filter  can  be  expressed  by  the  convolution 


x(t)  =  Re{[I(t)+jQ(t)]<!,“'ft} 

=  Re{X(.)ej"'r‘} 

(C-6a) 

y(t)  =  Re|[u(t)+jv(t)]e“'rl} 

=  Re{Y(t)e“"',}> 

(C-6b) 

Y(t)  =  f  du  X(u)h0(t— u). 

J  — OC 


(C-7) 


C.3  INTEGRALS  RELATED  TO  FILTERING  OF  THE  TRANSIENT  PHASE 

The  following  integrals  are  used  in  Section  3  of  the  text  in  connection  with  the 
effect  of  a  Gaussian  IF  bandpass  filter  on  the  transient  phase  difference  due  to  imperfect 
synchronism  between  transmitter  and  receiver  hopping  synthesizers.  The  integrals  pro¬ 
ceed  from  the  desire  to  calculate 

3(tj.  t2;  is,  <r,  y)  =  (<t42?t)-1  [  dzcos2iri/x  e~^x-7^  .  (C-8) 

Jti 

The  solution  begins  with  a  simple  translation  to  obtain 
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ft2-7  -r2/2<72 

<r 5( t ! ,  t2;  v,  <r,  7)  =  cos2ff»/7  dzcos27ri/x  e 

JM- 7 


ft2-7  -x2/2<r 

—  sin  2*1/7  dxsin2*i/xe  ' 

Jti-7 

=  COs2*I/7  •  <7^2*  3c(tj  —  7,  t2  — 7!  <7) 


.2/0-2 


—  sin2iri/7  •  (t^2tt S* ( t x  —  7,  t2  —  7;  t/,  <7). 
The  first  integral  defined  in  (C-9)  is 

rB 

( 

A 


3c(A ,  B;  v,  <7 )  =  ( <r ) _  1  [  dxcos2*i/x  e_r  ^2<7 


-  £  “*?*  J>* 


2„  -*2/2<r2 


n  =  0 
2t,2^2*n 


=  i  £  (~g^5')'''  {l(n+i;  BVlc')-,(n+\:  AV 2S)} 


n  =  0 


!T(n+i) 


(C-9) 


(C-lOa) 


(C-lOb) 


in  which  7(0;  x)  is  the  incomplete  gamma  function,  a  series  expansion  of  the  cosine 
function  was  used,  and  the  fact  that  (2 re)!  =  r(n-f|)  n!  22n/V*  was  also  used.  Now. 
recognizing  that  [17] 


7(w+^:  *2)  7(5:  a2) 

r(n+J)  r(J) 


erfx, 


(C-ll) 


we  replace  the  gamma  functions  in  (C-lOb)  with  their  values  for  n  =  0  to  obtain  the 
approximation 

3 C(A.  B:  v,  a)  «  e'2(7ri/a)2  •  |{erf  (B/vfy-erf (v4/<7>T2)}.  (C-12a) 

In  the  text,  we  have  v  —  kfAe ,  harmonics  of  the  fundamental  frequency  of  the  transient 
phase  difference  waveform,  and  <r  =  <rh  is  given  by  (C-2d).  Therefore,  since  the  factor 
involving  the  error  functions  in  (C-12)  is  a  most  unity,  the  value  of  3C  is  governed  by 
the  exponential  factor,  which  equals 


e-2(™<7)2=  e-7T(^#/W/f)2/2 


|Ho(*/a,)|. 


(C-12b) 


The  second  integral  defined  in  (C-9)  is 

3,(/4,  £;  v,  a)  =  («7>[Jjr)"1[  dxsin2*t /xe"1^ 

JA 


(C-13a) 
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-i  ^  (-i)"(2™)2n+i  rB,,__2n+i„-*72<T2 


-  HS)"  E  ]***“'< 

__  .  .,  2  2  2 \n  . 


iio  (2»*>) 

.rnf  (-2 fW)“ 


^  z_  „,r(„jf)  {r(»+l;  B>/2»’)-T(n+l;  j42/2<t2)}  (C— 13b) 


using  the  same  analytical  techniques  as  were  used  to  obtain  (C-lOb).  Recognizing  that 
the  incomplete  gamma  function  for  integer  first  arguments  can  be  expressed  by  [17] 


7(n-fl;  x2) 
r(n+l) 


(C— 14) 


we  replace  the  gamma  functions  in  (C-13b)  with  their  values  for  n  =  0  to  obtain  the 
approximation 


3,(A,  B;  i/,  a )  ss 


ttvct 

42 


(-2ir2i/2<72)n 

"W 1 


(C- 15) 


in  wrhich  the  series  is  the  confluent  hypergeometric  function: 

oc  (  o_2  2  2\n 

VO)  2  I?  /I.  3.  o_2  2_2\ 

r(»4) 

=  t,  e-W1F,(ii  i;2>W) 


1  i2  for  i/<r »  1. 


2(  Jri/<ry 

Thus  we  have  the  asymptotic  approximation 


1  (a-AW  -52/2<t2- 


(C-16a) 

(C-16b) 

(C-16c) 


(C— 17) 


Clearly,  the  factor  involving  the  exponentials  can  at  most  equal  ±1.  Therefore,  the  mag¬ 
nitude  of  (C-17)  is  controlled  by  the  first  factor.  Substituting  for  v  and  er,  we  find  that 
this  factor  equals 


1  _  1 
2ttv<t  l2w  2jr42(fc/aS/W  jf) 


(C-18) 
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